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FOREWORD 


The  research  described  here  has  been  a  group  effort  but  S.  A.  Sheffield 

si 

and  K.  M.  Ogllvle  have  been  the  major  contributors  through  their  dissertation 
efforts  which  are  summarized  In  Sections  IV  and  V.  Their  dissertations*  which 
are  available  from  the  Washington  State  University  Library*  contain  much  more 
detail  than  Is  given  here,  and  they  should  be  consulted  by  persons  planning 
research  In  this  field.  The  format  Is  not  uniform  throughout.  To  make  It 
so  would  have  Introduced  more  errors  than  I  care  to  contemplate.  However, 

It  Is  uniform  within  sections,  and  references  are  listed  at  the  end  of  each  . 
section. 

The  report  contains  a  good  deal  of  background  Information  on  shock-  • 

Induced  reactions.  This  may  be  useful  to  those  unfamiliar  with  the  field. 

Others  will  want  to  concentrate  on  Sections  IV  and  V.  The  report  was  assembled 
In  haste  and  flows  less  smoothly  than  I  would  like,  but  I  hope  It  will  serve  a 
useful  purpose. 

I  had  Intended  to  give  a  searching  examination  of  the  various  physical 
mechanisms  of  shock-molecule  Interactions  In  Section  VI,  but  It  turned  out  to 
be  too  time-consuming  and  to  speculative,  so  It  will  wait  until  another  time. 

The  report  was  assembled  by  Ns  Debra  Magglora,  who  performed  a  minor 
miracle  In  reordering  figure  numbers,  references,  etc.  Errors  remaining  are. 

I'm  afraid,  mine  rather  than  hers.  Also,  I  am  grateful  to  suggestions  and  comments 
made  by  P.  S.  P.  Wei  of  the  Boeing  Aerospace  Company. 


G.  E.  Duvall 
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I.  INTRODUCTION 
1.1  Fast  Reactions 

An  obvious  fact,  often  Ignored,  Is  that  chemical  reactions  and  continuum 
dynamics  are  essentially  and  unavoidably  coupled  In  principle.  Only  In  the 
limits  of  weak  waves,  very  slow  reactions  or  very  small  samples  are  they 
uncoupled. 

The  constitutive  relations  for  a  chemically  reacting  substance  consist 
of  an  equation  of  state  for  the  unreacted  substance,  an  equation  of  state  for 
the  reaction  products,  and  a  rate  equation  which  couples  mechanical  and  chemical 
parameters.  For  example 


Cp2p  +  Ja 

(la) 

Ap  +  Ba 

(lb) 

where  Cp  is  the  frozen  sound  speed  and  the  coefficient  J  depends  on  the 

rates  of  change  of  energy  and  volume  with  a.  J  may  be  positive  or  negative 

1  2 

depending  on  the  heat  of  reaction  and  volume  of  the  reaction  products.  * 

For  one  dimensional  linearized  plane  flow  Eq.  (la)  yields 


The  reaction  acts  as  a  source  term  for  the  pressure  waves.  A  more  complete 
formulation  would  show  that  these  waves.  In  turn,  heat  the  medium  and  this 
heat  may  cause  the  reaction  to  spread. 
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The  reaction  may  be  called  "fast"  If  the  characteristic  time  for 
completion  Is  less  than  the  largest  dimension  of  the  system  divided  by  the 
frozen  sound  speed.  The  strengths  of  waves  normally  required  to  produce 
reactions  are  large  enough  that  non-linear  effects  come  Into  play,  and  In 
such  cases  shock  waves  are  formed.  Thus  we  come  to  the  subject  of  this 
report  on  shock-induced  chemical  reactions.  It  Is  evident  from  this  brief 
discussion  that  there  Is  adequate  theoretical  reason  to  Investigate  the 
dynamics  of  chemical  reactions  In  material  media  even  without  practical 
considerations.  The  latter,  however,  exist  In  abundance.  In  fact,  current 
Interest  In  the  mechanisms  of  shock-induced  chemical  reactions  In  condensed 
matter  stems  from  two  practical  problems:  detonation  processes  In  high 
explosives  and  remarkable  chemical  changes  reported  by  Soviet  scientists 
In  materials  recovered  after  having  been  shocked  to  high  pressure. 

1.2  Explosions— A  Relation  between  Shocks  and  Chemical  Reactions 

Alfred  Nobel,  In  1863,  Introduced  the  detonator  to  Initiate  nitro¬ 
glycerine  and  showed  "the  first  Insight  Into  the  Important  role  of  a  strong 
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shock,  as  distinct  from  heating.  In  the  process  of  Initiation."  This  was 
the  beginning  of  the  modem  era  of  using  high  explosives  In  rock  breaking 
applications. 

During  this  same  period,  the  theory  of  shocks  was  developing,  primarily 
applied  to  gases.  The  work  by  Ranklne4  in  1869  and  Hugonlot5  In  1887 
concerning  the  statement  of  the  conservation  laws  for  a  steady  shock  resulted 
In  what  Is  now  known  as  the  Ranklne-Hugonlot  jump  conditions.  These  relations 
form  the  cornerstone  for  steady  shock  theory  today. 

Chapman®  In  1899  and  Jouguet  ^  1905  used  the  theory  of  shocks  to 

develop  a  theory  to  explain  the  steady  detonation  problem  and  thus  began 


to  link.  In  theory,  shocks  to  the  chemical  reaction  which  takes  place  In  a 

detonating  high  explosive. 

a 

In  1942,  von  Neumann  developed  a  theory  that  Incorporated  shock 
theory  Into  the  reacting  part  of  a  steady  detonation  wave.  He  theorized  that 
there  are  a  number  of  partially  reacted  Hugonlots  between  the  unreacted  and 
the  fully  reacted  Hugonlots.  Since  the  unreacted  explosive  cannot  be  expected 
to  react  before  the  shock  wave  hits  It,  he  argued,  the  shock  carries  It  to 
a  very  high  pressure  and  then  as  reaction  occurs,  the  locus  of  states 
traverses  the  partially  reacted  Hugonlots  and  finally  ends  up  at  the  fully 
reacted  steady  state  corresponding  to  the  Chapman-Jouguet  state. 

The  von  Neumann  model  Is  shown  In  Fig.  1  where  OC  Is  the  Ranklne-Hugonlot 
(R-H)  curve  of  the  unreacted  explosive  and  ABD  Is  the  R-H  curve  of  the  reaction 
products.  The  dotted  line  OC  Is  the  Rayleigh  line  connecting  the  initial 
and  shock  states  In  the  explosives.  The  assumptions  of  the  model  are  that 
1)  the  Inert  material  Is  shocked  to  the  state  C  so  rapidly  that  no  reaction 
occurs,  and  11)  that  the  flow  Is  one-dimensional  and  steady  In  coordinates 
attached  to  the  shock  front.  The  consequences  of  these  assumptions  are  that 
the  locus  of  states  starts  at  0,  rises  rapidly  to  C,  and  decreases  more  slowly 
to  the  "Chapman-Jouguet"  (C-J) point  B  as  the  reaction  progresses  to  completion. 

The  resulting  profile  Is  shown  schematically  In  Fig.  2.  The  region  to  the 
left  of  the  C-J  plane  Is  assumed  to  be  Inert  and  Is  often  called  the  (G.  I.) 

Taylor  wave.  It  expands  according  to  the  flow  equations  of  an  Inert  gas. 

This  model  was  described  Independently  by  Zeldovlch^and  Doering.10  It  ! 
Is  usually  called  the  "ZND"  model.  The  ZND  theory  Is  quite  successful  In  \ 


describing  the  relations  among  mechanical  variables  In  steady  detonation  of  j 

condensed  explosives.  Elaborate,  semi -empirical  computer  codes  exist  for 

relating  chemical  formulations  to  detonation  pressures  and  velocities11  In  such  j 

cases.  It  Is  now  believed  that  both  fundamental  assumptions  are  violated  to  j 

some  degree  In  some  Instances:  l.e..  In  some  cases  there  may  be  reaction  In  j 

1 

the  shock  front,  and  flow  may  be  non-steady  and  to  some  degree  mul tl -dimensional . 
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C-J,  ZNO  THEORY 

NO  REACTION  IN  SHOCK  FRONT 
COMPLETE  REACTION  IN  "REACTION  ZONE 
D*U+C  AT  END  OF  REACTION  ZONE 

C-J  HYPOTHESIS 


C-J  plant  Deton.  is  nowa 

C-J  Deflagration 


Fig.  2  Schematic  Representation  of  Pressure  Profile  Produced  by 
a  steady  ZND  Detonation. 


Of  more  Importance  than  these  questions  are  those  relating  to  the  two 
principal  practical  problems  of  detonation:  will  a  given  chemical  substance 
detonate  reliably  In  a  particular  configuration,  and  how  sensitive  Is  It  to 
accidental  Initiation?  These  are  usually  formulated  as  the  "critical"  or 
"failure  diameter"  and  "Initiation"  problems.  The  latter  Is  sometimes  called 
the  "DOT"  problem:  Deflagration  to  Detonation  Transition. 

1.2.1  Critical  Diameter  for  Detonation 
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Theories  of  the  critical  diameter  have  been  given  by  H.  Jones,  Eyrlng 
et  a]..,1*  Wood  and  Kirkwood,1®  V.  Rozlng  and  Iu.  B.  Khariton,16  Dremln,17 
Chalken  and  Edwards,  and  Apln  and  Bobylev.  Less  encompassing  contributions 
have  been  made  by  many  others.  Khariton,  Jones  and  Dremln  address  the  problem 
of  cooling  of  the  reaction  zone  near  the  boundary  of  the  detonating  charge 
by  the  rarefaction  which  originates  where  detonation  front  and  charge  boundary 
Intersect.  Eyrlng  et  aj[.  and  Wood  and  Kirkwood  developed  the  hydrodynamics 
of  curved  detonation  fronts  In  some  detail  and  applied  their  results  to 
finite  cylinders  with  homogeneous  first  order  reactions.  Apln  and  Bobylev  and 
Chalken  and  Edwards  have  considered  the  problems  of  grain  burning  In  powered 
and  heterogeneous  explosives.  Each  of  these  theoretical  developments  has 
had  some  success  In  correlating  data  on  diameter-dependence  of  detonation 
velocity  In  soft*  substances.  None  of  them  has  been  used  successfully  to 
predict  failure  diameter  In  a  substance  for  which  steady  detonation  has  not 
been  demonstrated.  Moreover,  each  of  them  appears  to  fall  In  some  critical 
way  when  applied  to  a  variety  of  detonating  systems. 

For  example,  the  diameter  dependence  of  detonation  velocity  In  liquid 
TNT  and  nltromethane  Is  fitted  by  the  theory  of  Eyrlng  et.  Ah ,  but  a  variety 
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of  high-density  heterogeneous  explosives  show  unrelated  behavior.  The 

Wood  and  Kirkwood  theory  relates  detonation  velocity  to  curvature  of  the 

wave  front  near  the  charge  axis.  Curvature  Is  difficult  to  measure,  but 

In  at  least  two  cases  It  appears  to  fit  data  for  heterogeneous  explosives 
20  21 

passably  well.  *  Enlg  and  Petrone  found  that  Dremln's  quenching  theory 

explains  the  temperature  dependence  of  detonation  velocity  In  nltromethane 

but  Titov  et  aK ,  who  measured  critical  diameters  In  cast  TNT  at  temperatures 

of  20.4,  77.4,  and  291 aK,  found  that  changes  with  temperature  had  the  wrong 

23 

sign.  If  Arrenhlus  kinetics  apply,  whatever  theory  Is  used. 

18 

Chalken  and  Edwards  have  gone  a  bit  farther  than  others  Inasmuch  as 
they  have  explicitly  Introduced  two  different  reaction  processes  In  a 
simplistic  theory  which  produces  a  pleasing  correlation  of  critical  diameter 
measurements  among  various  Group  1  and  Group  II  explosives;  but  their  theory 
Is  Incapable  of  predicting  critical  diameters  from  fundamental  physical, 
chemical  and  geometric  properties. 

The  near-success  of  these  theories  of  detonation  based  on  ad  hoc  assumptions 
about  the  reaction  process  and  the  role  of  the  shock  In  Initiating  reactions  has 
encouraged  workers  In  the  field  to  believe  that  just  a  slight  improvement  In  the 
descriptive  parameters  may  be  all  that  Is  required  to  enable  us  to  understand 
non-ideal  detonation  processes.  But  results  experienced  to  date  do  not 
justify  such  optimism  and  suggest,  by  Implication,  that  what  Is  required  Is 
a  more  fundamental  look  at  shock  waves  and  the  chemical  reactions  they  produce. 
Condensed  explosives  have  some  characteristics  which  tend  to  disqualify  them 
as  suitable  subjects  for  studies  of  this  kind  In  the  laboretory.  The  reacting 
molecules  are  normally  complex,  so  the  reactions  will  be  complex.  Reaction 
zones  In  explosives  which  detonate  In  laboratory  sizes  are  smell,  so  the 
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resolution  required  for  studies  In  the  reaction  zone  Is  Inordinately  large. 

And  the  destructive  action  of  explosives  complicates  the  experiments,  both 
In  preparation  of  specimens  for  study  and  In  containment  of  explosions. 

1.2.2  Initiation  of  Detonation 

Accidents  In  the  use  of  explosives  were  common  In  early  years.  Their 
causes  were  cloaked  In  mystery  and  attempts  to  clarify  and  categorize 
explosives,  handling,  and  safety  problems  were  manifold.  One  result  of 
this  concern  was  development  of  a  battery  of  tests  for  rating  the  safety  of 
any  given  explosive  formulation:  drop  tests,  thermal  tests,  spark  tests, 
friction  tests,  etc.  Operation  of  the  tests  was  not  well  understood.  None- 
the  less,  a  manageable  situation  had  been  reached  by  World  War  II,  though 
accidents  had  not  been  eliminated.  During  World  War  II  there  was  strong 
Interest  In  the  explosive  process.  In  new  explosives,  and  In  using  old 
explosives  In  new  ways.  The  scientific  comnunlty  was  brought  Into  the  field 
to  an  unprecedented  extent,  and  more  detailed  formulations  of  the  problems 
of  Initiation  and  detonation  were  developed  than  ever  before.24”28 

It  was  realized  a  generation  ago  that  experiments  In  which  Initiation 
of  detonation  was  produced  by  single  passage  of  a  shock  wave  offered  potential 
for  understanding  the  Initiation  process  that  Is  absent  in  drop  tests.  Various 
versions  of  the  shock  test  have  been  used,28  but  the  flying  plate  test  Is  most 
attractive  because  of  Its  relative  simplicity  of  Interpretation,  based  on 
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assumptions  of  one-dimensional  plane  flow. 

In  the  flying  plate  experiment  a  uniform  shock  wave  of  duration,  t, 
determined  by  the  flyer  thickness  Is  generated.  Amplitude  of  the  shock  wave 
Is  determined  by  speed  and  equation  of  state  of  the  flyer  end  by  the  equation 
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of  state  of  the  explosive  sample.  Observations  which  can  be  made  are  limited. 

Mean  velocity  of  the  shock  wave  In  the  sample  can  be  determined  as  a  function 

of  sample  thickness;  recently,  pressures  recorded  by  gauges  Imbedded  In 
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the  sample  have  been  reported.  *  Photographic  observations  can  be  made, 
and  If  the  sample  Is  transparent,  these  may  give  detailed  Information  about 
the  speed  and  local  character  of  the  shock  and  the  Induced  reaction.  The 
principal  results  of  such  experiments  have  been  to  define  values  of  mechanical 
parameters  which  are  required  to  Initiate  detonation  and  to  verify  the 
speculation  that  the  explosive  reaction  may  start  at  flaws  or  Inhomogeneities, 
or  that  It  may  (apparently)  be  a  bulk  phenomenon  labelled  "thermal  explosion." 
As  an  example  of  the  mechanical  parameter  criterion.  It  Is  found  that,  for 
some  explosives,  detonation  Is  Initiated  when  P^t  exceeds  a  characteristic 
value.  For  other  explosives  the  exponent  of  P  may  differ  from  two. 

Detailed  models  of  the  evolution  of  the  Initiating  reaction  Into  a 
stable  detonation  wave  have  been  built  for  comparison  with  plate  or  other 
experiments,  and  these  meet  with  varying  degrees  of  success.  Difficulties 
are  probably  due  to  the  use  of  relatively  simple  models  for  the  reaction 
process  and  for  the  enormous  variations  In  chemical  and  physical  properties 
of  explosives. 

In  attempting  to  define  the  initiation  problem  In  a  way  which  points 
toward  possible  solutions.  It  Is  useful  to  divide  the  process  Into  several 
reasonably  distinct  sub-processes.  No  such  classification  Is  unique,  but 
the  following  division  provides  a  framework  for  considering  various  aspects 
of  the  total  problem: 

1.  Microscopic  physical  conditions  required  to  start 
the  chemical  reaction. 
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-  -  — - - 

If  Initiation  Is  mechanical  or  thermal,  stress,  time 
and  temperature  must  be  prescribed.  The  Initiation  may  be 
heterogeneous  or  homogeneous,  though  given  the  nature  of 
condensed  matter,  apparently  homogeneous  reactions  are 
simply  heterogeneous  on  a  small  scale.  Critical  dimensions 
must  also  be  Included  In  this  determination:  The  volume 
of  material  reacting  must  exceed  some  critical  size  In 
order  for  the  reaction  to  grow;  or  there  must  be  adjacent 
sites  which  can  combine  to  produce  an  overall  stable  volume. 

If  the  Initiation  Is  not  mechanical  or  thermal,  as,  e.g., 
photoinitiation,  there  may  be  other  physical  factors  which 
are  Important,  like  bond  resonance.  For  thermally  Induced 
reaction,  the  reaction  parameters  and  their  variations  with 
pressure  or  stress  must  be  determined.  For  mechanical 
Initiation,  the  creation  of  Inhomogeneous  sites  must  be 
described. 

2.  Mechanisms  by  which  the  Initiation  conditions  are  established 
In  any  particular  Initiation  process. 

For  thermal  and  mechanical  processes  which  can  be  treated 
as  homogeneous,  where  thermal  explosion  theories  are  believed 
to  apply,  this  Is  a  relatively  simple  problem.  The  available 
energy  Is  simply  distributed  over  the  heated  or  shocked  or 
compressed  volume  In  a  suitably  uniform  way.  For  heterogeneous 
Initiation  the  problem  Is  much  more  difficult.  Gamow^*  treated 
Initiation  of  granular  explosives  by  assuming  that  mechanical 
energy  of  compression  was  Irreversibly  converted  to  heat 
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through  friction  and  that  this  heat  was  distributed  through¬ 
out  the  explosive  grains  In  a  manner  defined  by  the  heat 
conduction  equations.  Other  mechanisms.  Including  jetting, 
fracture,  adiabatic  compression  have  been  invoked  by  various 
workers . 

Af anas’ ev  and  BoboTev^  have  carried  out  a  thorough  and 

careful  analysis  of  drop  tests  and  concluded  that  the  evidence 

Is  overwhelming  that  the  Initiation  process  Is  heterogeneous. 

They  have  ventured  to  suggest  that  the  mechanism  responsible 

for  the  heterogeneity  Is  shear-induced  plastic  flow  or  slip. 

However,  they  have  not  taken  Into  account  the  possibility  that 

the  rate  law  at  high  temperatures  and  pressures  Is  quite 

00 

different  from  that  at  low  temperatures  and  pressures. 

Other  observations,  as  well  as  theoretical  considerations, 

certainly  suggest  very  strongly  that  Initiation  Is  hetero- 
29 

geneous. 

3.  Processes  by  which  the  initial  reactions  spreads  from  a  point 
or  set  of  points  and  becomes  a  progressing  volume  reaction  with 
a  well-defined  front,  l.e.,  a  deflagration. 

In  the  stage  before  the  reaction  becomes  self-supporting, 
ft  must  be  supported  by  the  Initiating  agency.  In  the  case 
of  mechanical  or  thermal  initiation,  this  means  that  supporting 
tmperature  and  pressure  must  be  maintained  until  the  rate  of 
release  of  energy  from  the  reaction  equals  or  exceeds  losses 
to  the  surrounding,  unreacted  material.  The  processes  of  growth 
of  the  reaction  during  this  critical  period  depend  entirely 
on  detailed  properties  of  the  explosive.  If  initiation  occurs 
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at  a  few  points,  well -separated,  growth  of  the  Individual 
reaction  volume  Is  controlling.  It  may  proceed  along  planar 
flaws.  It  may  be  quasi spherl cal ,  or  It  may  have  a  strongly 
randan  component  In  a  mixture  of  binder  and  subcrltlcal  grains. 
Pressure  Increments  resulting  from  Incremental  reaction  may 
or  may  not  be  important  In  promoting  the  reaction,  depending 
on  both  chemical  and  physical  properties.  If  Initiation  points 
are  cTosely  spaced  so  that  the  reaction  Is  quasi -homogeneous, 
the  reacting  volume  may  be  sufficiently  large  that  when  the 
reaction  has  proceeded  beyond  a  critical  point.  It  may  change 
from  subcrltlcal  to  supercritical  nearly  dlscontlnuously.  Con¬ 
struction  of  a  variety  of  theoretical  models  for  various  geo¬ 
metries  and  properties  may  shed  useful  light  on  this  step. 

Processes  involved  In  the  transition  from  deflagration  to 
detonation. 

Bowden^® and  Bobolev  and  Dubovlk®*"®®  have  mapped 

the  progress  of  reactions  In  small  samples  Initiated  In  Impact 

machines.  They  found  the  reaction  front  accelerating  smoothly 

to  a  speed  of  a  few  hundred  meters  per  second,  then  changing 

almost  dlscontlnuously  to  high  or  low  order  detonation.  Details 

depend  on  the  explosive.  Experiments  In  which  detonation  is 

Initiated  by  plate  Impact  on  unconfined  explosives  appear  to 

show  analogous  behavior  as  Indicated,  for  example,  by  gage 

records  like  those  reported  by  Kennedy  and  Nunzlato^7  and  by 

3B 

Wackerle,  Johnson,  and  Hal leek.  It  seems  likely  that  the 
development  of  accidental  explosions,  typified  by  the  detonation 
of  rocket  motors,  follows  a  similar  course,  though  such  cases 

are  hardly  subject  to  careful  experimental  study. 
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It  should  be  evident  from  the  preceding  description  that  the  technical 
problems  of  initiation  are  enormous  and  that  no  simplistic  analysis  Is 
likely  be  to  effective  in  describing  them.  Microscopic  processes  are 
central  to  the  entire  problem.  These  consist  of  chemical  changes  in 
the  molecules  produced  by  the  shock,  defects  which  are  produced  in  the 
shock  front,  and  response  of  the  explosive  molecules  to  stress  and 
temperature  in  their  local  environment. 

The  deflagration  to  deformation  transition  is  particularly  Important 
in  thermal  initiation.  Some  notes  on  the  logical  framework  for  this 
process  are  given  In  Appendix  A. 

One  of  the  characteristics  of  explosives  which  makes  microscopic  model* 
ling  difficult  Is  that  they  are  usually  rather  large  molecules,  so  reaction 
paths  are  many  and  complex.  Another  difficulty  Is  In  the  violence  and  speed 
of  reaction,  which  places  severe  restraints  on  the  kinds  of  measurements  which 
can  be  made.  It  Is  a  fundamental  precept  of  research  that  If  the  problem 
one  wishes  to  understand  Is  too  complicated  to  understand,  then  a  simpler  prob* 
lem  embodying  some  aspects  of  the  real  problem  is  studied.  This  suggests 
that.  Instead  of  studying  explosives  with  large  molecules,  one  study  explosives 
with  simple  molecules;  or,  even  better,  small  molecules  which  decompose  under 
the  action  of  shock  waves,  but  not  violently.  A  practical  benefit  from  such 

study,  besides  the  Insight  Into  explosive  processes,  may  be  an  understanding 
of  chemical  processes  In  shock  which  will  lead  to  the  development  of  methods 
for  synthesizing  new  materials. 
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1.3  Alteration  of  the  Chemical  and  Physical  Constitution  of  Materials 
by  Shock 

For  at  least  sixty  years  Investigators  have  been  shocking  samples  with 

40 

explosives  or  by  Impacting  and  examining  them  for  shock-induced  changes. 

The  Intensity  of  this  kind  of  Investigation  has  Increased  since  World  War  II. 

In  the  Soviet  Union  It  has  become  highly  refined  and  systematized.  Several 

hundred  papers  on  the  shock  chemistry  and  synthesis  of  condensed  organic 

and  Inorganic  substances  have  appeared  In  the  last  twenty  years.  Shock 

compression  has  come  to  be  regarded  In  the  USSR  as  a  useful  and  practical 

technique  for  synthesis  of  materials  having  unusual  properties ,  and  research 

41 

Is  continuing  at  a  strong  pace. 

There  Is  substantial  evidence  to  Indicate  that  chemical  changes  wrought 

by  shocks  are  not  explicable  In  terms  of  pressure  and  temperature,  but  that 

the  shock  front  Itself  Is  Involved.  Graham  has  called  this  the  "catastrophic 

42 

shock"  effect,  and  has  adduced  considerable  evidence  to  support  Its  reality. 

In  the  postwar  period,  also,  there  have  been  a  substantial  number  of 
dynamic  measurements  In  liquids  which  attest  to  shock-induced  changes  In 
their  chemical  or  physical  states.  In  these  experiments,  properties  such 
as  opacity,  conductivity,  shock  velocity,  temperature,  pressure,  etc.,  were 
measured.  Anomalous  results  In  some  materials  have  led  experimenters  to 
postulate  the  occurence  of  various  types  of  transitions  and  reactions. 
Interpretation  of  such  experiments  Is  frequently  confused  because  of  similar 
effects  produced  by  phase  transitions  and  chemical  reactions.  A  sumnary  of 
experiments  Is  contained  In  Table  I. 

Walsh  and  Rice43  (1957)  showed  that  carbon  tetrachloride  (CCl^)  becomes 
Increasingly  opaque  In  a  shock,  beginning  at  60  to  80  kbar,  and  becomes 
completely  opaque  at  130  to  170  kbar.  Doran  and  Ahrens  (1963)  state 
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that  CC14  Is  a  good  conductor  at  170  kbar.  01ck  'J  (1964)  Indicates  CC14 
begins  to  conduct  at  70  kbar  and  Is  a  good  conductor  above  120  kbar. 

Mitchell  and  Keeler46  (1968)  measured  the  conductivity  of  CC14  from  65  to 
160  kbar  and  found  It  continuously  Increased  by  seven  orders  of  magnitude 
over  this  range.  Hamann  et  al.4^"49  (1959,  1960,  1966)  measured  the 
conductivities  of  several  liquids  In  150  to  170  kbar  shocks  and  found  that 
some  were  good  conductors.  The  conductivity  of  water  was  measured  at  many 
pressures  up  to  130  kbar  and  found  to  continually  Increase*,  the  reason  for 
the  changes  was  postulated  to  be  self-1 onlzatlon. 

Dick50,51,59  (1968,  1970)  showed  that  benzene  (CgHg),  carbon  disulfide 
(CSg) ,  CC14,  and  liquid  nitrogen  all  have  transitions  of  some  kind  at  133, 

62,  165,  and  135  kbar,  respectively.  Kormer  et  al. 62,63  (1968)  Indicate 
water  becomes  opaque  (which  they  Interpreted  as  evidence  of  freezing)  In  a 
double  shock,  and  that  temperature  measurements  of  CC14  show  It  has  an 
anomalous  behavior  between  200  and  300  kbar,  which  Is  attributed  to  a  solid- 
liquid  phase  transition. 

Yakusheva,  Yakushev,  and  Dremln64  (1971)  showed  that  CgHg,  CC14,  acetone 
(C3Hg0),  dlchloroethane  (CgHgClj),  and  nltromethane  (CH3N02)  became  opaque 
at  pressures  above  135,  80,  100,  120,  and  85  kbar,  respectively.  They 
attributed  this  opacity  to  decomposition  and  the  formation  of  free  carbon. 
Afanasenkov  et  al.66  (1974)  measured  temperatures  of  carbon  disulfide  and 
benzene  shocked  above  the  transition  and  found  them  to  be  higher  (by  500aK  In 
the  case  of  carbon  disulfide)  than  calculated  for  shock  heating  only.  They 
also  calculated  P-V  Hugonlots  for  postulated  products  of  reaction  and  showed 
the  P-V  data  for  carbon  disulfide  and  benzene  above  the  transition  to  be  very 
close  to  calculated  values.  Yakushev,  Nabatov,  and  Yakusheva66  (1974)  measured 
shock  velocity,  particle  velocity,  light  absorption  (opacity),  electrical 
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conductivity,  dielectric  permeability,  and  shock  polarization  and  did 
recovery  experiments  on  acrylonitrile  (CjH^N).  They  observed  anomalous 
behavior  beginning  about  43  kbar  and  explained  It  by  saying  a  decomposition 
reaction  took  place.  The  recovery  experiments  supported  this  explanation. 
Dick57  (1975)  states  that  transitions  In  benzene  ring  materials  depend  on 
the  number  of  carbon  double  bonds,  with  lower  transition  pressures  associated 
with  larger  numbers  of  double  bonds. 

The  anomalous  behavior  In  liquids  has  been  explained  chiefly  by  three 
possible  mechanisms:  phase  transitions,  polymerizations,  and  chemical 
reactions.  It  Is  possible  that  all  of  these  have  been  Involved  In  some  of 
the  experiments  while  only  one  may  have  been  Involved  In  others,  but 
separating  the  effects  has  not  yet  been  possible.  Many  experimenters 

52 

have  looked  for  liquid-solid  phase  transitions,  but  only  Kormer  et  al . 
claim  to  have  observed  It.  The  idea  has  been  expressed  several  times 
(see  Ref.  43  for  example)  that  too  much  ordering  Is  required  for  this  type 
of  transition  to  take  place  In  the  microsecond  time  scales  of  a  shock 
experiment.  Chemical  reactions,  primarily  of  the  decomposition  or  pyrolysis 
type,  explain  many  of  the  transitions, 54-56 ,58  but  s^nce  uttle  fs  known  about 
this  process,  the  explanations  are  not  yet  completely  substantiated. 

However,  these  recent  papers5*"®6  Indicate  that  chemical  reactions  In  other 
than  high  explosives  can  take  place  In  a  shock  environment.  They  also  show 
that  reactions  yield  large  enough  changes  In  measureable  parameters  that 
conventional  shock  measurement  techniques  will  yield  accurate  Information 
about  properties  like  reaction  time  If  the  experiments  are  carefully  done. 
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II.  PROBLEMS  OF  SHOCK  EXPERIMENTATION 

The  results  of  experiments  to  be  described  in  this  report  are  probably 
of  greater  interest  to  those  outside  the  "shock  wave  community"  than  to  its 
members.  A  great  deal  of  work  has  been  done  on  chemistry  stimulated  by  shock 
waves  in  gases.  Very  little  has  been  done  on  condensed  matter.  The  condition 
of  shock  experiments  in  solids  and  liquids  are  quite  different  from  those  in 
gases  in  various  particulars.  Pressures  are  very  much  higher  and  temperatures 
are  much  lower.  Duration  of  the  shocked  state  is  limited  to  nanoseconds  or 
microseconds  compared  to  milliseconds  for  gases.  Very  Important  Is  the  fact 
that  there  is  no  "rigid  wall"  for  lateral  confinement  of  condensed  matter 
shocks.  This  implies  that  transverse  probing  of  the  state  behind  the  shock 
is  not  usually  feasible— existence  of  the  unique  state  of  uniaxial  strain 
which  makes  the  one-dimensional,  plane,  time-dependent  shock  easy  to  analyze 
Is  limited  by  the  Inexorable  Inward  progress  of  lateral  rarefactlve  waves 
generated  at  free  surfaces  by  passage  of  the  shock.  This  dictates  the  shape 
of  the  sample.  It  must  be  a  disk  of  thickness  small  compared  to  Its  lateral 
dimensions  and  shock  fronts  must  be  parallel  to  the  face  of  the  disk. 

Diagnostic  probes  must  sample  changes  in  the  direction  of  shock  propagation 
and  only  In  the  region  where  the  uniaxialstrain  state  is  retained  for  the 
duration  of  the  experiment. 

Probes  are  generally  of  two  types:  those  which  sample  time-varying 
conditions  In  a  plane  parallel  to  the  shock  front,  and  those  which  sample  a 
finite  volume  extending  a  significant  distance  In  the  direction  of  propagation. 
The  former  are  exemplified  by  free  surface  measurements,  quartz  gages, 
electromagnetic  velocity  (EMV)  gages  and  manganln  resistance  gages.  The 
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latter  are  required  for  resistivity,  piezoelectricity,  and  optical  transmission 
measurements.  Examples  of  EMV  gage  and  optical  transmission  measurements  are 
described  in  Sections  IV  and  V,  respectively. 

Generation  of  shock  waves  is  by  the  impact  of  explosion  or  of  a  projectile 
against  a  sample.  The  impact  must  be  planar  within  a  few  seconds  of  arc, 
and  diagnostic  Instruments  must  be  activated  within  the  few  hundreds  of 
nanoseconds  when  the  sample  is  in  the  shocked  state.  Picosecond  timing  has 
become  commonplace  in  various  kinds  of  pulsed  laser  experiments,  but  is 
yet  to  be  achieved  in  shock  experiments.  The  reason  is  that  "mechanical 
jitter"  due  to  uncertainty  In  positions  and  dimensions  of  physical  components 
normally  ranges  from  about  50  to  250  nsecs,  depending  on  the  particular 
experiment. 

The  experiments  are,  of  course,  single  shot  affairs.  A  few  microseconds 
after  impact,  the  event  and  the  experiment  are  over.  A  very  good  experimentalist 
with  a  good  machinist  and  some  additional  technical  help,  can  build  and  fire 
a  difficult  shot  about  every  two  weeks,  provided  the  design  does  not  change 
and  there  are  moderately  small  parametric  changes.  This  does  not  include 
analysis  time,  which  can  easily  add  another  week  and  a  half  per  shot,  or 
experiment  design  time. 

The  consequences  of  this  slow  rate  are  that  each  experiment  must  be 
designed  to  give  as  much  Information  as  possible  and  that  design  must  be 
careful  and  thorough,  if  anything  at  all  is  to  be  achieved. 

Experimentation  Is  further  complicated  by  the  generating  mechanism. 

Whether  It  Is  explosive  or  projectile  Impact,  the  presence  of  the  generator 
severely  restricts  access  to  the  Impact  face  of  the  cell.  This  has  little 
effect  on  In-plane  gage  measurements.  It  presents  serious  problems  for 
volume-dependent  measurements.  In  some  cases  It  Is  possible  to  get  around 
these  restrictions,  as  can  be  seen  in  Section  V. 
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There  Is  a  further  problem  In  doing  new  kinds  of  shock  experiments. 

There  Is  considerable  lore  and  technical  art  Involved  In  shock  experimentation. 
Uhen  a  new  problem  arises,  like  spectroscopy,  It  takes  some  time  and  creative 
skill  to  develop  the  necessary  arts.  To  do  spectroscopy,  for  example, 
suitable  light  sources  must  be  developed,  windows  which  remain  transparent 
under  shock  must  be  found,  luminescence  of  structural  materials  must  be 
examined,  etc. 

Furthermore,  It  Is  necessary  to  develop  an  equation  of  state  for  each 
substance  or  mixture  studied.  In  experiments  where  shocks  are  created  by 
projectile  Impact,  the  principal  controllable  parameter  which  determines 
shock  strength  Is  projectile  velocity.  From  this  it  is  necessary  to  calculate 
pressure  and  temperature  In  the  sample  as  a  function  of  time.  This  is  done 
by  combining  the  equations  of  continuum  mechanics  with  equations  of  state 
for  each  significant  substance  in  the  experiment:  impactor,  front  of  cell, 
sample,  back  of  cell.  Equations  of  state  for  condensed  materials  are  not 
usually  well  known,  and  creating  one  from  available  data  Is  a  mixture  of  art 
and  science.  Since  pressure  In  a  shock  wave  can  be  obtained  directly  from 
the  jump  conditions,  pressures  calculated  from  these  synthesized  equations 
of  state  will  be  reasonably  accurate.  Temperatures,  however,  are  not 
normally  measured  in  shock  experiments,  so  there  Is  essentially  no  check 
on  temperature  values  obtained  from  the  synthetic  equations  of  state. 
Consequently,  the  calculated  temperature  Is  more  In  the  nature  of  an  ordering 
parameter  than  a  true  temperature.  This  Is  an  unfortunate  circumstance  which 
can  be  corrected  In  time*,  for  the  present.  It  can  only  be  noted. 

One  of  the  Important  data  sets  used  to  synthesize  an  equation  of  state 
Is  the  R-H  curve  for  ambient  Initial  conditions.  These  exist  for  many  materials 
When  they  do  not,  additional  experiments  are  required  to  produce  one. 
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III.  CHOICE  OF  SAMPLE  MATERIALS 


It  should  be  evident  from  the  comments  of  the  preceding  section  that 
Investigation  of  reactions  In  a  given  substance  or  mixture  of  substances  will 
require  a  major  Investment  of  time  and  funds.  This  means  that  careful 
thought  must  be  given  to  sample  selection. 

There  should  be  good  reason  to  believe  that  (1)  the  system  selected 
will  react  under  pressure  and  temperature  conditions  which  can  be  reached 
wlthavallable  experimental  facilities,  (11)  Its  study  will  lead  to  enlighten¬ 
ment  on  some  chemical  or  physical  principal,  (111)  the  scale  of  reaction, 
l.e.,  time  and  space,  lies  In  the  accessible  range  of  experimentation, 

(1v)  safety  hazards  can  be  handled  adequately.  In  our  case,  which  was  a 
pioneering  venture,  there  appeared  to  be  two  candidates,  based  on  study 
of  the  literature  cited  in  Table  I,  Section  1.3:  acrylonitrile  and  liquid  CS 
Table  I  shows  the  comparison  for  CS2  and  acrylonitrile.  The  data 
for  the  two  materials  are  about  the  same  except  in  the  molecular  structure 
and  density.  CS2  has  a  far  simpler  structure  and,  therefore,  the  simpler 
decomposition  reaction  mechanism. 


TABLE  :I  — Comparison  of  Carbon  Disulfide  and  Acrylonitrile 
as  Candidates  for  Further  Experiments 


Material 

Hugonlot 

Anomaly 

(kbar) 

Particle  Velocity 
at  which  Anomaly 
Occurs  (mm/usec) 

Density 

g/cm3 

Molecular 

Structure 

High 

Explosive 

Carbon 

Disulfide 

(cs2) 

62 

1.35 

1.26 

S«C«S 

no 

Acrylo¬ 

nitrile 

<c3h3n) 

43 

1.40 

0.806 

H  R 

i  i 

H-C-C-CSN 

no 

1  2 
As  Indicated  previously,  Walsh  and  Rice  and  later  Cook  and  Rogers 

did  Hugonlot  measurements  on  CSg.  Neither  of  them  found  any  evidence  of 

reaction  although  one  point  measured  by  Rice  and  Walsh  was  apparently  In 

the  reacted  region. 

3  4 

R.  Dick  *  did  a  large  group  of  shots  on  CSg  using  explosive  plane 
wave  generators  to  provide  the  shocks  and  found  the  anomalous  behavior  to 
begin  at  62  kbar  but  could  not  measure  a  two-wave  structure  (see  Fig.  3.1). 

He  also  could  not  say  definitely  what  type  of  reaction  was  occurring,  although 
he  seemed  to  favor  the  conversion  to  the  "bl ack ,  solid  CS^1  first  discovered  by 

c 

Bridgman  In  static  experiments. 

lysne®  made  two  low  pressure  shots  on  CS2  In  which  he  measured  the 
pressure  and  shock  velocity,  primarily  to  see  If  the  data  indicated  that  the 
shock  velocity-particle  velocity  Hugonlot  extrapolated  to  sound  speed.  He 
found  no  Indication  of  anomalous  behavior  In  the  0  to  11  kbar  pressure  range. 

Butcher  et  al.7  did  a  number  of  static  high  pressure  tests  at  elevated 
temperatures  on  CS2  to  produce  and  study  the  “black  solid  CSg”  Bridgman  had 
discovered  earlier.  From  their  data  and  Bridgman's  data,  they  constructed 
a  diagram  that  Is  reproduced  In  Fig.  3.2.  They  Indicate  the  “black  solid  CS2" 
Is  formed  only  In  a  relatively  small  region  as  Indicated  In  the  figure  and 
that  at  temperatures  and  pressures  above  this,  decomposition  occurs.  The 
regions  marked  decomposition  Is  where  the  shock  results  of  Dick  (with  the 

reaction  occurring)  are  located. 

a  7  3  4 

Afanasenkov  et  al .  combined  the  work  of  Butcher  et  al.  and  Dick  * 

to  do  some  Interesting  calculations.  They  assumed  a  decomposition  reaction 

to  carbon  and  sulfur  was  taking  place  and  combined  the  carbon  and  sulfur 

Hugonlots  using  the  stoichiometric  ratio  to  calculate  a  Hugonlot  for  the 
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Temperature, 


Pressure,  kbar 

Fig.  3.2  —Temperature-Pressure  Diagram  Constructed 
by  Butcher  et  al.  in  Ref.  7.  The  cusp  observed  by  R.  Dick 
occurs  well  up  into  the  region  marked  "Decomposition." 


products.  This  curve  went  right  through  Dick's  data  above  the  anomalous 
behavior.  They  also  measured  the  temperature  by  a  brightness  technique 
and  found  It  to  be  3000°K,  about  500°K  higher  than  their  calculations 
Indicate  It  should  be.  Using  both  the  reacted  Hugonlot  data  and  the  temperature 
measurement  they  argued  that  a  decomposition  reaction  was  taking  place. 

A  rather  lengthy  study  on  acrylonitrile  was  done  by  Yakushev,  Nabatov, 
g 

and  Yakusheva  In  which  they  observed  a  cusp  In  the  Hugonlot  at  43  kbar, 
based  on  shock  velocity  measurements.  They  also  used  an  electromagnetic 
particle  velocity  gage  and  measured  what  appears  to  be  a  two-wave  structure 
but  they  said  the  scatter  In  the  data  was  so  great  that  It  was  not  possible 
to  "establish  unambiguously  the  nature  of  the  second  discontinuity  In  the  mass 
velocity."  The  measurement  reproduced  was  12  mm  Into  the  liquid.  Among 
several  other  tests,  recovery  experiments  were  done  where  the  acrylonitrile 
was  shocked  to  240  kbar  and  then  recovered.  Upon  opening  they  discovered 
the  evolution  of  large  amounts  of  ammonia  with  a  black  velvet  powder  remaining 
which  was  Identified  by  x-ray  diffraction  as  hexagonal  graphite.  Using  this 
data  they  speculated  the  decomposition  reaction  Is 

CH2-CH-CN - W  3C  +  NHj. 


This  work  on  acrylonitrile  Is  the  first  good  work  done  on  non-explosives  for 
the  purpose  of  Identifying  a  decomposition  reaction  taking  place  In  a  shock. 

Acrylonitrile  Is  such  a  complicated  molecule  that  the  reaction  mechanism 
could  Involve  several  tens  of  intermediate  reactions  before  the  C  and  NH3  are 
formed,  but  this  work  shows  rather  conclusively  that  a  decomposition  reaction 
Is  taking  place  at  some  point  during  the  shocking  process. 
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Carbon  disulphide  decomposition  has  been  extensively  studied  because 
of  the  simplicity  of  its  linear  triatomic  molecule.1®"1^  The  kinetics  of 
CS2  decomposition  In  the  vapor  phase  has  been  studied  using  shock  tube  and 
optical  spectroscopy  techniques.  The  two  principal  works  are  by  Gaydon 
and  colleagues,  who  measured  absorption  spectra,11  and  Arnold  and  colleagues, 
who  measured  emitted  light.  In  the  former  work,  absorption  was  measured 
for  CS^,  CS  and  S2.  The  latter  two  species  appeared  as  the  former  disappeared, 
and  the  results  were  compatible  with  the  assumption  that  reaction  was  occurring 
according  to 


CS  ♦  s2 


(1) 


1 3 

Arnold  et  al .  measured  light  emission  throughout  the  visible  spectrum. 
They  felt  their  results  to  be  incompatible  with  Eq.  (1)  and  proposed 


C$2  +  CS2 


CS3  +  M 


CS3  +  CS 


CS2  +  S  +  M 


(2) 

(3) 


They  also  found  an  Induction  time  preceding  the  reaction  for  temperatures  below 
2600°K,  whereas  Gaydon  et  al.11  found  no  evidence  of  Induction.  Both  groups  of 
workers  agreed  that  the  reaction  was  second-order  In  CS2  concentration,  with 
reaction  constants 


k 


2 


,0M.f  #-43.7kcal/RT 


Gaydon  et  al. 


k  •  1011*53  e*45.2kcal/RT^ 


(4) 


Arnold  et  al. 


(5) 


In  experiments  where  CS2  vapor  has  been  excited  by  electron  Impact  or 
UV  light.  It  appears  to  be  common  to  observe  spectra  of  CS  and  S  but  not  S2> 

This  Is  peculiar  because  the  heat  of  formation  of  S2  relative  to  S  under 
standard  conditions  Is  about  -34  kcal/mol.  Extrapolation  of  the  decomposition 
rate  for  the  reaction  of  Eq.  (2)  to  the  temperature  and  concentration  produced 
by  a  shock  wave  In  liquid  CSg  leads  to  a  reaction  rate  of  approximately  10'Vusec 
whereas  the  shock  experiments  of  R.  01 ck  (see  Table  I,  Section  1.3)  suggest 
that  It  may  be  about  three  orders  of  magnitude  faster  than  this.  With  this 
much  guidance  from  past  work.  It  Is  evident  that  a  search  for  CS  In  the 
reaction  products  and  a  measure  of  Its  rate  of  appearance  could  be  useful 
In  unravelling  the  basic  kinetics. 

With  these  considerations,  it  was  decided  to  work  on  liquid  CS2»  The 
reasons  can  be  summarized  as  follows.* 

1.  It  has  a  linear  tri-atomic  molecule  and  an  extensive  literature 
on  Its  chemistry  exists. 

2.  Spectroscopic  studies  of  the  kinetics  of  gas  phase  decomposition 
have  yielded  specific  but  controversial  reaction  models. 

3.  Single-shock  Hugonlot  measurements  have  been  made  which  show  a 
cusp  In  the  neighborhood  of  50  kbars,  which  Indicates  chemical 
decomposition  or  a  phase  change. 

4.  Existing  Information  on  the  static  phase  diagram  suggests  that  the 
observed  reaction  Is  a  decomposition. 

5.  The  density  of  liquid  CS2  is  great  enough  that  the  pressures 
and  temperatures  required  for  reaction  can  be  produced  In  a 
double-shock  experiment  with  the  WSU  gas  gun. 
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6.  Absorption  lines  In  CS,  which  appear  to  be  a  most  likely  Interme¬ 
diate  species,  lie  within  the  spectral  range  of  the  time-resolving 
spectrograph  available  to  us. 

7.  CS^  can  be  handled  In  our  shock  facility  without  serious  hazard. 

The  experimental  program  was  planned  in  two  phases.  The  first  was  a  set 
of  mechanical  measurements  Intended  to  map  out  the  regions  of  shock  pressure 
and  temperature  where  reaction  occurred  and  to  determine  reactions  rates  under 
the  various  conditions  achieved.  The  second  was  to  be  a  spectroscopic  study 
which  would  enable  us  to  Identify  reaction  products  and  follow  their  growth 
and  decay  as  reaction  proceeded.  The  first  phase  was  successfully  completed 
and  Is  described  in  Section  IV.  The  second  phase  has  turned  up  some  unusual 
and  unexpected  phenomena  and  Is  still  In  progress.  It  Is  described  In  Section  V. 
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IV.  REACTION  RATES  IN  CS2* 


The  goals  of  this  study  were 

1.  To  determine  as  much  as  possible  about  reaction  mechanisms 
and  kinetics  from  conventional  shock  experiments. 

2.  To  observe  the  evolution  of  a  two  wave  shock  structure  in 
a  chemically  reacting  system. 

3.  To  create  an  equation  of  state  suitable  for  pressure  and 
temperature  calculations  in  strong  shocks  for  off- 
Hugoniot  states . 

4.  To  determine  how  reaction  rate  is  affected  by  changes  in 
the  pressure-temperature  state. 

4.1  Plan  of  Experiments 

Although  it  is  believed  that  the  anomaly  in  the  Hugoniot  recorded  by 
R.  Dick  and  the  associate  behavior  reported  here  result  from  a  decomposition 
reaction,  the  terms  "transition"  and  "reaction"  will  be  used  interchangeably. 

Shock  parameters  were  measured  by  imbedding  several  electromagnetic 
particle  velocity  gages  in  the  liquid  sample  so  that  both  particle  velocity 
and  shock  arrival  times  could  be  continuously  recorded  during  the  times  of 
Interest.  From  these  records  shock  velocity  can  also  be  determined. 

The  particle  velocities  at  which  the  CS^  undergoes  transition  are 
typically  1.35  nsn/ysec  for  a  single  shock  and  greater  for  a  double  shock. 

The  light  gas  gun  available  at  WSU  cannot  produce  1.35  mm/usec  in  the  CS^  in 
a  single  shock,  so  a  double  shock  experiment  was  designed.  The  experimental 
design  is  illustrated  in  Fig.  1.4.  The  target  was  placed  in  a  magnetic  field 

*  Note  that  in  all  of  Section  IV,  Tables,  Figures,  and  Equations  are  labelled 
in  a  consistent  fashion,  but  the  system  does  not  conform  to  that  used  in 
the  other  sections. 


during  the  shot  so  it  was  necessary  to  use  sapphire,  an  insulator,  when  a 
high  impedance  material  was  required. 

As  can  be  seen  from  Figs.  1.4a  and  1.4c,  the  sapphire  faced  projectile 
drives  a  shock  into  the  PMMA  which  then  traverses  the  CS2  as  a  relatively 
low  pressure  wave.  When  the  shock  reaches  the  sapphire  back  plate,  a  high 
pressure  wave  (state  3')  is  driven  back  into  the  already  shocked  CS2-  The 
experiment  can  be  designed  so  the  high  pressure  wave  is  above  the  transition 
pressure  and  a  reaction  is  produced. 

States  1,  2,  3,  and  3',  which  are  reached  at  various  times  during  the 

experiment,  are  labeled  in  Figs.  1.4b  and  1.4c.  These  designations  for  the 
particular  state  they  represent  will  be  used  in  the  remainder  of  this  thesis 
to  maintain  consistency. 

Gages  were  placed  at  the  PMMA-CSj  interface,  in  the  middle  of  the 
CS2>  and  at  the  CS^sapphire  interface.  The  character  of  outputs  from  these 
gages  can  be  inferred  from  Fig.  1.4c.  Figure  1.5a  illustrates  what  the 
particle  velocity- time  profiles  might  be  expected  to  look  like.  Corresponding 
pressure-time  profiles  are  shown  in  Fig.  1.5b,  although  these  profiles  cannot 
easily  be  constructed  from  the  particle  velocity-time  data,  and  no  attempt 
has  been  made  to  do  so.  Pressures  attained  after  each  jump  were  calculated 
from  the  jump  conditions,  however. 

Notice  the  first  reflected  wave  and  the  second  reflected  wave  have 
been  labeled  Che  wave  and  R2  wave,  respectively,  in  Fig.  1.5.  These 
designations  will  be  used  throughout  the  remainder  of  this  section  for  the 
sake  of  brevity. 

Analysis  of  the  gage  records  yields  unrsacted  Hugoniot  data  (state  1) 
as  well  as  information  which  can  be  used  to  calculate  states  2,  3,  and  3'. 


(a)  Cross-section  of  experiment 
Projectile  Target 


Front  Sample 


Distance 


Fig.  1.4.— Experimental  Sketch  and  Shock  Diagrams  for  tha  CS_ 
Experiments.  (a)  shows  tha  axparimaneal  layout,  i.a.  a  sapphire  impactor 
hitting  a  CS2  call  with  a  PMMA  front  and  a  sapphire  back,  (b)  shows  the 
direction  of  tha  shocks  produced  and  tha  various  states  achieved. 


Pressure 


P-Up  Diagram 


Sapphire 


Fig.  1.4  (continued),  (c)  shows  the  states  achieved 
in  pressure-particle  velocity  space.  Notice  the  relaxation 
that  takes  place  from  state  3'  to  state  3  as  the  reaction 
proceeds . 


(a)  Particle  Velocity- time  profiles 


Fig.  1.5.~ Wave  Profiles  Expected  from  a  Three  Gage  Experi¬ 
ment.  (a)  shows  particle  velocity  waveforms  to  be  measured  in  the 
experiment,  (b)  shows  corresponding  pressure  profiles  although 
these  are  not  measured  or  calculated  in  this  study.  The  first  and 
second  reflected  waves  have  labeled  the  R^  wave  and  R2  wave, 
respectively. 
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An  unexpected  noise  associated  with  the  reaction  was  measured  as  part  of  the 


CS^-sapphire  gage  data. 

Reaction  rates  can  be  estimated  from  five  separate  pieces  of  data: 
particle  velocity  decay,  noise  decay,  Induction  time  to  the  start  of  reaction, 
initial  reflected  shock  decay  (R^  wave  decay)  and  wave  structure  of  the 
evolving  Rj  wave  (R^  wave  structure) .  ■  All  the  data  yield  only  estimates 
because  the  exact  nature  of  their  relation  to  the  reaction  rate  is  not  yet 
established  or  because  estimates  had  to  be  made  for  certain  pieces  of  needed 
information. 

Equations  of  state  for  the  unreacted  CS^  and  also  the  reaction  pro¬ 
ducts  were  developed  so  temperature  calculations  and  one-dimensional  wave 
propagation  calculations  could  be  made  to  simulate  the  experiments.  These 
added  some  insight  as  to  what  was  happening  during  the  reaction. 
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E  -  E  -  kP  +  P  )(V  -  V)  (energy)  (2.3) 

o  2  oo 

where  u  is  particle  velocity,  U  is  shock  velocity,  Q  is  density,  V  is 
P  * 

specific  volume  (1/p) ,  E  is  specific  internal  energy  and  P  is  pressure. 

The  subscript  "o"  refers  to  the  condition  ahead  of  the  shock,  and  unsubscripted 
variables  denote  conditions  behind  the  shock.  The  material  is  assumed  to  be 
in  an  equilibrium  condition  at  each  of  these  states. 

The  experiments  outlined  in  Section  1.4  involve  several  waves  (see 
Fig.  1.5).  The  first  shock  propagates  into  a  material  at  rest,  bringing  the 
material  to  state  1  (see  Fig.  1.4).  The  first  reshock  wave  (^  wave)  propa¬ 
gates  in  the  opposite  direction  and  takes  the  material  from  state  1  to  state 
2.  The  second  reshock  wave  (R2  wave)  propagates  behind  the  wav®  and  takes 
Che  material  from  state  2  to  state  3.  Appropriate  jump  conditions  for  each  of 
these  waves  are  tabulated  in  Table  2.1.  It  should  be  pointed  out  that  the 
algebraic  signs  of  the  shock  velocities  and  particle  velocities  must  be  cor¬ 
rect.  For  example,  from  Fig.  1.4b,  shock  velocity  Ufl  is  positive  since  the 
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where  up  is  particle  velocity,  Us  is  shock  velocity,  o  is  density,  V  is 
specific  volume  (1/p),  E  is  specific  internal  energy  and  P  is  pressure. 

The  subscript  "a"  refers  to  the  condition  ahead  of  the  shock,  and  unsubscripted 
variables  denote  conditions  behind  the  shock.  The  material  is  assumed  to  be 
in  an  equilibrium  condition  at  each  of  these  states. 

The  experiments  outlined  in  Section  1.4  involve  several  waves  (see 
Fig.  1.5).  The  first  shock  propagates  into  a  material  at  rest,  bringing  the 
material  to  state  1  (see  Fig.  1.4).  The  first  reshock  wave  (R^  wave)  propa¬ 
gates  in  the  opposite  direction  and  takes  the  material  from  state  1  to  state 
2.  The  second  reshock  wave  (Rj  wave)  propagates  behind  the  R^  wave  and  takes 

the  material  from  state  2  to  state  3.  Appropriate  jump  conditions  for  each  of 

» 

these  waves  are  tabulated  in  Table  2.1.  It  should  be  pointed  out  that  the 
algebraic  signs  of  the  shock  velocities  and  particle  velocities  must  be  cor¬ 
rect.  For  example,  from  Fig.  1.4b,  shock  velocity  is  positive  since  the 
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TABLE  2.1. — Jump  Conditions  Applied  to  Experiments  in  This  Study 


initial  shock  moves  from  left  to  right  but  the  R^  wave  moves  from  right  to 
left  and  so  would  be  negative.  The  particle  velocities  are  all  positive 

as  shown  in  Fig.  1.4c. 

Difficulties  can  arise  in  applying  the  jump  conditions  to  the  R^  and 

R2  waves  since  they  have  their  origin  in  a  single  unsteady  wave  that  evolves 

into  a  two  wave  structure.  If  the  structure  is  not  fully  developed  and  steady, 

the  jump  conditions  in  Table  2.1  do  not  apply  exactly.  In  that  case  the 

relations  among  physical  parameters  can  be  approximated  in  a  progressing  wave 

by  a  multiple  gage  analysis  developed  by  Fowles  and  Williams  1  and  refined  by 
2  3 

others.  *  Such  an  analysis  has  not  been  used  in  the  present  work  because 
it  is  not  evidently  warranted  and  because  Interest  lies  primarily  in  the 
reflected  wave  region  where  techniques  for  multiple  gage  analysis  have  not 
been  developed.  Instead,  the  jump  conditions  have  been  applied  in  every  case 
to  propagating  regions  of  rapid  change  in  particle  velocity  and  it  is  believed 
the  errors  are  not  large  and  do  not  significantly  influence  the  interpretation 
of  records. 

4 

The  time  dependent  conservation  relations  can  be  written  as 

3u  3p 

+  3^  *  0  (momentum)  (2.13) 

|f  +  up(^r>  +  *  0  (mass)  (2.14) 

dt  "  '^dt*  +  dt  (energy )  (2.15) 

where  dQ  is  the  heat  added  per  unit  mass  and  d/dt  ■  3/3t  +  u  (3/3x')  is 

P 

the  convective  derivative.  These  are  most  tractable  when  solved  numerically 
in  the  form  of  difference  equations  in  a  one-dimensional  wave  propagation  code. 
This  is  the  approach  which  has  been  taken  to  simulate  the  experiments. 


Equations  of  State  for  Unreacted  CS„  and  Reacted  Products 

The  jump  conditions  provide  three  equations  among  the  five  unknown 

variables,  P,  u  ,  U  ,  V,  and  E  behind  the  shock  front.  If  two  variables  are 
P  * 

measured  (as  is  the  case  in  this  study)  the  mechanical  state  and  the  P-V-E 

thermodynamic  state  are  determined. 

A  relation  between  P  and  V  for  all  possible  shocks  running  into  a 

state  P  ,V  is  called  a  Rankine-Hugoniot  or  Hugoniot  P-V  curve.  In  what  fol- 
o  o 

lows  it,  along  with  other  relationships  among  shock  parameters  (P-u^,  Ug-Up, 
T-P,  etc.),  will  be  called  simply  the  Hugoniot.  If  it  is  known,  it  provides 
a  fourth  equation,  so  that  the  other  four  parameters  can  be  determined  if  one 
of  the  five  is  measured. 

The  Hugoniot  curve  is  a  single  locus  of  points  on  the  three  dimen¬ 
sional  equation  of  state  surface.  In  order  to  produce  a  numerical  simulation 

based  on  Eqs.  (2.13)  -  (2.15),  a  complete  equation  of  state  is  required. 
Moreover,  since  time-dependent  reactions  will  be  investigated,  a  model  for 
transitions  from  initial  to  final  state  is  required.  Since  temperature  is 
expected  to  play  an  important  role  in  the  reaction  process,  it  must  be  calcu¬ 
lated.  There  is  not  enough  information  about  the  thermodynamic  behavior  of 

\ 

CS  and  its  decomposition  products  to  determine  exact  equations  of  state,  so 
thermodynamic  calculations  will  be  based  on  an  equation  of  state  synthesized 
from  the  limited  data  which  are  available,  following  examples  set  by  Andrews ,5 
Hayes , ^  and  others.?*?  The  remainder  of  this  section  is  devoted  to  a 
description  of  the  equations  of  state  used  here.  Reference  in  is  a  useful 
compendium  of  CS2  properties  and  references. 
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Unreacted  CS„  Equation  of  State 


The  unreacted  CS^  equation  of  state  for  this  study  was  based  on 
Bridgman's^  isotherm  to  12  kbar  and  shock  Hugoniot  data.  Specific  heat  at 
constant  volume  (Cy)  and  the  Gruneisen  parameter  (F)  were  taken  from  the 
thermodynamic  data  at  P  ■  P  ■  1  atmosphere  and  T  ■  TQ  *  293°K.  Functional 
relations  for  Cy  and  r  were  developed  so  that  unreacted  Hugoniot  data  could 
be  matched. 

The  particular  form  chosen  for  the  unreacted  CS^  equation  of  state 

12 

includes  a  Murnaghan  equation  to  describe  the  isothermal  P(V)  relation. 


((%»  .  l} 


(2.16) 


where  K_  is  the  isothermal  bulk  modulus  at  ?  and  T  ,  V  is  the  initial 
To  o  o  o 

specific  volume,  and  N  is  an  adjustable  parameter.  The  isothermal  bulk 
modulus  is  obtained  by  differentiating  Eq.  (2.16): 


*T  * 


(2.17) 


was  assumed  to  have  the  form 


cv  -  cvo(a  +  b(r))  • 

o 


(2.18) 


Since  Cy  ■  CyQ  when  V  ■  V  ,  b  •  1-a  and  Eq.  (2.18)  can  be  written  as 


S  ’  So(m  +  (1"a)V  } 

o 


(2.19) 


There  are  certain  compatibility  relations  between  Cy  end  F  that 

13 

must  be  satisfied  to  maintain  thermodynamic  consistency.  Fowles  developed 


the  relation. 


(2.20) 


When  Eq.  (2.19)  is  substituted  in  Eq.  (2.20)  and  the  result  integrated,  the 
resultant  relation  for  f  is, 


(1-a)  ln(T/T  ) 

r  -  ro  +  - - - —  +f(V)  (2.21) 

1  +  a(y-  -  1) 

where  is  zero  pressure  I\  Tq  is  the  starting  temperature,  and  f(V) 

is  an  arbitrary  function  of  V.  The  final  form  for  r/V  which  was  used  is, 

(1-a)  ln(T/T  ) 
o 

r  v 

(r)v0  (1  +  a(f  - 1)} 

o 

where  a,  c,  and  d  are  adjustable  constants. 

Equations  (2.16),  (2.19),  and  (2.22),  with  initial  values  and  a  sec 
of  constants  a,  c,  and  d,  are  enough  to  calculate  pressure,  volume  relations 
on  the  Hugoniot  curve.  This  was  done  for  various  values  of  a,  c,  and  d,  and 
the  set  which  gave  the  best  fit  was  chosen.  These,  with  initial  values  of 
thermodynamic  variables,  are  given  in  Table  2.2. 


c(l  -  v  ) 


(ia)v  | 

'V  ;  o 
o 


(2.22) 


These  relationships  provide  a  description  of  a  complete  equation  of 
state  surface  for  unreacted  CSg'  Values  for  each  of  the  thermodynamic  vari¬ 
ables  may  be  calculated  by  integrating  on  the  surface  to  the  desired  point. 
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TABLE  2. 2.~ Unreacted  CS^  Equation  of  State  Parameters 
(Reference  State,  P0  *  lbar,  T0  *  293°K) 


Parameter 


Value 

0.624  x  10"5  Mbar  cn3/g9K 

0.01176  Hbar  cm3/g 

1.98  x  10'5  Mbar  cm3/g°K 

0.7918  cm3/ g 

2.108  g/cm3 

0.012  Mbar 

293°K 

7.1 


a  2.4 

c  1.5 


d 
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Reacted  Produces  Equation  of  State 


Obtaining  an  equation  of  state  for  the  reacted  products  was  more  dif¬ 
ficult  than  for  the  unreacted  CS^,  primarily  because  the  exact  nature  of  the 

products  is  not  known.  It  has  been  assumed  that  the  products  are  a  stoichio¬ 
metric  mixture  of  carbon  and  sulfur,  0.1578  gram  of  carbon  and  0.8422  gram  of 
sulfur  per  gram  of  products.  This  assumption  allows  one  to  estimate  what 
some  of  the  equation  of  state  parameters  might  be  in  order  to  have,  a  starting 
point. 

The  form  of  the  equation  of  state  is  much  the  same  as  that  for  the 
unreacted  CS^,  namely,  the  Murnaghan  form  for  the  isotherm  (Eq.  2.16),  a 
constant  C^,  and  a  constant  T/V. 

Dick's  data 14  above  the  transition  region  provided  the  information  to 
which  the  reacted  product  equation  of  state  was  matched.  The  basic  intent  of 
the  matching  procedure  was  to  match  P-V-E  states  obtained  from  Dick's  data  and 
the  energy  jump  condition  with  P-V-E  states  calculated  by  the  reacted  product 
equation  of  state. 

The  final  parameters  determined  for  use  in  the  reacted  product  equa¬ 
tion  of  state  are  listed  in  Table  2.3. 


TABLE  2.3. — Reacted  Product  Equation  of  State  Parameters 
(Reference  State,  PQ  *  lbar,  Tq  ■  293°K) 


Parameter 


Value 

1.0  x  10"5  Mbar  cm3/g°K 

0.0 

0.911  x  lO-5  Mbar  cm3/g°K 

3 

0.480  cm  /g 
1.87  g/cm3 
0.100  Mbar 
293°K 


N 
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Chemical  Thermodynamics 

According  Co  Callen,  ^one  can  postulate  a  fundamental  equilibrium 
relation  in  which  internal  energy  is  a  function  of  entropy,  volume  and  mole 
number  of  the  constituents.  For  a  system  of  CS^  reacting  to  form  sulfur  and 
carbon,  the  fundamental  relation  is 


E  -  E(S,V,NrN2,N3) 


(2.26) 


where  -  Ngg  ,  N2  »  Ng,  and  Ng  -  N 
the  usual  definitions  of  temperature, 
obtain 


g.  One  can  differentiate  this  and  use 
pressure  and  chemical  potential  to 


3 

d!  -  TdS  -  PdV  +  ^  U.dNi  (2.27) 

i-1 

where  is  the  chemical  potential  of  the  iC^  species.  The  mole  numbers 

are  related  through  the  stoichiometry  of  the  reaction  so  that  a  reaction 
coordinate  can  be  defined  as 


d3. 

dN„ 

dN, 

— i  - 

_ 2  , 

3 

vl 

v2 

v3 

(2.28) 


where  x  is  the  reaction  coordinate  and  vi  is  the  stoichiometric  coef¬ 
ficient  (negative  cor  reactants  and  positive  for  products).  This  can  be  used 
in  Eq.  (2.27)  to  obtain 


dE 


TdS  -  PdV  +  dx 


I 

i-1 


U  V 

i  i 


(2.29) 
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where  Che  encire  relation  has  been  divided  by  a  factor  of  proportionality 
which  is  a  constant  so  E,  S,  and  V  are  now  specific  quantities  and  x  is 
dimensionless. 


.15 


The  same  kind  of  development  can  be  carried  out  for  the  Gibbs  free 
energy  representation  to  yield 


dG 


-SdT  +  Vd?  +  dx 


i  v. 


i-l 


(2.30) 


where  G  is  specific  Gibbs  free  energy.  Sometimes  the  summation  term  is 
called  by  another  name 

3 

“t,p  ■  Z  Vi  • 

i-l 

so  Eq.  (2.30)  becomes 


dG  -  -SdT  +  VdP  +  £G_  dx  .  (2.32) 

T,  Jr 

Gibbs  free  energy  is  often  used  as  a  criterion  for  determining  whether  or  not 
a  reaction  at  a  constant  temperature  and  pressure  proceeds;  since,  as  one  can 
see  from  Eq.  (2.30),  dG  is  directly  related  to  dx.  Callen  has  shown  that 
for  a  material  at  constant  temperature  and  pressure  the  Gibbs  free  energy  is 
minimized  at  equilibrium.  This  means  that  if,  at  constant  temperature  and 
pressure,  the  Gibbs  free  energy  of  the  products  is  less  than  that  of  the 
reactants,  the  reaction  tends  to  proceed. 

The  definition  of  Gibbs  free  energy  is 


G  -  E  +  PV  -  TS  .  (2.33) 

In  normal  chemical  reactions  the  PV  part  of  G  is  very  small  (E  is  on  the 

3  3 

order  of  0.02  Mbar  cm  /g  and  PV  is  on  the  order  of  0.000001  Mbar  cm  /g).  As 


pressure  increases,  say  to  10  kilobars,  the  PV  term  becomes  on  the  order  of 
0.01  Mbar  cm3/g.  It  is  easy  to  see  that  in  shock  experiments  very  different 
things  can  happen  to  G  since  both  temperature  and  pressure  increase, 
whereas  in  normal  chemistry  only  temperature  changes. 

The  Gibbs  free  energy  of  formation  can  be  calculated  from  the 
enthalpy  of  formation  and  the  entropy  at  standard  conditions.  The  values 
for  CS^  and  the  products  (carbon  and  sulfur)  are 

AG  (CS,)  -  0.00857  Mbar  cm3/g 

298  L 

and 


3 

AG,  (Carbon  +  Sulfur)  *  0.0  Mbar  cm  /g 

E298 

which  indicates  that  CS^  is  in  metastable  equilibrium  with  respect  to  its 
decomposition  products,  i.e.,  it  tends  to  decompose  at  room  temperature.  The 
rate  is  so  slow  that  it  would  take  years  to  produce  measurable  amounts, 
however.  It  is  interesting  to  note,  that  many  of  the  liquids  chat  react  in 
a  shock  are  metastable  at  standard  conditions. 


Since  the  CS^  is  being  shocked  to  high  pressure  and  temperature 
states  in  the  experiments,  it  is  of  interest  to  calculate  the  Gibbs  free 
energy  for  both  CS^  and  the  reacted  products  to  see  if  the  products  have  a 
lower  Gibbs  free  energy.  This  has  been  done  using  the  equations  of  state 
previously  discussed.  In  both  the  principal  Hugoniot  and  the  reshock 
Hugoniot  regimes  of  Interest,  the  Gibbs  free  energy  of  the  products  is 
always  less  than  that  of  CS2.  Although  this  analysis  was  made  at 
temperatures  and  pressures  along  the  Hugoniot  (which  change  as  reaction 
proceeds) ,  it  does  demonstrate  that  from  an  equilibrium  thermodynamic 
standpoint,  CS2  tends  to  decompose  into  its  products  in  the  state  regions 


of  interest. 


The  only  reaction  rate  law  which  will  be  considered  here  is  that 
for  a  simple  decomposition: 

dx/dt  *  A  exp  (-E/RT) 

where  x  is  the  function  reacted,  E  is  activation  energy,  and  A  is  frequency 
constant.  In  measurements  to  be  reported  here,  the  simplest  rate-related 
variable  is  a  relaxation  time, 

T  •  1/x  *  exp(E/RT)/A 

Graphs  of  In  i  verses  1/Tca^c  will  be  used  to  estimate  E  and  A. 

Reaction  Rate  Data  from  Evolution  of  Two  Wave  Structure 

The  general  theory  of  kinetics  of  shock  induced  phase  transitions  has 

4,5,6,16 

been  addressed  in  detail  by  several  researchers.  Much  of  the 

formalism  applies  directly  to  the  kinetics  of  chemical  reactions.  Evolution 
of  the  two  wave  structure  that  results  from  a  reacting  material  is  illus¬ 
trated  in  Fig.  2.6.  This  figure  represents  "snapshots"  of  the  particle 
velocity-distance  profile  at  several  times  during  the  evolutionary  period. 
According  to  Hayes, wave  profiles  measured  during  the  evolutionary  period 
provide  at  least  four  different  pieces  of  Information,  each  related  to  the 
reaction  rate.  The  four  are:  decay  of  the  first  wave  (called  wave  decay), 
delay  in  forming  the  second  wave  (called  R^  wave  delay) ,  particle  velocity 
transients  at  the  loading  Interface,  and  the  structure  between  the  two  waves 
(called  R^  wave  structure).  It  is  obvious  chat  the  reaction  rate  in  a 
chemically  reacting  wave  is  changing  from  a  relatively  slow  rate  at  the 
beginning  to  a  very  high  rate  toward  the  end  so  that  these  different 
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observables  would  be  expected  to  yield  different  rates  associated  with 
different  times  during  the  reaction.  The  wave  delay  is  difficult  to 
measure  and  will  be  ignored  in  the  following  discussion. 


Wave  Decay 


A  relation  has  been  developed  to  relate  the  particle  velocity  decay 
in  the  wave  to  the  reaction  rate  at  the  Lagrangian  position  where  the 
reactive  wave  starts.  The  relation  is 


(— ) 
'3th 


2 

pa  . 

o  s. ,3u. 

U  J(3t;i 


(2.40) 


-  P(£)xl  -  AE($)X 

x 


w*lere  P0  is  am*>ient  density  of  the  material.  ag  is  frozen  sound  speed. 

Kg  is  frozen  bulk  modulus ,  &V  and  AE  are  changes  in  specific  volume  and  internal 

energy  experienced  by  a  mass  element  in  going  from  the  initial  to  the  reacted 

state.  Subscripts  refer  to  a  position  immediately  behind  the  shock  front  and 

h  to  an  arbitrary  Lagrangian  position.  U  is  shock  velocity  and  T  is  the 

18 

Gruneisen  parameter. 

Suffice  it  to  say  that  this  evaluation  can  be  made  only  at  h  »  0, 

n  a 

where  h— o  1-8  n*8iisibl*  and  is  the  initial  particle  velocity 

decay  in  the  R^  wave.  Since  this  can  be  estimated  from  the  measured  data 

and  all  other  parameters  can  be  calculated  from  the  equation  of  state  of 

3x 

unreacted  CSj  (since  x  •  0  in  the  wave  front) ,  an  estimate  for  (,g^')h^j  can 

be  made. 
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Fig.  2.6.—  Snapshots  of  Che  Particle  Velocity 
Profile  for  an  Evolving  Two  Wave  Structure  in  a  Reac¬ 
ting  Material.  This  illustration  applies  to  the  single 
shock  case.  0  ^  is  the  initial  particle  velocity  as 
the  wave  begins  to  propagate. 


Loading  Interface  Particle  Velocity  Transients 


The  parcicle  velocity  measured  at  the  downstream  CS^-sapphire  inter¬ 
face  (see  Fig.  1.4c)  will  decrease  from  3'  to  3  as  the  reaction  progresses. 

A  relaxation  time  (t)  for  the  reaction  to  be  completed  can  be  obtained  from 
the  particle  velocity-time  history.  The  reciprocal  of  this  time  provides  an 
estimate  of  the  average  value  of  the  reaction  rate  at  this  interface, 

x  =  i  .  (2.41) 


Wave  Structure 


Quite  obviously,  the  reaction  rate  varies  a  great  deal  from  the  top 
of  the  wave  to  the  top  of  the  R2  wave,  starting  out  slowly  and  then 
increasing  rapidly  as  the  temperature  increases.  The  slope  of  the  wave 
front  is  most  likely  where  the  major  part  of  the  reaction  occurs;  it  there¬ 
fore  provides  another  estimate  of  the  reaction  rate.  If  one  multiplies  this 
slope  by  the  reciprocal  of  the  change  in  velocity. 


Au 
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(2.42) 


an  estimate  of  reaction  rate  results.  Obviously  this  assumes  all  the  reac¬ 
tion  occurs  in  the  R^  wave  front  and  it  would  therefore  be  expected  to  yield 
a  low  upper  limit  on  the  reaction  rate. 


Particle  Velocity  Measurements  with  Electromagnetic  Gages 


The  electromagnetic  gage  consists  of  a  conductor  placed  perpendicular 

to  a  magnetic  field.  As  the  conductor  moves,  it  cuts  the  field  lines  and 
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acts  essentially  as  a  voltage  source  that  obeys  the  relation 

£  ■  Biv  x  10  4  (2.43) 

where  £  is  induced  voltage  in  volts,  B  is  the  magnetic  field  strength  in 
gaus,  i  is  the  length  of  the  conductor  in  millimeters,  and  v  is  the  con¬ 
ductor  velocity  in  ram/jjsec. 

There  are  several  things  one  must  keep  in  mind  when  using  these  gages 
in  a  shock  environment: 

1.  Effects  of  gage  material  and  thickness  on  the  material 
being  measured. 

2.  Movement  of  leads  at  the  edges  of  the  experiment. 

3.  Mutual  Inductance  between  gages  in  multiple  gage  experiments. 

Each  of  these  will  be  briefly  discussed  in  the  following  subsections. 

Gage  Material  and  Thickness 

Since  the  electromagnetic  gage  is  an  in-material  gage,  Che  shocking 
up  process  between  the  gage  and  the  material  is  of  interest.  Figure  2.7  is 
a  representation  of  this  shocking  up  process  between  a  high  impedance  gage 
and  a  low  impedance  medium.  Notice  Che  rounding  in  the  front  that  continues 
past  the  gage  and  the  perturbation  traveling  into  the  already  shocked  region 
on  the  left.  If  the  gage  is  thin  enough  the  rounded  wave  will  rapidly  shock 


Time 


a; 

3 

CO 

</) 

03 

v. 

CL 


Fig.  2.7. — Electromagnetic  Gage  Shock-Up  Process. 

(a)  shows  the  waves  generated,  (b)  shows  che  pressure-particle 
velocity  states  attained,  (c)  is  a  snapshot  of  the  pressure  profile 
showing  the  rounded  wavefront  proceeding  past  the  gage  and  the  per¬ 
turbation  moving  to  the  left  in  the  already  shocked  material.  These 
diagrams  apply  to  a  high  impedance  gage  in  a  low  impedance  medium. 
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up  as  the  small  x*aves  overtake  the  front  and  the  perturbation  traveling  to 
the  left  will  be  depleted  by  the  rarefactions  from  behind,  thus  leaving  the 
medium  in  the  desired  state. 


When  dealing  with  temperature  sensitive  reactions,  one  must  be  cog' 
nizant  of  the  possibility  that  reaction  may  occur  in  the  high  pressure  per¬ 
turbation  moving  to  the  left,  thereby  causing  a  much  larger  perturbation 
chat  affects  the  whole  experiment. 

The  ideal  gage  would  be  one  that  matches  the  medium  in  impedance. 

The  alternative  is  to  make  it  very  thin,  i.e. ,  a  few  microns  chick. 

Lead  Movement 

The  leads  connecting  the  conductor  to  the  recording  circuit  must 
come  out  parallel  to  the  magnetic  field  (or  else  in  a  position  so  chat  it 
does  not  cut  the  field  as  it  moves)  to  maintain  a  constant  effective  conduc¬ 
tor  length.  As  the  experiment  progresses,  the  outside  edges  of  the  experi¬ 
ment  lose  one-dimensionality  and  it  is  conceivable  the  leads  begin  to  spread. 
This  spreading  would  manifest  itself  on  the  records  as  a  gradual  increase  in 
the  particle  velocity,  which  should  remain  constant.  These  type  of  results 
have  been  observed  but  as  yet  no  one  has  attempted  to  quantify  the  effect. 

The  records  from  the  experiments  in  this  study  indicate  this  was  not  large 
enough  to  worry  about. 

Multiple  Gage  Mutual  Inductance 

Using  very  simple  relationships  for  the  mutual  inductance  between 
two  gages  in  a  multiple  gage  experiment,  the  worst  case  was  found  to  be  two 
gages,  one  on  top  of  the  other,  separated  by  a  layer  of  the  material  being 
shocked.  For  gages  2  mm  apart,  the  mutual  Inductance  was  calculated  to  be 
0.025  x  10’6  henries.  Voltage  is  generated  between  these  two  circuits  only 


when  Che  current  in  one  or  che  ocher  is  changing  since 


£  ■  M  “  (2.44) 

where  £  is  Che  induced  volcage,  M  is  che  mutual  inductance  and  i  is  Che 
current.  Assuming  a  current  change  of  20  milliamps  in  10  nanoseconds  leads 
to  an  induced  volcage  of  50  millivolts.  This  would  last  only  while  the  cur¬ 
rent  was  changing  and  should  therefore  be  manifest  as  spikes  on  the  records. 
Some  of  these  have  been  observed,  but  they  do  not  alter  che  data  appreciably 


4.3  EXPERIMENTAL  PROCEDURE 


General  Comments  on  Experimental  Facility 

All  experiments  were  performed  on  the  44  foot  long,  4  inch  bore 
light  gas  gun  operated  by  the  Shock  Dynamics  Laboratory  at  '.Washington  State 
University.  The  characteristics  and  performance  of  the  gun  are  documented 
in  Ref.  20.  It  is  capable  of  delivering  a  projectile  at  velocities  ranging 
from  0.2  to  1.35  mm/ysec,  although  the  higher  velocities  are  somewhat  more 
difficult  to  attain  because  the  projectile  must  withstand  the  harsh  environ¬ 
ment  of  high  pressure  (6000  PSI)  while  having  a  mass  of  only  half  a  kilogram. 

-3 

The  target  area  in  the  gun  is  evacuated  to  less  than  10  Torr  prior 
to  firing  to  eliminate  the  air  cushion  which  would  otherwise  form.  This 
area  is  large  enough  to  accommodate  an  electromagnet  for  experiments 
requiring  particle  velocity  measurements.  Adjustments  can  be  made  on  the 
target  assembly  so  the  projectile-target  tilt  is  controlled  to  less  than  a 
milliradian. 

Projectile  velocity  is  measured  as  a  result  of  the  projectile  shor¬ 
ting  4  sets  of  velocity  pins  just  prior  to  impact.  Signals  from  these 
shorting  pin  pairs  are  routed  through  a  pulse  shaping  network  to  both  a  time 
interval  counter  and  an  oscilloscope.  Velocity  measurements  on  the  order 
of  to. 57.  are  possible  with  this  system. 

The  instrumentation  facility  consists  of  nine  Tektronix  oscilloscopes 
(one  7844,  two  7904* s,  two  454’ s,  and  four  585 ’s)  which  can  all  be  used  at 
one  time  if  necessary.  The  7000  series  scopes  have  interchangeable  ampli¬ 
fiers  (7A13  and  7A19)  and  interchangeable  time  bases  (7B70  and  7B71). 

Signals  from  the  experiment  are  carried  from  the  target  chamber  to 


the  instrumentation  rack  via  four  RG213U  coaxial  cables  isolated  from 
ground. 


CS^  Cell  Design  and  Fabrication 

The  design  and  fabrication  of  the  cell  to  contain  the  CS2  evolved 
as  experimentation  progressed,  primarily  due  to  problems  encountered  and  a 
desire  to  obtain  as  much  information  as  possible  from  each  experiment. 

Design,  fabrication,  and  problem  areas  will  be  addressed  in  the  next  two 
subsections . 

Design  and  Construction  of  the  Cell 

During  the  course  of  experimentation,  five  different  cell  designs 

were  used.  The  differences  are  outlined  in  Table  3.1.  All  cells  were 
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constructed  from  layers  of  polymethylmethacrylate  (PMMA)  which  were  lapped, 
polished,  and  epoxied  together.  PMMA  was  chosen  because  it  is  transparent, 
easy  to  work,  with,  and  a  reasonably  good  shock  impedance  match  to  CS^, 
particularly  at  higher  pressures.  It,  like  all  plastics,  is  a  viscoelastic 
material  which  gives  rise  to  steady  rounded  shocks  propagating  through  it. 

A  pictorial  breakdown  of  the  parts  which  make  up  a  cell  of  the  C 
type  is  shown  in  Fig.  3.1.  The  gages  were  mounted  on  kapton  sheets  prior  to 
assembly.  Layers  were  epoxied  together  using  Epon  815  epoxy  allowing  approx¬ 
imately  24  hours  for  each  epoxy  joint  to  cure.  Over  12  separate  epoxy 
operations  were  required  in  the  construction  of  the  more  complicated  cells 
such  as  the  C  type. 

Sapphire  back  plates  were  single  crystal  sapphire,  unoriented,  but 

22 

most  likely  in  the  60°  orientation.  They  were  purchased  this  way  to  keep 
the  costs  down  since  oriented  z-cut  sapphire  costs  approximately  two  times 
more. 
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TABLE  3.1. — CS2  Cell  Configurations 


Front  Face 
Material 

Back  Face 
Material 

dumber 
of  Gages 

Placement 

of 

Gages* 

Designation 

F 

MF  M 

MB 

B  BS 

A 

PMMA 

PMMA 

2 

X 

X 

B 

PMMA 

Sapphire 

2 

X 

X 

C 

PMMA 

Sapphire 

3 

X 

X 

X 

D 

PMMA 

Sapphire 

4 

X 

X 

X  X 

E 

PMMA 

Sapphire 

4 

X 

X 

X 

X 

*  F  *  Front  gage  at  CS2-PMMA  interface. 

MF  *  Gage  in  CS^  approximately  0.7  mm  from  front  face. 

M  *  Gage  in  CS^  approximately  1.0  mm  from  front  face. 

MB  ■  Gage  in  CS^  approximately  1.3  mm  from  front  face. 

B  *  Back  gage  at  CS^-Sapphire  interface  approximately  2mm  from 

front  face. 

BS  *  Gage  between  two  layers  of  sapphire  in  the  back  plate. 


Lead  wires 


Sapphire  back 
plate 


Exploded  view 
of  cell 


'PMMA  plate  for 


PMMA-sapphire 
cell  back 


Extra  pour  of 
epoxy  here  to 
strengthen 
epoxy-sapphire 
bond 


Fig.  3.1. — Exploded  View  of  a  Three  Gage  CS^  Cell.  The  cell  is 
assembled  by  epoxying  che  various  layers  together. 
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Compatibility  of  the  CS^  with  cured  Epon  815  epoxy,  kapton  and  P’CIA 
was  checked  by  placing  these  materials  in  CS^  tor  several  days  and  then 
looking  carefully  for  visual  signs  of  degradation.  They  were  all  found  to 
be  compatible,  although  no  attempt  to  quantify  the  results  was  made.  One 
epoxy,  Dow  DER332  resin  with  Hysol  3561  hardener,  was  found  to  be  incom¬ 
patible  wich  CS^. 

23  24 

The  electromagnetic  gages  were  designed  by  Fowles  and  Roller  *  of 

25 

this  laboratory  and  are  purchased  mounted  on  a  50  micron  thick  kapton 
backing  which  is  later  removed.  The  gages  are  5  microns  thick  copper  with 
the  leads  and  end  0.5  mm  wide.  Effective  length  (approximately  6.5  mm)  is 
determined  by  taking  an  average  of  the  inside  and  outside  length  of  the  end. 
The  leads  are  25.4  mm  long.  There  is  a  center  tap  to  be  used  with  differ¬ 
ential  amplifier  measurements,  if  desired.  Prior  to  assembly  in  the  cell, 
the  gages  were  epoxied  onto  7  micron  thick  kapton  to  make  it  easier  to  handle 
them  and  also  to  provide  support  for  the  gages  suspended  in  the  CS^. 

Access  to  the  cell  for  filling  was  provided  by  two  holes  drilled  in 
the  side  and  into  the  sample  area  in  the  PMMA  back  plate  ring.  Soft  copper 

tubes ^  1.53  mm  outside  diameter  by  0.56  mm  inside  diameter  were  inserted  in 
the  holes  but  not  into  the  sample  cavity.  These  were  epoxied  in  at  the  time 
the  cell  was  potted  into  the  target  plate.  Soft  copper  was  used  because  it 
is  non-magnetic  and  easy  to  seal  by  crimping  after  the  liquid  has  been 
inserted. 

Figure  3.2  is  a  sketch  of  the  finished  ta.-get  assembly  ready  to  be 
filled  with  CS2.  A  number  of  construction  details  h»ve  been  covered  only 
briefly  and  several  others  (such  as  attaching  leads  .  nd  cables  to  the  gages, 
building  up  the  leads  with  solder,  potting  the  cell  in  the  target  ring, 
testing  the  finished  cell,  etc.)  have  been  skipped  entirely. 


There  were  a  number  of  characteristics,  other  than  those  specified 
in  Table  3.1,  which  changed  from  cell  to  cell.  These  differences  are 
detailed  in  Table  3.2  for  each  of  the  14  cells  used  in  the  CS^  shots. 

Difficulties  Encountered 

A  number  of  problems  were  encountered  during  the  course  of 
experimentation  and  changes  were  made  to  correct  the  deficiencies  on  later 
cells.  Table  3.3  is  a  listing  of  the  problems,  reasons  for  the  problems, 
and  the  solution  developed  to  overcome  it. 

Projectile  Design  and  Preparation 

In  these  experiments,  the  primary  function  of  the  projectile  body 
was  to  hold  the  sapphire  impactor  in  such  a  fashion  that  it  would  impact 
the  cell  with  a  minimum  amount  of  tilt  as  illustrated  in  Fig.  3.3.  Because  of 
the  nature  of  the  experiments,  stringent  constraints  were  put  on  the  pro¬ 
jectile  design.  Since  the  experiment  was  carried  out  in  a  magnetic  field, 
the  projectile  must  be  nonconductive  or  at  least  have  a  nonconductive  nose 
on  it.  The  high  velocities  necessary  to  achieve  the  desired  state  in  the 
CS^  required  the  projectiles  be  light  but  at  the  same  time  able  to  withstand 
the  extremely  harsh  environment  of  being  accelerated  by  6000  psi  helium. 

Projectile  Design 

Three  different  projectile  designs  were  tried  which  are  illustrated 
in  cross-sectional  views  in  Fig.  3.4.  The  all  nylon  projectile  (Fig.  3.4a) 
was  abandoned  because  it  did  considerable  damage  to  the  magnet  interior  when 
it  shattered.  The  aluminum-PMMA  projectile  (Fig.  3.4b)  was  only  good  for 
low  velocity  shots  because  it  was  too  heavy.  Nylon-syntatic  foam  projec- 
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Copper  Tubes 


4  Each 
RG58A/U 
A  Cables 


Epoxy 


CS*  Cel 


PMMA^ 
Tubing 


Epoxy 

Filled 


IISMUlig 


Section  A-A 

Fig.  3.2.— Finished  Four  Gage  Target  Assembly.  The  g-'.ges 
are  visible  from  the  back  of  the  target  so  the  gage  leads  can  be 
visually  aligned  .  to  make  them  parallel  to  the  magnetic  field. 
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tiles  (Fig.  3.4c),  patterned  after  an  SRI  design,  were  used  for  nearly 
all  the  shots . 

Each  projectile  had  a  brass  trigger  pin  shorting  strip  (1.5  mm  thick 
by  10  mm  wide,  covering  120°  of  arc)  positioned  on  the  outside  edge  as  shown 
in  Fig.  3.4.  Apparently  the  projectile  did  not  rotate  appreciably  during  its 
travel  down  the  barrel  since  the  strip,  with  only  one  exception,  shorted  the 
velocity  pins. 


Projectile  Preparation 


In  order  to  minimize  tilt  at  impact,  the  projectile  front  was  lapped 
and  checked  for  perpendicularity  with  the  sides.  The  maximum  acceptable 
tolerance  was  0.01  milliradians  out  of  perpendicularity.  This  process  was 
completed  prior  to  the  installation  of  the  sapphire  impactor  because  of  the 
difficulty  in  lapping  sapphire. 

Sapphire  used  for  the  impactor  was  of  the  same  type  used  for  the  cell 

22 

back  place — unoriented,  single  crystal.  Dimensions  for  the  impactors,  along 
with  other  projectile  properties,  for  each  shot  are  listed  in  Table  3.4. 

Figure  3.5  shows  the  impactor  before  installation  (grooves  were  cut  to  give 
the  epoxy  something  to  key  on)  and  after.  It  was  epoxied  into  the  projectile 
with  Epon  815  epoxy.  Notice  the  rubber  0-ring  which  was  placed  in  the  bottom 
of  the  cutout  to  hold  the  sapphire  firmly  against  a  lapped  flat  steel  plate 
which  was  weighted  (about  10  Kg)  to  make  it  sit  on  the  lapped  projectile  face 
while  the  epoxy  cured.  Measurements  indicated  the  impactor  was  parallel  with 
the  projectile  face  to  within  0.06  milliradian. 

Projectile  mass  was  a  big  consideration  and  the  nylon-syntactic  foam 
projectiles  were  found  to  be  the  lightest  ones  yet  designed  for  this  gun. 
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Velocities  of  1.2  mm/ysec  were  easy  to  attain  with  1.30  mm/usec  possible  (and 
achieved  in  another  set  of  experiments).  The  mass  of  each  of  the  projectiles 
used  in  this  study  is  given  in  Table  3.4. 


(a)  Nylon  projectile  configuration 


TABLE  3.4. —  Projectile  Properties 
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(a)  Sapphire  before  installation 


Projectile  Face 

^  ^-Sapphire  Impactor 


(b)  Sapphire  impactor  glued  in  place 


Fig.  3.5. — Sapphire  Impactor  Installation  in 
Projectile.  The  rubber  0-ring  was  used  to  keep  the  impac¬ 
tor  flush  with  the  projectile  face  during  the  epoxy  opera¬ 
tion.  Slots  were  made  in  the  impactor  for  the  epoxy  key  on. 


4.4  RESULTS 


Results  of  the  experiments  on  CS2>  along  with  calculations  based  on 
the  equations  of  state  introduced  in  Section  4.2,  will  be  present  in  this 
chapter.  Also  included  are  the  results  of  calculations  using  a  one-dimensional 
wave  propagation  code  to  simulate  some  of  the  experiments.  A  complete 
discussion  of  both  the  calculated  and  experimental  re^ilts  will  be  presented 
in  the  next  chapter. 

General  Shot  Data 

Sixteen  usable  experiments  were  performed  in  the  course  of  this  study; 
14  with  CS  ,  one  with  H-0  and  one  with  a  cell  made  entirely  of  PMMA.  The 
general  information  for  each  of  these  shots  is  listed  in  Table  4.1.  Particle 
velocity  waveforms  were  time-correlated  and  then  plotted  so  the  various  infor¬ 
mation  could  be  obtained  for  the  initial  shock  and  the  two  reflected  shocks. 
Time-correlated  data  for  shots  77-106,  77-070,  and  77-020  are  shown  in  Figs.  - 
4.2,  and  4.3,  respectively.  Notice  in  Fig.  4.1  the  particle  velocity  decrea'r- 
as  each  of  the  two  reflected  shocks  passes  the  gage  plane.  Also  net:..e 
fiducial  marks  that  were  used  to  time-correlate  the  data. 

Two  data  reduction  problems  will  be  discussed  at  this  -  . 

the  cells  were  filled,  the  intermediate  kapton  gage  planes  wr-.'.#  .  * 
so  final  gage  positions  were  considerably  different  fro*  »♦  •«*  , 
the  cell  was  constructed.  Because  of  this  solidly 
gages  were  used  to  determine  the  initial  cho  n  «•  . 
intermediate  gage  positions  were  back  a.  ..  . - » 
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TABLE  4.1. — General  Shot  Information 
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'Positioned  by  using  shock  velocity  computation. 


Correlated  Particle  Velocity  Waveforms  for  Shot  77-020.  F  and  B  are  the 
id  the  back  (CS  -PUMA)  gage  signals,  respectively. 


The  second  problem  is  readily  observable  from  Fig.  4.1  and  4.2,  ehae 


of  different  scopes  yielding  a  different  particle  velocity  waveform  even 
though  they  were  hooked  to  the  same  gage.  Obviously  a  number  of  things  might 
cause  this  result:  bad  calibrations  from  scope  to  s.:ope,  different  kinds  of 
scopes  recording  the  same  pulse  differently,  D-C  calibration  techniques  being 
used  to  calibrate  even  though  pulses  are  being  recorded,  base  line  shifts,  etc. 
It  is  likely  that  all  of  these  are  involved  to  a  greater  or  lesser  extent  to 
make  the  particle  velocity  amplitudes  vary  by  up  to  10Z.  Two  mock  shots  were 
carried  out  using  a  pulse  generator  to  simulate  the  waveforms  in  regular 

experiments.  The  results  were  essentially  identical  to  the  shot  data, 
i.e.,  voltages  of  the  pulses  varying  by  up  to  10 Z.  The  decision 
was  then  made  to  average  the  particle  velocity  amplitude  data  and  use  the 
mean  as  the  measured  value.  Data  obtained  by  doing  this  on  the  CS^,  E?0, 
and  PMMA  shots  are  consistent  with  results  of  other  experiments,  as  will 
be  discussed  later. 

Unreacted  Hugoniot  Data 

In  each  experiment  the  initial  shock  is  nonreactive  and  therefore 
yields  an  unreacted  Hugoniot  point.  Fourteen  new  Hugoniot  points  between  17 
and  SO  kbar  were  generated  in  this  study  and  are  tabulated  in  Table  4.2. 
Maximum  and  minimum  values  of  particle  velocity  for  state  1  are  tabulated 
along  with  the  mean  value  which  was  used  in  the  calculations  to  determine  P^, 
V1§  and  E^.  One  particle  velocity  for  each  waveform  was  used  in  the  averaging 
procedure.  The  appropriate  value  for  the  intermediate  gage  waveforms  was 
apparent  because  the  waveforms  were  essentially  flat  on  top.  CS^-PMMA  gage 


TABLE  4.2. — CS  Unreacted  llugonlot  Data,  State 


'Hay  have  been  reacting  in  first  shock;  rcsltock  pressure  above  elastic  limit  of  sapphire. 


values  ware  harder  co  pick  because  of  the  viscoelastic  nature  of  the  PMHA.  A 
value  just  above  the  major  part  of  the  wave  front  rounding  was  used.  Figure 
4.4  illustrates  the  waveform  position  of  the  particle  velocities  used,  A 
being  the  value  for  the  PMKA-CS^  gage  and  B  being  the  value  for  the  inter¬ 
mediate  gages. 

Shock  velocity  was  calculated  from  the  front  and  back  gage  data  as 
previously  mentioned.  Since  both  shock  velocity  and  particle  velocity  were 
measured  quantities,  no  other  information  was  required  to  completely  deter¬ 
mine  the  P-V-E  values  of  state  1.  It  was  unnecessary  to  know  the  PMMA 
Hugoniot  exactly  or  to  use  the  projectile  velocity  in  obtaining  this  state. 

Unreacted  CS„  Reshock  States 

1  ■  -c 

Reshock  states  in  the  CS^  were  obtained  from  the  measured  data  in 
two  ways.  For  the  shots  that  did  not  immediately  react,  the  initial  particle 
velocity  measured  by  the  CS2-sapphire  gage,  along  with  the  known  sapphire 
Hugoniot,  yield  the  reshock  state  3'  pressure  and  particle  velocity  (Fig. 
1.4c).  These  data  for  six  of  the  shots  are  tabulated  in  Table  4.3  along  with 
state  1  data.  Note  that  shots  77-003,  77-020,  and  77-044  do  not  appear  in 
reshock  tables  because  they  had  PMMA  back  plates.  Also  included  in  the  table 
are  calculated  reshock  states  based  on  the  unreacted  CS^  EOS.  The  calcula¬ 
tions  were  made  using  the  Initial  state  P^  as  the  reference  point  for  the 
reshock  calculation.  The  final  state  (31)  was  obtained  graphically  from  the 
point  where  the  left  going  reshock  Hugoniot  crossed  the  right  going  sapphire 

Hugoniot  as  shown  in  Fig.  1.4c.  The  sapphire  Hugoniot  used  in  the  calcula- 

28 

tlons  was  that  measured  by  Barker  and  Hollenback, 

U#  -  1.119  +  1.0  u  (cm/usec)  .  (4.1) 


18 


Fig.  4.4. — Positions  of  State  1  Particle  Veloc¬ 
ities  on  Waveforms.  A  was  the  position  chosen  for  the 
front  (PMMA-CS2)  gage  signal  and  B  was  the  position 
chosen  for  the  intermediate  gage  signals.  Note  that 
the  intermediate  gage  signals  are  essentially  flat- 
topped  while  the  front  gage  signals  show  structure  due 
to  the  viscoelastic  nature  of  PUMA. 
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The  second  method  used  to  obtain  reshock  states  was  to  use  the 


measured  state  2  data  for  the  R^  wave.  Shock  velocity  was  determined 
graphically  from  a  t-x  diagram.  Particle  velocity  was  taken  as 
the  mean  of  those  measured  at  the  top  of  the  wave;  one  value  for 
each  scope  record  was  used  in  the  mean  calculation.  Pressure  and  specific 
volume  were  calculated  using  Eqs.  (2.7)  and  (2.3)  with  state  1  that  given  in 
Table  4.3.  Table  4.4  contains  the  resulting  FL  wave  data.  It  is  believed 
that  these  data  are  not  as  reliable  as  those  obtained  by  the  first  method 
because  of  wave  evolution  and  also  because  of  inaccuracies  in  the  particle 
velocity  measurements  which  are  compounded  each  time  the  next  calculation  is 
made.  Also  included  in  Table  4.4  are  calculated  reshock  states  made  by  using 
the  unreacted  CS^  EOS  with  P^  as  the  reference  pressure  and  P^  as  the  final 
pressure.  Temperatures  are  only  approximate,  but  no  errors  can  be  assigned 
because  of  EOS  uncertainties. 

Reacted  Product  Data 

There  are  two  methods  to  obtain  reacted  product  information  from  the 
experimental  measurements.  One  is  to  use  the  relaxed  particle  velocity 
(state  3  in  Fig.  1.4c)  from  the  CS2~sapphire  gage  data  and  the  other  is  to 
calculate,  from  measured  shock  and  particle  velocities,  state  3  which  lies 
behind  the  R2  wave.  As  previously  indicated,  the  first  method  gives  the  most 
reliable  information  because  it  is  a  direct  measurement. 

State  3  particle  velocities  from  the  CS2-3apphire  gage  are  tabulated 
in  Table  4.5  along  with  the  pressure  calculated  using  the  sapphire  data 
(Eq.  4.1)  in  the  momentum  jump  condition.  Also  Included  in  this  table  are 
state  3  values  calculated  from  the  measured  R2  wave  data.  The  shock  velocity 
was  determined  graphically.  Pressure  and  specific 


91 


TABLE  4.4. — Calculated  and  Measured  llnreacted  CS  R  Wave,  State  2,  Reshock  States 


Gage  broken  prematurely ,  no  R1  wave  data  measured. 


TABLE  4.5. — Measured  State  3  for  Reacted  Products 


volume  were  calculated  using  Eqs.  (2.10)  and  (2.11)  along  with  state  2 


values  determined  by  the  R^  wave  analysis. 

Relaxation  Data  from  CS,,-Sapphire  Gage 

Measurements  at  the  CS^-sapphire  gage  yielded  some  interesting  data 
concerning  the  reaction.  Figure  4.5  is  an  idealized  drawing  of  a  gage  record 
showing  the  various  features  and  data  available.  At  the  time  of  reaction 
onset  a  high  frequency  noise  was  measured  that  had  a  decay  time  longer  than 
the  particle  velocity  decay  time.  Because  of  the  noise  it  was  impossible  to 
determine  the  exact  nature  of  the  particle  velocity  decay.  Data  for  this 
decay  were  estimated  by  drawing  a  curve  through  the  noise  at  what  appeared 
to  be  the  mean  position.  Noise  decay  time  was  also  estimated  but  it  was  dif¬ 
ficult  to  tell  the  end  of  the  noise  from  normal  noise  on  the  approximately 
100  millivolt  signals.  For  many  of  the  shots  the  induction  time  was  so  small 
that  it  was  not  resolvable  in  the  10  to  20  nsec  risetime  of  the  gage. 

Waveforms  for  some  of  the  CS^-sapphire  gage  measurements  are  contained 
in  the  data  records  of  Sec.  4.1,  but  the  details  cannot  be  determined  from 
them.  The  measurements  and  estimations  were  made  on  expanded  waveforms. 
Information  from  this  exercise  is  tabulated  in  Table  4.6  along  with  calculated 
temperatures  (T^,)  and  pressures  (P^,)  from  Table  4.3  as  a  reference. 


Particle  Velocity  Relaxation  in  the  R^  Wave  Front 

As  indicated  in  Sec.  4.2,  the  initial  particle  velocity  decay  rate 

in  the  R,  wave,  Qu  /3t)_,  is  needed  to  make  the  calculation  to  estimate  the 
1  P  5 

reaction  rate.  This  can  be  obtained  by  plotting  changes  in  particle  velocity 
as  the  wave  propagates.  It  starts  out  with  a  small  value  (state  3')  which 
increases  as  it  travels  because  the  wave  is  left-traveling  in  a  medium  moving 
to  the  right.  Figure  4.6  is  a  plot  of  R^  wave  front  particle  velocity  versus 
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Time 


Fig.  4.5.— Idealized  CS2~Sapphire  Gage 
Record,  tj  is  Che  induction  time,  tN  is  Che 
noise  decay  time,  tu  is  the  particle  velocity 
decay  time,  and  Aup  Is  the  particle  velocity 
decrease. 


TABLE  4.6. — Reaction  Data  from  CS.-Sapphire  Gage 


Estimate  baaed  on  calculated  state  3'  particle  velocity  (see  Table  4.3)  as  the  initial  state. 
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Pig.  4.6. — Rj  Wave  Decay  for  Shot  77-107.  The  initial  part  of  the  decay 
ve  is  estimated  as  represented  by  the  dashed  line.  The  value  chosen  for  (iu  /3t) 
this  shot  is  2.S  mm/psec^. 
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cine  for  shot  77—107.  The  nethod  used  to  estimate  che  desired  decay  race  is 
illuscrated  in  che  figure.  A  decay  race  obcained  by  using  che  inicial  par- 
Cicle  velocity  and  che  parcicle  velocity  of  the  closest  gage  to  che  back  is 
obviously  too  small.  If  one  estimates  the  shape  of  the  curve,  it  is  apparent 
that  two  to  three  times  this  rate  is  more  reasonable.  Table  4.7  is  a  tabu¬ 
lation  of  the  various  particle  velocities  and  the  associated  times  after 
reflection  from  the  sapphire,  along  with  what  is  considered  a  reasonable 
estimate  of  (ou^/3t)^  at  the  CSj-sapphire  interface.  Only  the  shots  with 
intermediate  gages  were  examined. 

Rise  Time  of  the  Wave 

The  rise  time  of  the  R£  wave  yields  an  estimate  of  reaction 
rate.  A  review  of  the  waveforms  indicates  the  R£  wave  front  is  usually 
rounded  at  both  the  cop  and  bottom  of  the  wave,  so  only  the  middle  slope, 
near  the  inflection  point,  was  considered  in  the  analysis.  It  is  also 
apparent  the  wave  fronts  were  not  always  smooth  so  some  judgment  was 
required  in  deciding  what  to  use. 

Results  of  this  analysis  are  listed  in  Table  4.8  along  with  the  reac¬ 
tion  races  calculated  from  the  average  value  of  Au  /At  divided  by  the 

P 

particle  velocity  difference  from  the  top  to  che  bottom  of  the  wave.  Note 
that  some  of  the  waves  steepen  as  they  travel  while  others  spread  out.  Any 
change  implies  that  the  wave  is  not  steady  and  that  the  jump  conditions  used 
in  calculating  states  2  and  3  do  not  strictly  apply.  There  are  not  enough 
data  to  say  anything  about  a  trend  in  this  area.  Some  waveforms  exhibited 
more  than  one  slope,  which  may  indicate  more  than  one  reaction  in  the  wave. 

Miscellaneous  Experiments 

Several  experiments  were  run  to  determine  the  nature  of  certain 
problems  and  anomalous  effects.  One  shot  (77-104)  consisted  of  a  PMMA  cell 
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TABLE  4.8.— R  Wave  Front  Reaction  Rate  Data 
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with  PMMA  between  the  gages;  another  shot  (77-105)  consisted  of  a  regular 
cell  filled  with  distilled  water;  and  finally  two  shots  (77-094  and  77-095) 
had  gages  between  two  pieces  of  sapphire  to  look  at  the  noise  generated  by 
the  reaction.  These  shots  are  discussed  in  the  following  subsections. 

PMMA  Shot 

Shot  77-104  consisted  of  a  regular  cell  with  a  sapphire  backplate 
but  without  the  cutouts  for  the  liquid,  i.e.,  PMMA  layers  were  between  the 
gages.  The  objectives  of  this  shot  were,  first,  to  determine  if  any  extra¬ 
neous  electrical  signals  were  generated  in  the  PMMA  that  might  influence  the 
gage  outputs,  and  second,  to  examine  the  form  of  the  wave  transmitted  through 
the  PMMA.  In  the  CS^  shots,  the  CS^-PMtlA  gages  exhibited  initial  shock  wave¬ 
forms  which  are  rounded,  a  characteristic  of  viscoelastic  materials.  3y  the 
time  the  shock  reaches  the  middle  CS^  gage,  all  the  rounding  is  essentially 
gone.  This  change  was  attributed  to  shocking  up  in  the  liquid  but  a  check 
was  desired. 

No  extraneous  electrical  signals  were  generated  and  rounded  wave¬ 
forms  were  recorded  by  each  gage,  even  when  the  wave  reflected  from  the 
sapphire  back  plate  at  a  relatively  high  pressure.  This  substantiates 
the  shocking  up  idea.  Hugoniot  data  generated  by  this  experiment  agree 
well  with  other  data. 


H20  Shot 

In  order  to  make  sure  the  reaction  data  obtained  in  the  CS^  shots 
were  not  caused  by  the  experimental  setup,  one  cell  was  filled  with  distil¬ 
led  water,  which  was  felt  to  be  a  nonreacting  material  (although  there  was 


some  possibility  a  phase  change  might  occur) , 


No  anomalous  gage  response  or  two-wave  structure  was  observed.  A 


Hugonlot  point  was  generated  that  is  consistent  with  other  measured  data. 
This,  along  with  the  PMMA  shot  results,  substantiate  the  utility  of  using 
this  method  for  obtaining  good  Hugoniot  data  that  does  not  depend  on  other 
parameters  such  as  projectile  velocity  or  interface  material  Hugoniot. 

Sapphire  Gage  Shots 

Two  shots  (77-094  and  77-095)  were  made  in  which  an  electromagnetic 
gage  was  sandwiched  between  two  pieces  of  sapphire,  which  were  then  used  for 
Che  back  place  of  a  CS^  cell.  The  purpose  for  doing  this  was  twofold:  to 
measure  the  particle  velocity  decay  wave  transmitted  into  the  sapphire  with¬ 
out  the  noise  superimposed  on  it,  and  to  see  if  the  noise  was  an  electrical 
signal  propagating  through  the  whole  region  rather  than  just  noise  at  the 
CS^-sapphire  gage. 

Waveforms  from  the  shots  indicate  the  noise  is  an  electrical  signal 
propagating  through  the  whole  experiment,  starting,  apparently,  at  the  time 
reaction  starts.  In  fact,  on  these  and  most  later  shots,  the  noise  was 
detected  on  all  the  signals  at  the  same  time  although  the  perturbations  are 
quite  small  on  the  higher  level  signals. 

Obtaining  a  noise-free  record  of  the  particle  velocity  decay  was  not 
possible  because  the  noise  was  still  being  recorded  at  the  time  the  shock 
arrived  at  the  gage  approximately  0.15  usee  later. 

Rarefaction  Speed  Measurements 

An  unexpected  piece  of  data  was  measured  as  a  result  of  the  sapphire 
inpactor  and  PMMA  cell  front  being  relatively  thin  in  the  early  experiments. 
A  rarefaction  from  the  interfaces  traveled  through  the  cell  during  the  time 
the  CS2  was  at  the  Initial  pressure,  P^.  By  constructing  a  time-distance 


diagram  it  was  possible  eo  determine  the  speed  of  Che  beginning  of  the  rar¬ 
efaction  (which  is  simply  the  sound  speed  at  state  1)  on  three  separate 
experiments.  The  data  for  these  shots  are  shown  in  Table  4.9. 


TABLE  4.9.— CS ^  Sound  Speed  Measurements 


Shot  No. 

h 

(kbars) 

State  1  Data 

u.i 

(mm/usec) 

V1 

(cm3/g) 

Sound  Speed 

C1 

(mm/usec) 

77-003 

30.8 

2.81 

0.548 

3.03 

1.08 

77-044 

43.4 

3.22 

0.530 

3.62 

1.12 

77-055 

20.5 

2.45 

0.579 

2.73 

1.11 

Sound  speeds  obtained  in  this  manner  can  be  used  to  calculate  the 
isentropic  bulk  modulus  and  also  the  Gruneisen  parameter  at  state  1  condi¬ 
tions.  The  isentropic  bulk  modulus  can  be  calculated  from 

Ks  V^3V^S  "  T  *  (4.2) 

The  Gruneisen  parameter  can  be  obtained  by  solving  the  P-V  differential 
equation  of  the  Hugoniot  for  T  to  obtain 


r 


2[v2(— >  +  c2i 

1  vd\r  H  J 

v>(§)h<Vv>  +  (p-V> 


(4.3) 


Values  for  the  slope  of  the  P-V  Hugoniot  were  calculated  using  the  unreacted 
CS^  EOS.  Results  of  this  analysis  are  listed  in  Table  4.10. 

Computer  Simulation  of  Experiments 

One-dimensional  wave  propagation  code  calculations  using  the  POD 
code  available  at  WSU  were  made  to  determine  if  the  waveforms  recorded  by 


Che  particle  velocity  gages  could  be  reproduced,  thereby  leading  to  a  better 

understanding  of  the  process.  POD  is  a  general  purpose  code  using  a  second 

order  finite  difference  numerical  integration  of  the  flow  equations.  It 

30 

utilizes  the  von  Neumann-Richtmeyer  artificial  viscosity  to  prevent  numer¬ 
ical  instabilities  at  the  shock  front. 


TABLE  4. 10. --Calculation  of  Kg  and  T  from  Sound  Speed 


Shot  No. 

<dP/dV)„“ 

H  2 

(cm/ysec) “ 

Ks 

(mbars) 

r 

r/v 

(g/cm3) 

77-003 

-0.460 

0.168 

2.08 

3.79 

77-044 

-0.675 

0.247 

1.66 

3.13 

77-055 

-0.296 

0.129 

2.01 

3.47 

Calculated  using  unreacted  CS^  EOS. 


Andrews  5  wrote  an  equation  of  state  subroutine  which  allows  a  mate¬ 
rial  to  react  according  to  a  time  dependent  reaction  rate  relation.  This 
subroutine  requires  two  ocher  subroutines  be  written  to  calculate  the  equa¬ 
tion  of  state  properties  of  the  reactant  and  products.  The  equations  of 
state  described  in  Section  4.2  were  used  for  this  purpose. 


Input  Data 

To  simulate  the  experiments,  it  was  necessary  to  provide  mechanical 
P-V  equations  of  state  for  PMMA  and  sapphire.  The  form  of  the  P-V  relation 
is 


A  +  By  +  CU2  +  DU3 


(4.4) 


where  U  ■  -rr  -  1.  Parameters  used  for  PMMA  and  sapphire  are  listed  in 


Table  4.11.  PMMA  has  an  unnatural  relation  because  there  was  no  desire  to 


model  its  viscoelastic  behavior  but  rather  to  provide  a  form  that  gives  the 
proper  state  in  CS^  for  a  given  projectile  velocity. 


TABLE  4.11. — PMMA  and  Sapphire  EOS  Constants  for  POD 


A 

B 

C 

D 

0o, 

Material 

(mbars) 

(mbars) 

(mbars) 

(mbars) 

(g/cm3) 

PMMA 

0.0 

0.047 

-0.1388 

0.5525 

1.185 

Sapphire 

0.0 

4.99 

4.637 

17.57 

3.985 

The  equation  of  state  parameters  used  for  unreacted  CS£  and  the 
reacted  products  are  those  of  Tables  2.2  and  2.3,  respectively. 

Shots  Modeled 

Calculations  were  made  to  simulate  a  number  of  shots  ranging  from  no 
reaction  in  the  reflected  wave  to  reaction  in  the  initial  wave.  The  following 
shots  were  modeled:  77-055,  77-103,  77-107,  77-071,  77-106,  and  77-045. 


Reaction  Rate  Law  Used 

A  two-step  rate  law  based  on  the  complex  Arrhenius  rate  law  of 
Table  2.5  was  used  in  the  calculations.  The  rate  law  below  1125°K  was 


dx 

dt 


1.3  x  10 


15 


-72,500/RT  ,  -1 

e  psec 


(4.4) 


and  above  was 


dx 

dt 


1300 


,-10,400/RT  UMC-1 


(4.5) 


A  maximum  rate  was  established  by  limiting  the  change  in  x  to  0.05 
for  each  time  step  in  order  to  maintain  stability.  This  was  not  a  serious 
constraint  because  rates  up  to  50  usee**1,  were  still  possible. 
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Results  of  Calculations 


Calculated  waveforms  did  not  simulate  the  measured  waveforms  very  well 
because  a  great  deal  of  relaxation  was  calculated  behind  the  first  wave  that 
was  not  observed  in  the  experiments.  Time  limitations  made  it  impossible  to 
determine  if  this  relaxation  was  an  artifact  of  the  calculations  owing  to 
interface  instabilities  or  if  other  undiscovered  problems  were  involved.  For 
this  reason,  none  of  the  wave  structure  results  will  be  presented  here.  It 
is  presently  felt  that  the  wave  structure  between  the  R^  and  R^  waves  may 
hold  some  very  important  information  regarding  intermediate  reactions,  etc. 
that  cannot  yet  be  interpreted  because  of  lack  of  understanding. 

There  are  several  important  pieces  of  information  that  result  from 
the  calculations  which  are  rather  insensitive  to  the  wave  structure.  Temper¬ 
ature  and  volume  changes  resulting  from  the  reaction  are  available  from  the 
output  of  the  code  calculations.  Temperature  change  information  is  summarized 
in  Table  4.12  and  volume  change  information  in  Table  4.13. 


TABLE  4.12.— Temperature  Change  Information  from  POD  Calculations 


Simulated 

Shot  No. 

T1 

(QK) 

T3 ' 

no 

T3 

(°K) 

Temperature 

Change 

<VV 

("K) 

77-055 

635 

879 

873 

no  reaction 

77-103 

766 

1079 

1650 

571 

77-107 

857 

1187 

1720 

533 

77-071 

808 

1139 

1680 

541 

77-106 

960 

1295 

1770 

475 

77-045 

1000 

1521 

1820 

299* 

^Reaction  in  the  initial  shock  was  apparent. 


TABLE  4.13.— Volume  Change  Information  from  POD  Calculations 


Simulated 
Shot  Mo. 

Y 

(cnr/g) 

v3 

(cm3/g) 

Volume 

Change 

(V3.-V3) 

(cm^/g) 

Percent 

Change 

77-055 

0.505 

0.503 

no  reaction 

— 

77-103 

0.480 

0.392 

0.088 

18.3 

77-107 

0.470 

0.381 

0.089 

18.9 

77-071 

0.474 

0.387 

0.087 

18.4 

77-106 

0.461 

0.372 

0.089 

19.3 

77-045 

0.447 

0.364 

0.083 

18.6 
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4.5  DISCUSSION  OF  RESULTS 


This  chapter  comprises  a  discussion  of  the  results  of  the  previous 
chapter  in  comparison  with  other  published  data  and  with  various  ideas 
developed  in  Sect.  4.2.  Discussion  will  include  applicability  of  the  equa¬ 
tions  of  state  developed  for  the  unreacted  CS^  and  reacted  products,  reaction 
rate  analysis,  enlightenment  from  code  calculations,  and  speculation  regarding 
the  noise  and  reaction  mechanism. 


Unreacted  CS^  Hugoniot  Data 

Fourteen  Hugoniot  points  for  the  unreacted  CS2  principal  Hugonict  have 

been  generated  in  this  study  and  are  plotted  along  with  data  of  Cook  and 

31  32  14  33 

Rogers,  Lysne,  Dick,  and  Walsh  and  Rice  in  U  -u  and  P-V  spac*  i  • 

s  p 

Fig.  5.1.  Also  plotted  is  the  calculated  Hugoniot  based  on  the  unrea^ed 
CS2  EOS.  It  is  apparent  that  the  data  generated  in  this  study  are  in  agree¬ 
ment  with  the  results  of  ocher  experiments. 

Dick’s  data  reported  in  Ref.  34  were  based  on  a  1968  revision  of  the 

Hugoniot  used  to  determine  the  states  in  the  aluminum  buffer  place  in  front 
of  the  CS2>  The  effect  of  changing  Che  aluminum  Hugoniot  from  its  earlier 
form  is  indicated  by  six  points  in  Fig.  5.1  with  arrows  from  three  original 
points  to  the  later  data. 

Dick  estimated  the  cusp  in  the  Hugoniot  to  occur  at  approximately 
62  kbars.  The  data  plotted  in  Fig.  5.1  suggest  that  large  deviations  from 
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Fig.  5.1. — Comparison  of  Unreacted  CS2  llugoniot  Data.  Lines  are  unreac 
(a)  is  the  Ua-un  llugoniot  data. 


Fig.  5.1.— (Continued),  (b)  is  tha  P-V  Hugoniot  Data. 


the  calculated  curve  begin  around  50  kbars.  Shot  77-045  had  an  anomalous 


CS2~sapphire  gage  response  that  appears  to  have  been  caused  by  a  two-wave 
structure  indicating  that  reaction  occurred  in  the  initial  shock.  This  shot 
had  state  1  pressure  of  48.3  kbars  (which  would  be  in  error  if  the  wave  were 
reacting).  Shot  77-044,  with  a  state  1  pressure  of  43.8  kbars,  shows  no 
evidence  of  reaction.  Based  on  this  information  it  is  estimated  that  the 
cusp  occurs  in  the  46  to  54  kbar  range.  This  would  also  seem  to  indicate 
that  a  single  shock  study  of  the  cusp  is  possible  on  the  WSU  gas  gun. 


Unreacted  CS^  EOS 

The  calculated  equation  of  state  fits  the  data  very  well  up  to  about 
45  kbars  at  which  point  the  data  seems  to  fall  off  in  the  direction  cor¬ 
responding  to  a  transition. 

Sound  speed  data  lead  to  calculations  of  isentropic  bulk  modulus  and 
Gruneisen  parameter.  These  data  have  been  plotted  along  with  the  curve  that 
represents  the  unreacted  CS^  EOS  model  in  Fig.  5.2.  The  bulk  modulus  data 
lie  very  close  to  the  model  curve.  T/V  is  in  the  same  region  as  the  model 
curve  but  appears  to  have  a  good  deal  of  scatter  or  else  to  be  a  very 
complicated  function. 

These  data,  along  with  the  good  fit  to  the  Hugoniot  data,  suggest 
that  the  equation  of  state  is  quite  good  in  the  principal  Hugoniot  region. 

One  can  infer  from  this  that  the  temperature  calculations  are  reasonable 
although  Kg  and  r/V  are  not  the  best  indicators.  Accurate  temperature  mea¬ 
surements  are  still  needed. 


Calculated 


Pressure,  kbars 
(b)  Ks-P 

x-Data  based  on  sound 
speed  measurements 


»  C.l 


Calculated 


0  10  20  30  40  50 

Pressure,  kbars 

Fig.  5.2.— Comparison  of  T/V  and  K,  Valuta.  Solid 
lints  rtprtstnc  Che  calculated  values  based  on  the  models 
used  in  the  unreacced  CS2  EOS.  Individual  x's  ere  values 
calculated  from  sound  3peed  measurements.  Plots  are  made 
versus  state  1  pressures. 
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Unreacted  CS»  Reshock  Data 

A  comparison  has  been  made  between  the  reshock  data  calculated  from 
the  unreacted  CS^  EOS  and  measured  values.  Two  comparisons  can  be  made,  one 
at  state  3'  (the  initial  reshock  state)  and  the  other  at  state  2  after  pas¬ 
sage  of  the  R^  wave. 

Figure  5.3  is  a  comparison  of  state  3'  as  measured  by  the  CS^-sapphire 
gage  with  calculated  values  for  the  six  shots  where  data  were  available 
Note  that  three  sets  of  points  are  essentially  on  top  of  each  other 
and  the  calculated  values  are  below  the  measured  values  in  the  other 
three  sets. 

Since  the  calculated  points  are  based  on  an  Initial  shock  pressure, 

P,  ,  several  possibilities  for  these  differences  exist.  Either  the  EOS  is 
wrong  or  the  state  1  pressure  is  wrong  or  the  state  3r  particle  velocity  at 
the  CS2*sapphire  interface  is  wrong.  These  six  points  have  been  flagged  in 
Fig.  5.1.  Notice  that  the  shots  showing  good  agreement  (77-055,  77-062,  and 
77-071)  all  lie  close  to  the  calculated  Hugoniot  line,  indicating  the  EOS 
gives  accurate  reahock  states.  Shots  77-103  and  77-050  are  also  close  to 
the  calculated  Hugoniot  so  the  difference  is  attributed  to  a  high  u^,  ®«a~ 
surement.  Shot  77-094  does  not  lie  very  close  to  the  calculated  Hugoniot  so 
its  difference  is  probably  caused  by  a  relatively  Inaccurate  state  1  measure¬ 
ment. 

The  state  3'  measurements  provide  an  independent  check  on  both  state  1 
measurements  and  the  calculated  reshock  Hugoniot.  The  agreement  above  indi¬ 
cates  that  both  measurements  and  calculations  are  reasonably  accurate. 

The  other  comparison  to  be  made  of  measurements  to  reshock  calcu¬ 
lations  is  at  state  2,  which  occurs  at  the  top  of  the  wave.  Figure  5.4 
provides  a  comparison  of  calculated  and  measured  P-V  values  for  state  2  from 
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Measured 

Calculated 


Fig.  5.3. — Comparison  of  Calculated  and  Measured  State  3'  Information.  Measured  values  were 
obtained  from  the  CS2~sapphiro  gage  before  reaction  was  observed.  Calculated  values  were  based  on 
the  pressure  at  state  1  and  obtained  graphically  by  determining  the  point  where  the  left  gotng  resliock 
ilugoniot  contacted  the  right  going  sapphire  llugoniot  centered  at  zero  pressure. 
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Fig.  5.4.-- ’Comparison  of  Calculated  and  Measured 
State  2  Information.  Measured  points  were  determined  from 
measured  R^  wave  data.  Calculated  points  were  obtained  by 
making  reshock  calculations  centered  at  Pj^  based  on  the 
unreacted  CS2  EOS  with  P2  (the  same  as  for  the  measured 
values)  as  the  final  pressure. 
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Table  4.4.  As  can  be  seen,  all  the  data  are  to  the  left  of  the  principal 
Hugonlot,  which  is  as  it  should  be  but  the  measured  states  are  to  the  left 
of  the  calculated  states  with  difference  varying  widely.  This  seems  to 
Indicate  the  measured  data  are  less  accurate,  possibly  being  influenced  by 
the  use  of  the  particle  velocities  several  times  in  the  jump  conditions  to 
obtain  the  states. 

It  should  also  be  remembered  that  the  R^  wave  is  evolving  as  can  be 
seen  from  Table  4.7.  Values  of  P  and  V  in  state  2  have  been  determined  from 
the  jump  conditions,  which  do  not  strictly  apply  to  a  changing  wave.  This 
undoubtedly  has  also  led  to  some  inaccuracies. 


Temperature  Calculations  Using  Unreacted  CS^  EOS 

Calculated  temperatures  have  been  plotted  in  temperature-pressure 

35 

space  in  Fig.  5.5  with  the  T-P  diagram  developed  by  Butcher  et  al.  shown  for 

comparison.  Figure  5.5a  shows  the  principal  Hugonlot,  the  reshock  Hugonlots 

for  four  shots  and  a  line  representing  the  locus  of  possible  state  3'  reshock 

states.  Notice  that  the  principal  Hugonlot  above  20  kbars  and  all  the 

reshock  end  states  are  well  up  into  the  region  labeled  "decomposition  to 

35 

carbon  and  sulfur"  by  Butcher  et  al.  This  would  seem  to  be  strong  evidence 
that  the  reaction  is  a  decomposition  reaction  rather  than  a  polymerization 
to  the  black  CS^  solid  as  indicated  by  Dick  and  Kusubov.^ 

Figure  5.5b  is  part  of  the  same  plot  with  the  measured  R^  wave  states 
plotted.  The  region  of  the  cusp  has  also  been  indicated  as  being  between 
46  and  54  kbars  as  discussed  previously.  Two  regions  have  been  shaded  in 
the  figure.  If  the  transition  were  strictly  pressure  determined,  the  R^ 
wave  states  would  be  expected  to  lie  in  the  vertical  region.  If  it  were 
temperature  dependent,  the  states  would  be  expected  to  be  in  the  horizontal 
region,  as  most  of  them  are.  Only  the  data,  for  the  low  pressure  nonreactive 
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Fig.  5.6.— Reacted  Produce  Hugonioc  Information.  The  lines 
are  calculated  Hugoniots  based  on  the  reacted  product  COS,  one  cen¬ 
tered  at  293*K  and  the  other  representing  states  achieved  after  reac¬ 
tion.  Data  points  obtained  in  this  study  were  obtained  from  measured 
R  wave  data  (state  3)  in  a  reshock  situation.  Data  of  other  experi¬ 
menters  were  obtained  in  single  shock  experiments. 


shoes  and  chat  for  shoe  77-045  (which  has  previously  been  mentioned  as 
probably  reacting  in  the  initial  wave)  are  outside  this  region.  This  indi¬ 
cates  the  reaction  is  more  closely  associated  with  temperature  than  pressure 
and  supports  the  hypothesis  that  a  chemical  reaction  occurs  in  the  shock. 

This  analysis  is  based  on  steady  waves  and  the  waves  in  the  experiments  are 
evolving  so  it  is  impossible  to  say  that  the  R^  wave  states  never  would  evolve 
into  the  vertical  region.  However,  a  review  of  Table  4.7,  along  with  the 
various  waveforms,  indicates  chat  the  greater  portion  of  the  evolution  has 
occurred  in  the  high  pressure  shots  by  the  time  the  R^  wave  passes  the  inter¬ 
mediate  gages. 

The  scatter  of  data  is  somewhat  bothersome  since  it  had  been  hoped 
that  the  R^  wave  temperatures  would  lie  on  a  single  line,  but  this  may  be 
due  to  the  fact  that  the  steady  wave  structure  has  not  had  time  to  establish 
itself  in  the  approximately  1  mm  of  travel  available  in  the  CS^  cell.  It 
would  be  very  helpful  if  experiments  with  thicker  CS^  cavities  could  be  done 
to  better  determine  the  steady  wave  profiles.  The  inaccuracies  discussed 
earlier  in  the  R^  wave  states  should  also  be  kept  in  mind. 

Reacted  Product  Data  , 

Comparison  of  the  measured  state  3  data  to  chat  published  and  also 
to  the  reacted  product  EOS  will  be  made  in  this  section.  In  addition,  com¬ 
parison  of  the  R2  wave  data  to  the  relaxed  CS^-sapphire  gage  data  will  be 
made  and  then  a  discussion  concerning  the  nature  of  recovered  reacted  pro¬ 
ducts  will  be  presented. 

Reacted  Product  EOS  Comparison 

Figure  5.6  is  a  plot  of  all  the  reacted  prodvc  Hugoniot  data  avail- 
14 

able  in  the  literature  along  with  the  data  generated  in  this  study.  The 


curve  to  the  left  of  the  data  is  the  Hugoniot  for  the  reacted  products 


centered  at  an  initial  temperature  of  293°K  and  one  atmosphere  pressure. 

The  relative  position  of  the  two  curves  is  as  one  might  expect  because  of 
the  high  temperature  along  the  Hugoniot  at  reaction  temperature. 

The  amount  of  scatter  in  all  the  data  is  considerable  but,  neverthe¬ 
less,  the  points  cluster  about  the  calculated  curve.  The  considerable  scatter 
in  Che  R?  wave  data  of  this  study  may  be  the  results  of  particle  velocity 
error  multiplication  and  unsteady  waves  as  has  been  previously  discussed. 

Measured  State  3  Data  Comparison 

State  3  data  were  also  obtained  from  the  relaxed  particle  velocity 
recorded  by  the  CS^-sapphire  gage.  Since  this  is  a  direct  measurement,  one 
would  expect  ic  to  be  more  accurate  than  the  wave  data  obtained  from 
measured  shock  and  particle  velocity  in  conjunction  with  the  R^  wave  data. 
Comparisons  between  these  two  data  sets  are  made  in  Fig.  5.7  for  the  eight 
shots  where  this  comparison  was  possible.  As  can  be  seen  the  R,  wave  data 
are  sometimes  higher  and  sometimes  lower  than  the  other  data.  The  two  points 
are  quite  close  together  on  shots  77-107  and  77-095  and  one  might  infer  that 
the  P-V  states  from  R£  wave  data  are  accurate.  These  points  have  been 
singled  out  on  Fig.  5.6  and  are  closer  chan  most  to  the  calculated  Hugoniot, 
indicating  Chat  the  calculation  is  reasonable. 

Nature  of  Recovered  Reacted  Products 

Calculations  were  made  to  estimate  the  temperatures  one  might  expect 
by  shocking  293°K  reacted  products.  Shocking  reacted  products  to  100  and 
160  kbars  would  be  expected  to  cause  temperatures  to  increase  to  600  and 
1000*K,  respectively.  For  reshocked  unreacted  CS^  with  a  100  kbar  pressure 
at  state  3',  the  temperature  is  about  1250“K.  Add  to  this  another  500  degrees 
for  the  reaction  input  and  the  reacted  product  temperature  would  be  about 
1700°K.  Upon  releasing  the  pressure,  the  residual  temperature  of  the  reacted 


produces  would  be  on  the  order  of  1000°K.  All  the  known  polymers  (black  CS^, 


polymerized  CS  and  C.^  polymer)  decompose  well  below  this  temperature 
(450°K,^  470°K,^  and  360°K^‘‘).  It  seems  reasonable  that  one  should  recover 
only  carbon  and  sulfur  in  recovery  experiments. 

The  temperatures  are  not  to  be  considered  extremely  accurate  but  they 
do  serve  to  point  out  that  residual  temperatures  are  high  and  would  be  expec¬ 
ted  to  greatly  influence  the  nature  of  the  recovered  material. 

Reaction  Rate  Data 

Five  different  kinds  of  data  have  been  analyzed  for  the  purpose  of 
obtaining  information  about  the  rate  law  which  governs  the  reaction  process. 
Induction  time,  particle  velocity  decay  and  noise  decay  data  were  all  obtained 
from  the  gage  record  at  the  CS^-sapphire  interface.  P^  wave  decay  and  R.,  wave 
structure  data  were  obtained  from  an  analysis  of  the  time-correlated  wave 
forms.  These  will  be  discussed  in  the  following  subsections. 

R^  Wave  Decay  Calculations 

wave  decay  data  were  reported  in  terms  of  particle  velocity  decay 
in  Chapter  4  but  were  not  reduced  to  initial  reaction  rate  data.  The  first 
term  in  the  numerator  of  Eq.  (2.40)  is  negligible  in  these  experiments  so  the 
equation 
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was  used  to  calculate  x.  Calculations  were  made  using  Eq.  (5.1)  and  the 
calculated  conditions  of  state  3’,  along  with  the  particle  velocity  derivative 


values  from  Table  4.7.  Table  5.1  contains  the  results  of  this  calculation 

2 

along  with  estimates  of  AE,  AV,  and  K^V/U  which  were  used. 

Reaction  Rate  Information  Analysis 

The  various  reaction  rate  values  are  summarized  in  Table  5.2  for  each 
shot.  The  temperatures  are  those  obtained  from  reshock  calculations  at 
state  3’.  This  is  probably  not  an  accurate  temperature  during  the  reaction 
because  energy  is  released  and  the  temperature  increases  several  hundred 
degrees  with  the  reaction  rate  likewise  increasing.  In  the  absence  of  more 
meaningful  data  it  was  used  so  that  reaction  rate  versus  temperature  plots 
could  be  made. 

All  the  points  were  plotted  on  Fig.  5.9  for  comparison  purposes.  The 
particle  velocity  decay  shows  the  highest  reaction  rate  while  the  noise  decay 
the  lowest.  The  R^  wave  decay  and  R^  wave  structure  points  generally  lie 
between  the  two.  Notice  that  the  particle  velocity  decay  and  the  &2  wave 
structure  data  appear  to  indicate  two  different  rates,  one  below  1125°K  and 
one  above. 

Induction  and  Relaxation  Time  Analysis 

Table  5.2  also  contains  information  for  the  induction  time  and  par¬ 
ticle  velocity  and  noise  relaxation  times.  These  have  been  plotted  as  the 
logarithm  of  the  time  versus  reciprocal  temperature  in  Fig.  5.10.  Note  that 
the  induction  time  points,  along  with  the  particle  velocity  points,  indicate 
a  highly  temperature  dependent  low  temperature  rate  much  the  same  as  the 
reaction  rate  data. 
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TABLE  5.1. — Calculation*  of  Reaction  Rate  from  R,  Wave  Information 
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Fig.  5.8. — Arrehenius  Plots  of  the  Reaction  Rate  Data,  (a)  is 
from  particle  velocity  decay,  (b)  is  from  noise  decay.  Temperature  is 
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Fig.  5. 9. —Comparison  of  Arrhenius  Plots.  All  four  estimates 
of  the  reaction  rate  are  shovn.  Note  chat  a  low  temperature,  as  well 
as  a  high  temperature,  rate  lav  is  indicated. 
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Fig.  5.10. — Induction  Time  and  Relaxation  Time  Plots,  tj  is 
the  induction  time,  tUpis  the  particle  velocity  relaxation  time  and 
t>j  is  the  noise  relaxation  time.  Temperatures  are  at  state  3'  (T„ , ) 


lace  Lav  Parameters 


Activation  energy  (E)  and  frequency  factor  (A)  can  be  obtained  from 
the  various  plots  by  assuming  a  zeroth  order,  complex  Arrhenius  rate  law  of 
the  form 
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(5.2) 


Relaxation  time  data  can  also  yield  estimates  of  these  parameters  as  dis¬ 
cussed  in  Sec.  4.2. 

Values  for  these  parameters  obtained  from  the  various  plots  are 
tabulated  in  Table  5.3.  Obviously,  there  are  considerable  differences  in  the 
values  as  one  might  expect  from  the  way  the  rates  were  estimated.  They  do 
indicate,  however,  the  difference  between  the  two  rate  laws  and  also  the 
magnitude  of  the  parameters  one  might  expect. 

If  one  uses  the  particle  velocity  decay  time  constants  to  compare  the 
behavior  of  the  initial  wave  in  shots  77-045  and  77-106  as  to  reaction  rate, 
an  interesting  correspondence  to  the  measured  data  results.  Table  5.4  is  a 
listing  of  the  values  used  and  the  calculated  reaction  rate  in  the  initial 
wave.  It  indicates  that  evidence  of  reaction  should  have  been  observable  in 
77-045  and  not  in  77-106.  As  indicated  earlier,  77-045  had  a  CS^-sapphire 
gage  record  that  only  seems  explicable  by  the  initial  wave  having  a  two-wave 
structure  as  a  result  of  reaction.  Shot  77-106  showed  none  of  this  anomalous 
behavior.  This  would  seem  to  Independently  corroborate  the  low  temperature 
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TABLE  5.4.- 

-Reaction 
Shock  in 

Rate  Calculations 
Shots  77-045  and 

for  the  Initial 
77-106 

Shot  No 

P1 

.  kbars 

T1 

(°K) 

E 

(cal/mole) 

A 

(usee 

• 

X 

(Usee  1) 

77-045 

48.3 

1058 

72,500 

1.3  x  1015 

1.37 

77-106 

38.5 

908 

72,500 

1.3  x  1015 

0.005 

All  the  sets  of  rate  law  parameters  are  unphysical  and  therefore 
this  whole  analysis  must  be  considered  highly  speculative.  It  does  serve 
to  point  out  that  a  simple  overall  rate  law  which  describes  the  reaction 
accurately  is  not  possible.  It  also  indicates  that  many  experiments  (as 
is  normally  the  case  in  reaction  kinetic  experiments)  must  be  carried  out 
before  useful  and  believable  reaction  rate  laws  can  be  established. 

Comments  Regarding  Induction  Time 

Induction  times  have  been  observed  in  shock  tube  decompositions  of 
37  38 

gaseous  CS^  and  SO^.  These  observations  have  been  explained  as  the 
time  required  to  reach  a  critical  concentration  of  a  particular  species. 

In  view  of  this,  it  seems  probable  that  the  observation  of  an  induction 
time  in  the  experiments  in  this  study  is  a  strong  indicator  of  a  multistep 


reaction  mechanism. 


Electrical  Noise 


As  has  been  briefly  touched  on  in  Chapter  4,  electrical  noise  was 
recorded  during  the  experiment.  The  following  discussion  relates  what  is 
presently  known,  along  with  some  conjecture  concerning  the  nature  of  it.  It 
is  possible  that,  once  understood,  this  noise  phenomenon  might  become  a 
valuable  experimental  tool  in  diagnosis  of  certain  chemical  events. 

Gage  Signals  with  Noise 

The  noise  was  best  measured  by  the  CS2~sapphire  gage  because  of  the 
low  voltage  output  of  this  gage  and  greater  sensitivity  to  electrical  noise. 
Two  typical  records  are  shown  in  Fig.  5.11:  shot  77-070  which  showed  no 
induction  time  and  shot  77-071  which  did  show  one.  The  noise  is  large  at 
first,  decays  with  time  and  has  a  frequency  of  about  150  MHz. 

Noise  Appearance 

Noise  appearance  corresponds  with  the  reaction  as  evidenced  by  the 
following  observations. 

1.  Particle  velocity  decay  appears  at  the  same  time  as 
noise  decay. 

2.  A  two-wave  structure  is  observed  in  the  same  experiments 
that  noise  is  observed. 

3.  The  reaction  rate  law  based  on  noise  decay  time  is  con¬ 
sistent  with  other  rate  laws. 

4.  Noise  was  not  observed  in  the  PMMA  experiment  even  though 
the  PMMA  reshock  pressure  was  in  the  same  range  as  in  four 
CS2  experiments  where  noise  had  been  observed. 
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Fig.  5.11.—  Selected  CS2*-Sapphire  Gage  Data  Showing  Noise. 
Waveform  in  (a)  shows  an  induction  time.  Both  were  recorded  at  sweep 
rate  of  50  nanosec/div. 


V 


Nature  of  the  Noise 


The  nature  and  cause  of  the  electrical  signal  remains  a  mystery 
although  a  little  information  has  been  generated.  As  mentioned  earlier, 
a  gage  sandwiched  between  two  pieces  of  sapphire  recorded  the  noise  at  the 
same  time  as  the  CS ,,-sapphire  gage,  long  before  the  shock  arrived.  The  noise 
was  also  observed  as  a  small  perturbation  on  the  other  gage  outputs  at  a  cor¬ 
responding  time.  This  indicates  the  noise  is  an  electrical  signal  that  per¬ 
vades  the  whole  region. 

It  has  been  suggested  that  the  noise  may  be  an  interaction  between 
ions  formed  by  the  decomposition  and  the  magnetic  field.  A  very  puzzling 
characteristic  is  that  the  signal  appears  to  hit  a  maximum  and  then  decay 
which  suggests  a  one-time  occurrence.  If  the  noise  were  a  result  of  the 
continuing  reaction,  it  should  have  been  observed  during  the  whole  experiment 

Another  explanation  appears  to  be  that  of  a  runaway  first  reaction 
(thermal  explosion)  that  results  in  a  large  number  of  intermediate  species 
forming  very  rapidly  at  the  sapphire-CS2  interface.  Then  as  the  reaction 
settles  down  to  a  more  or  less  steady  process,  the  concentration  of  these 
intermediate  species  drops  rapidly  to  a  much  lower  steady  state  value  in  the 
reacting  wave.  The  one-time  occurrence  of  the  noise  appears  to  be  another 
indicator  of  a  multiple  step  reaction  mechanism. 

Discussion  of  POD  Calculations 

Although  it  was  not  possible  to  accurately  simulate  the  measured 
waveforms,  it  was  possible  to  obtain  information  regarding  the  temperature 
and  volume  change  resulting  from  the  reaction.  As  indicated  in  Table  4.11,  a 
temperature  change  of  about  570®K  was  calculated  for  shot  77-103  which  had  a 
temperature  of  1080°K  at  state  3'.  The  temperature  change  for  shot  77-106 
with  a  state  3'  temperature  of  1300°K  was  about  470°K.  Shots  77-107  and 
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77-071  were  between  these  two  values.  This  indicates  a  temperature  increase 
of  about  500°K  should  be  expected  as  a  result  of  the  reaction.  It  is 


interesting  to  note  that  C^S^  static  high  pressure  experiments  showed  a 

39 

temperature  change  of  about  330°K  when  decomposition  was  observed. 

Table  4.12  indicates  a  volume  change  on  the  order  of  18.5”  results 
from  the  reaction.  R.  Dick  ^  estimated  11’/,  based  on  his  experiments. 

Calculations  to  simulate  shots  77-045  and  77-106  showed  reaction  in 
the  initial  wave  of  the  former  but  not  the  latter,  as  was  discussed  in 
Sec.  5.3.4.  Calculated  particle  velocities  at  the  CS^-sapphire  interface  for 
shot  77-045  were  much  like  those  measured  but  stability  problems  obscured  the 
exact  nature  of  the  calculated  waveform,  making  it  impossible  to  make  a 
definitive  statement. 

Further  details  of  this  work  are  contained  in  Ref.  40. 
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4.6  SUMMARY 


Unreacted  CS^ 

Fourteen  unreacted  CS^  Hugoniot  points  were  generated  in  this  study 
using  an  electromagnetic  gage  technique,  most  of  them  from  Initial  wave  data 
in  reshock  experiments.  Since  both  shock  velocity  and  particle  velocity 
were  measured,  these  points  are  independent  of  other  assumptions.  They  cor¬ 
relate  well  with  other  measured  data. 

Verification  of  the  reasonableness  of  the  unreacted  CS2  EOS,  which 
was  based  on  a  Murnaghan  isotherm  and  developed  forms  for  and  17 V,  came 
as  a  result  of  the  good  fit  to  the  Hugoniot  data  and  also  from  a  surprisingly 
good  comparison  with  Kg  and  I7V  calculated  from  sound  speed  measurements  on 
three  shots.  This  equation  of  state  was  fit  to  a  universal  liquid  Hugoniot 
in  order  to  get  the  correct  parameters,  so  the  agreement  with  the  Hugoniot 
measurements  also  demonstrates  the  applicability  of  the  universal  liquid 
Hugoniot  to  liquid  CS2- 

Comparing  the  calculated  unreacced  CS,,  EOS  Hugoniot  with  the  measured 
data  in  this  study  and  from  other  studies,  indicates  the  two  begin  to  deviate 
in  the  range  of  46  to  54  kbars.  One  shot  with  an  initial  wave  amplitude  of 
about  48  kbars  showed  definite  indications  of  reaction  occurring  in  the 
initial  shock.  Two  shots  at  38.5  kbar  and  one  at  43.4  kbar  do  not  show  any 
indications  of  initial  shock  reaction,  which  further  supports  the  existence 
of  a  cusp  in  the  46  to  54  kbar  region. 


Measured  reshock  states  compared  quite  well  with  calculated  states 
but  the  former  had  quite  a  lot  of  scatter,  ascribed  mainly  to  particle  velocity 


error  multiplication  in  the  jump  condition  calculations.  Reshock  temper¬ 
ature  calculations  showed  the  dramatic  effect  of  reshoclcing  on  the  state  of 
the  material,  with  reshock  temperatures  being  several  hundred  degrees  below 
single  shock  temperatures  at  the  same  pressure.  This  made  Lt  much  easier  to 
control  the  state  3'  temperatures. 

A  temperature-pressure  diagram  was  plotted  and  all  the  principal 

Hugoniot  states  above  20  kbais  are  in  the  region  which  has  been  design." t.’d  by 
35 

Butcher  et  al.  as  a  decomposition  region.  All  the  reshock  states  in  this 
study  were  also  in  that  region. 

It  was  observed  in  the  measurements  that  a  two-wave  structure  devel¬ 
oped  in  the  high  pressure  wave  reflecting  off  the  sapphire  back  place.  This 
two-wave  structure  development  corresponded  to  observations  at  the  CS7- 
sapphtre  gage  of  particle  velocity  decay  and  also  to  an  unexplained  electrical 
noise  signal,  both  of  which  were  delayed  by  an  induction  time  at  the  lower 
reshock  pressures. 

When  the  wave  states  were  plotted  on  the  temperature-pressure  dia¬ 
gram,  it  was  observed  that  they  fell  into  the  region  one  would  expect  for  a 
temperature  controlled  rather  than  a  pressure  controlled  reaction.  Some 
uncertainty  exists  in  this  statement  due  to  the  fact  that  the  R^  wave  was 
evolving  during  the  experiments  and  the  statement  is  strictly  true  only  for  a 
steady  wave.  Wave  forms  for  the  higher  pressure  shots  indicate,  however, 
that  evolution  was  nearly  complete  by  the  time  the  wave  reached  the  middle 
gage. 

Thermodynamic  calculations,  using  the  equations  of  state  developed 
for  unreacted  CS2  and  the  reacted  products,  indicate  that  CS2  is  metastable 
with  respect  to  a  mixture  of  carbon  and  sulfur  at  standard  conditions  as  well 
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as  ac  principal  Hugonioc  and  reshock  Hugoniot  states  attained  in  the  experi¬ 
ments.  It  seems  reasonable  that  the  shock-induced  temperature  rise  would 
cause  Che  reaction  rate  to  increase  to  such  an  extent  that  CS^  decomposition 
to  carbon  and  sulfur  could  very  well  occur  in  the  experiments. 

The  CS^  metastability ,  temperature  dependent  R^  wave,  induction  times 
and  noise  all  add  credibility  to  the  hypothesis  that  the  CS2  is  undergoing  a 
chemical  reaction  in  the  shock. 

Reacted  Product  Data  and  EOS 

Dick's  data  above  100  kbars  were  essentially  on  a  reacted  product 
Hugonioc,  so  they  were  used  to  fit  a  reacted  product  EOS  based  on  a  Murnaghan 
isotherm  and  a  constant  Cv  and  T/V.  It  is  interesting  that  the  specific 

3 

volume  was  treated  as  a  fitting  parameter  and  an  initial  value  of  0.48  cm  /g 

resulted.  This  correlates  quite  well  with  the  density  of  a  compacted  mixture 
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of  carbon  and  sulfur  at  0.46  to  0.48  cm  /g. 

R2  wave  states  were  calculated  by  substituting  the  measured  data  in 
the  appropriate  jump  conditions.  In  P-V  space  these  data  clustered  around 
the  calculated  EOS.  Relaxed  particle  velocity  data  from  the  CS^-sapphire 
gage  compares  quite  well  with  the  R^  wave  data  but  demonstrates  the  inac¬ 
curacies  in  the  R2  wave  data  calculations. 

Wave  propagation  code  calculations  indicate  the  volume  change  brought 
about  by  the  reaction  in  reshock  experiments  is  about  18.5%.  The  reaction- 
precipitated  temperature  change  depends  on  the  conditions  at  state  3’.  For  a 
T3,  of  1080*K,  the  change  is  570°K  while  at  1300°K  it  is  about  470aK. 

All  the  known  polymers  of  CS2,  CS,  and  C3S2  decompose  below  500aK, 
which  is  much  below  the  residual  temperatures  expected  for  the  reacted  pro¬ 
ducts.  The  high  residual  temperature  results  not  only  from  the  differences 


in  compress ibilicy  between  the  CS^  and  the  reacted  products  but  also  from 
the  approximately  500°K  temperature  increase  due  to  the  reaction.  This 
indicates  that  it  would  be  very  difficult  to  recover  anything  but  carbon  and 
sulfur  from  recovery  experiments  even  if  polymers  were  formed  in  the  shock. 

Reaction  Rate  Data 

Five  different  types  of  measurements  and  estimations  were  used  in  an 
effort  to  obtain  information  about  the  ’•eaction  rate.  They  were  induction 
time,  particle  velocity,  noise  decay,  wave  decay  and  R^  wave  structure 
data.  The  reaction  rate  data  ranged  from  2.5  per  Uscc  to  27  per  usee  depen¬ 
ding  on  the  particular  type  of  data  and  shot  being  examined. 

Plots  were  made  of  the  logarithm  of  either  reaction  rate  or  relax¬ 
ation  time  versus  reciprocal  temperature,  where  the  temperature  was  that  of 
state  3'.  The  particle  velocity  decay  gave  the  highest  reaction  rates  and 
the  noise  decay  the  lowest  reaction  rates,  with  the  other  two  in  beween. 

Two  different  rate  laws  were  suggested  by  the  induction  times,  the  particle 
velocity  decay  races  and  the  wave  structure  rates.  The  low  temperature 
rate  law  (below  1125°K)  was  highly  temperature  dependent  while  the  high 
temperature  rate  law  (above  1125°K)  varied  considerably  less  with  temperature 
but  nevertheless  gave  quite  high  reaction  rates. 

i 

Actual  rate  law  parameters  were  determined  for  the  zeroth  order,  com¬ 
plex  Arrhenius  rate  law.  The  low  temperature  rate  had  activation  energies 
from  37,200  to  72,500  cal/mole  and  frequency  factors  of  8.7  x  10^  to 
1.3  x  lO1^  usee'1  depending  on  the  particular  data  used.  The  high  temper¬ 
ature  rate  law  had  activation  energies  from  7000  to  13,200  cal/mole  with 
frequency  factors  of  102  to  2370  usee  Obviously  these  Arrhenius  parameters 
are  non-physical  and  the  reason  for  this  is  unknown.  It  should  be  remembered, 
however,  chat  all  the  race  law  estimates  were  made  without  knowing  the  true 
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relationship  between  the  observables  and  the  reaction  rate,  t.’hen  spectro- 
graphic  measurements  of  intermediate  species  are  made,  it  may  be  possible  to 
determine  which  of  the  shock  observables  yields  the  best  estimate  of  reaction 
rate. 

The  low  temperature  rate  law  was  corroborated  very  nicely  by  shot 
77-045  that  apparently  reacted  in  the  initial  wave.  According  to  the  particle 
velocity  decay  low  temperature  rate  law  it  would  have  had  a  reaction  rate  of 
about  1.4  usee  *  in  the  initial  wave,  so  in  the  approximately  0.5  sec  that 
it  traveled  from  the  front  of  the  cell  to  the  back,  a  substantial  amount  of 
reaction  would  have  taken  place.  Two  shots  at  slightly  lower  projectile 
velocities  showed  no  evidence  of  reaction  and  the  calculated  rate  was  0.005 
usee  \  so  they  would  not  be  expected  to  show  evidence  of  reaction. 

Induction  times  were  observed  on  some  of  the  lower  pressure  shots. 

37  38- 

Experimenters  studying  the  decomposition  of  gaseous  CS^  and  SO^  in  shock 
tubes  have  also  observed  induction  times  which  they  attributed  to  the  time 
required  for  a  particular  molecular  species  to  attain  a  critical  concen¬ 
tration. 

3oth  die  induction  time  and  the  two  reaction  rates  would  appear  to 
substantiate  the  hypothesis  of  a  chemical  reaction  occurring  and,  moreover, 
suggest  that  a  multiple  step  reaction  mechanism  is  involved. 

Reaction  Associated  Noise 

Electrical  noise  was  observed  on  the  CS2*sapphire  gage  and  was  found 
to  correlate  with  the  particle  velocity  decay,  evolution  of  a  two-wave 
structure,  and  reaction  rate  data.  Zt  was  therefore  assumed  to  be  a  direct 
result  of  the  reaction  taking  place,  coupled  perhaps  with  the  magnetic 
field. 
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This  noise  was  found  to  pervade  the  whole  region,  have  a  frequency 


of  approximately  150  MHz  and  be  a  one-time  phenomenon  that  decays.  Although 
very  little  else  is  known,  this  appears  to  support  the  idea  of  a  chemical 
reaction  producing  some  kind  of  charged  species  in  the  shock. 

Reaction  Mechanism 

The  information  from  this  study  does  not  provide  any  information 
about  the  reaction  mechanism  except  to  indicate  by  the  two  step  rate  law, 
the  induction  time,  and  the  noise  that  a  multiple  step  reaction  mechanism 
may  be  involved. 
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V.  TIME-RESOLVED  SPECTRA  OF  SHOCKED  LIQUID  CS2 
5.1  Introduction 


If  one  had  free  choice  of  spectral  range  to  be  studied  in  attempting 
to  arrive  at  an  understanding  of  shock-induced  chemical  reactions,  he  would 
probably  choose  to  study  vibrational  and  rotational  levels  in  the  infrared. 
Severe  experimental  constraints  imposed  by  the  shock  environment  make  this 
a  poor  second  choice  at  present.  Problems  of  visible  and  near-ultraviolet 
spectroscopy  appeared  to  be  sufficiently  challenging  to  engage  our  full 
attention  at  the  outset  of  this  program.  The  techniques  used  to  obtiin 
spectra  in  this  region  and  results  to  date  are  described  here. 

Carbon  disulphide  is  transparent  in  the  visible  but  strongly  absorbing 
in  portions  of  the  UV  spectrum  (Fig.  1).  Kleman^  describes  absorption  bands 

O 

in  gaseous  carbon  disulfide  between  3300  and  4300  A.  Those  bands  to  the  red 


side  of  about  3660  A  are  too  weak  to  be  observed  in  the  time  resolved  spectra 

for  the  current  experiments  when  the  samples  are  at  ordinary  temperature 

and  pressure.  However,  they  may  be  significant  as  temperature  and  pressure 

are  increased.  This  possibility  will  be  discussed  later. 

Kleman's  analysis  of  the  spectrum  divides  it  into  several  systems 

which  correspond  to  different  electronic  excitations.  He  labels  these  the  R, 

S,  and  V  systems,  and  he  calls  some  unassigned  bands  the  U  bands-  Jungen, 

2 

Malm,  and  Merer  continue  to  use  Kleman's  R,  S,  and  V  system  notation  and 
they  Identify  another  electronic  upper  state  and  label  transitions  to  this 
state  as  the  T  system.  This  system  Includes  Kleman's  U  bands.  The  lower 
state  for  each  of  these  systems  is  the  linear  molecule  ground  state  with 
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various  levels  of  the  vibrational  bending  mode,  v^,  excited.  Transitions 
from  excited  levels  of  the  vibrational  stretching  modes,  and  v^»  were  not 
observed  by  Kleman. 

O 

The  R  system  extends  from  4300  to  3275  A.  The  upper  state  for  these 
transitions  is  a  bent  molecule  with  bending,  v^,  and  symmetric  stretching, 

O 

v^,  vibrations  excited.  The  S  system  extends  from  3700  to  3350  A.  The 
bending  mode  is  the  only  excited  vibrational  upper  state  observed  in  the 
S  system.  The  spectral  extent  of  each  system  is  shown  in  Fig.  1. 

The  I  and  V-sys terns  are  overlapping  in  the  region  between  3400  and 
2900  £  with  the  T-system  predominating  in  the  red  end  of  this  region. 

Jungen  et  al.  identify  the  T  and  V-systeas  as  two  Renner-Teller  components 
of  a  transition  to  a  configuration  which  is  the  bent  molecule  analog  of  a 
*Au  configuration  in  a  linear  molecule.  The  syuaetries  of  these  two  states 
are  *A2  for  the  T-system  and  for  the  V-system  The  absorption  in  the 
T  and  V-systems  is  much  stronger  than  in  the  R  and  S-systems. 

Douglas  and  Zanon  4  describe  and  partially  analyze  a  system  of  bands 
which  they  observe  between  2300  X  and  1850  £.  They  call  them  the  2100  £ 
bands.  As  with  the  R,  S,  T,  and  V-system s,  the  2100  £  bands  are ' attributed 
to  transitions  to  a  bent  excited  state.  This  system  is  much  stronger  than 
any  of  the  R,  S,  T,  or  V-systema. 

Hot  bands  are  features  in  the  absorption  spectrum  of  a  substance 

which  increase  in  strength  when  the  temperature  is  increased.  They  have 

12  4 

considerable  effect  on  the  carbon  disulfide  spectrum.  ’  ’  Kleman 
recorded  various  regions  of  the  spectrum  at  295°K  and  580°K. 
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LIQUID 


show  the  absorbance 


AC  295°K  Che  long  wavelengch  end  of  che  spectrum  was  extremely  weakly 

absorbing  compared  co  che  shore  wavelengch  end.  Ac  Che  long  end,  che  spec- 

Crum  was  recorded  Co  4100  £  using  a  sample  densicy  of  9.1  g  of  sample  per 

square  cencimecer  of  che  lighe  beam.  BuC  ac  che  shore  wavelengch  end,  3100  £, 
-3  2 

only  1.2x10  g/cm  were  required.  The  racio  of  chese  sample  densicies, 

3 

7x10  ,  gives  a  rough  idea  of  Che  large  difference  in  absorption  beeween  Che 
R-system  and  che  T  and  V-sysCems. 

When  Che  cemperaCure  was  increased  Co  580°K,  a  cross  seccional  sample 
-3  2 

densicy  of  only  2.5  x  10  g/cm  was  required  co  obtain  spectra  from  around 
4100  £  out  co  4300  £.  This  would  indicaCe  thac  ac  580°K  Che  hoc  bands  of  che 
R-system  are  becoming  comparable  in  strengch  Co  Che  T  and  V-sysCems  ac 
295°K. 

The  encire  T-sysCem  is  composed  of  hoc  bands  in  Jungen's  analysis  of 
che  T  and  V-sysCems.  This  is  because  dipole  transitions  from  che  linear 

ground  scace  co  a  state  with  che  A-  symmetry  of  che  T  state  is  symmetry 
8  * 

forbidden.  So  ic  is  necessary  Co  excice  a  bending  mode  in  Che  ground  scace 
Co  see  any  dipole  transitions  Co  che  T  scace.  The  T-system  comprises  nearly 
che  encire  red  end  of  che  T-V  absorpeion  region;  so  enhancing  che  strength 
of  che  T-system  by  increasing  Che  temper acure  would  make  the  apparent  edge 
of  che  T-V  syscem  stove  co  che  red. 

Douglas  and  Zanon  have  observed  ehac  many  of  Che  bands  Co  che  red 
side  of  2150  £  in  Cheir  2100  £  system  are  hoc  bands  since  Chey  see  a  decrease 
in  Che  strength  of  Che  absorptions  when  Che  sample  is  cooled.  These  bands 
appear  in  progressions  showing  che  energy  differences  between  Che  vibrational 
bending  modes  of  che  ground  aCaCe. 

In  each  of  Che  above  systems  of  transitions,  thermally  excited 
bending  vibrations  of  che  linear  ground  scace  significantly  enhance 


transitions  to  the  bent  upper  states.  These  hot  bands  contribute 
significantly  to  the  absorption  spectrum  of  carbon  disulfide,  particularly  on 
the  red  ends  of  the  various  absorption  regions. 

Liquid  carbon  disulfide  has  two  major  absorption  features  within  the 
spectral  range,  2500  to  7000  $,  accessible  in  these  experiments.  The  shorter 
wavelength  feature  is  a  very  strong  absorption  with  its  red  edge  between 
about  2590  and  2660  X  depending  on  the  thickness  of  the  sample  cell.  This 
region  will  be  referred  to  as  absorption  region  3.  The  blue  edge  of  this 
region  is  outside  the  spectral  range  available.  The  longer  wavelength 
feature,  region  A,  is  a  moderately  strong  absorption  region  between  about 
2900  &  and  3580  1  when  seen  through  a  0.15  mm  thick  cell.  In  the  thinnest 
cells,  <1  ym,  this  region  passed  enough  light  to  cause  some  exposure  of  the 
film  even  at  the  wavelength  of  maximum  absorption,  which  was  about  3200  X. 

The  regions  of  absorption  in  liquid  carbon  disulfide  correspond 
roughly  with  the  regions  of  maximum  absorption  in  the  gas  phase  spectra  as 
shown  in  Figure  1.  However,  significant  differences  are  noted.  The  most 
striking  of  these  is  that  both  absorption  regions  extend  farther  to  the  red 
than  their  gas  phase  counter  parts  if  one  ignores  the  very  weak  R  and 
S-systems.  Region  B  extends  to  between  2590  £  and  2660  X  while  region  A 
extends  to  about  3580  X,  about  180  X  to  the  red  side  of  the  strong  T-V 
systems  in  the  gas  phase  spectra. 


-  v 


Gas  phase  studies  of  CS2  decomposition  indicate  that  the  CS  molecule 
is  one  of  the  earliest  appearing  decomposition  products.  The  CS  molecule 
has  an  absorption  band  in  the  gap  between  absorbing  regions  A  and  B 
shown  in  Fig.  1.  It  was  our  optimistic  believe  at  the  start  of  this 
program  that  passage  of  a  shock  would  cause  transmission  in  this  gap  to 
fade  away,  thereby  enabling  us  to  determine  the  rate  of  formation  of  CS. 

These  expectations  were  misdirected  because  of  the  appearance  of  other 
effects  which  are  perhaps  even  more  interesting,  as  will  be  related. 

Sheffield's  experiment,  described  in  Section  IV,  was  designed  so  that 
the  sample  would  be  thrust  into  the  reacting  state  by  the  first  reflected 
shock.  Cell  thickness  was  great  enough  to  allow  it  to  remain  in  this  state 
for  the  better  part  of  a  microsecond  and  also  great  enough  that  EMV  gages 
could  be  suspended  in  the  cell.  It  had  two  major  disadvantages  for 
spectroscopy.  First  was  that  the  front  of  the  cell,  being  PMMA,  was  opaque 
in  the  UV.  Second  was  that  cell  thickness  was  very  great  compared  to 
absorption  lengths  in  liquid  CS^.  The  latter  condition  would  make  it 
difficult  to  unravel  any  time-dependent  spectral  changes  occurring  in  the 
region  behind  the  shock  front.  These  considerations  made  it  necessary  to 
redesign  the  cell  in  the  manner  shown  in  Figs.  2-4.  In  Fig.  2,  A  is  the 
impacting  sapphire  which  strikes  the  front-face  sapphire,  B,  at  t  *  0.  D  is 
the  CS^  sample.  The  rear  face  of  the  cell,  E,  is  also  sapphire.  Pressure 
in  D  rings  up  to  its  final  value  through  states  2,  3,  4,  5,  etc.,  shown  in 
Figs.  3  and  4.  In  Fig.  4,  successive  pressure  states  at  the  C-D  boundary  are 
shown  by  the  solid  line;  mean  pressure,  as  seen  by  a  light  ray  normal  to  the 
cell,  is  indicated  by  the  dotted  line.  Duration  of  states  3,  5,  etc.  can 
be  adjusted  by  varying  the  thickness  of  D.  This  permits  considerable 
flexibility  in  varying  pressure  and  temperature. 

.1^ 
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Fig.  3  Pressure  States  in  CS,  resulting  from  the  impact  shown  in 


Fig.  1.2. 
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The  cell  of  Fig.  2  represents  a  compromise  which  is  not  altogether 
satisfactory.  The  final  state  is  reached  through  a  succession  of  shocks, 
which  may  not  produce  the  same  effects  as  a  single  strong  shock.  The 
temperature  reached  when  the  system  rings  up  to  its  final  pressure  is 
closer  to  that  of  adiabatic  compression  then  to  a  single  shock.  If  transit 
time  through  the  CS^  is  large  compared  to  time  resolution  of  the  recording 
instrument,  the  effects  of  successive  shocks  can  be  followed,  thereby 
increasing  the  amount  of  information  obtained  in  a  single  experiment. 

These  have  come  to  be  called  "thick  cell  experiments."  In  thin  cells  the 
transit  time  is  small  compared  to  time  resolution  and  the  consequences  of 
"sudden"  applications  of  pressure  and  temperature  can  be  recorded.  So 
the  experiments  described  have  some  shortcomings,  but  they  open  a  new 
regime  for  study  and  are  useful.  They  should  not  represent  the  last  work 
in  experimental  design. 

In  the  experiments  described  in  the  following  sections  the  basic 
experimental  procedure  consisted  of  impacting  a  cell  containing  the  liquid 
with  a  high  speed  projectile  to  drive  a  shock  wave  into  the  liquid.  Simultan¬ 
eously,  a  high  intensity  pulsed  light  source  would  illuminate  the  cell  and 
the  light  transmitted  through  the  sample  would  then  be  analyzed  by  a 
dispersion  unit  and  a  streak  camera  to  provide  a  time  resolved  spectrum  of 
the  transmitted  light. 

Experiments  were  conducted  using  two  variations  of  this  idea.  Both 
variations  utilized  a  sapphire  sample  cell  with  a  thin  front  or  impact  side 
and  a  thicker  back.  A  sapphire  impactor  would  strike  the  front  of  the  cell 
sending  a  shock  wave  into  the  liquid  in  the  cell  cavity.  In  the  first 
variation,  referred  to  as  reflection  experiments,  light  from  a  xenon  flash- 
lamp  was  focused  into  the  back  of  the  cell.  The  light  was  to  pass  through 
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ch«  cell  and  sapphire  iapactor  and  be  reflected  by  an  aluminum  mirror  vapor 
plated  on  the  back  of  the  iapactor.  The  light  would  then  pass  through  the 
iapactor  and  cell  a  second  time  as  shown  in  Figure  5.  This  arrangement  had 
the  advantage  of  keeping  the  airror  free  of  daaage  by  the  shock  wave  until 
after  the  experiment.  It  was  found  that  the  reflection  from  the  first 
sapphire-liquid  interface  in  the  light  path  was  a  problem  because  it  reached 
the  streak  camera  along  with  the  desired  light  signal  chat  had  passed  through 
Che  cell.  To  get  around  this  problem  a  second  variation  was  developed  in 
which  the  xenon  light  source  was  mounted  behind  the  impactor  inside  the 
projectile  as  shown  in  Figure  6.  These  will  be  referred  to  as  transmission 
or  flashlamp-in-projectile  experiments . 
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Figure  5  —Optical  Path  f or  Reflection  Experiments .  This  cross 
sectional  sketch  of  the  projectile  and  target  shows  the  intended  path  of  Che 
light  beam  during  a  reflection  type  experiment.  A  lens*  noc  shown,  focuses 
an  image  of  the  flash lamp  onto  the  impactor 's  mirror. 
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Figure  6  —Optical  Path  for  Transmission  Experiments .  In  trans- 
mission  experiments,  the  flashlamp  is  mounted  in  the  projectile  so  that 
light  passes  only  one  direction  through  the  sample. 
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5.2  APPARATUS  AND  PROCEDURES 


These  experiments  were  also  carried  out  in  the  Shock  Dynamics  Laboratory 
at  Washington  State  University  using  the  4  inch  gas  gun  facility  augmented  by 
a  Beckman  and  Whitley  model  1500  spectrograph,  a  Beckman  and  Whitley  model 
339B  rotating  mirror  streak  camera,  suitable  light  sources,  and  other  optical 
components.  These  are  described  below. 

A.  THE  CAMERA  ASSEMBLY 

The  Beckman  and  Whitley  model  339B  streak  camera  and  model  1500 
dispersion  unit  bolt  together  into  a  single  camera  assembly.  This  camera 
assembly  is  mounted  above  the  target  tank  and  looks  vertically  down  into 
the  target  tank.  An  expendable  mirror  placed  about  20  cm  behind  the  target 
directs  light  from  the  target  up  into  the  camera  assembly  during  both 
reflection  and  transmission  experiments.  The  camera  assembly,  optics, 
and  other  components  of  the  spectrographic  system  are  described  in  the 
following  sections,  which  refer  to  Fig.  1. 

Dispersion  Unit  Optics 

The  dispersion  unit  optical  path  starts  with  an  11.4  cm  diameter, 

46  cm  focal  length,  objective  mirror  which,  by  way  of  a  relay  mirror  and 
quartz  vacuus  window,  focuses  an  image  of  the  target  onto  a  slit  which  is 
adjustable  in  two  dimensions.  The  slit  dimensions  provide  control  of  both 
spectral  and  temporal  resolution.  Light  passing  through  the  slit  is 
reflected  by  a  relay  mirror  to  a  10.2  cm  diameter,  41  cm  focal  length 
collimating  mirror.  The  collimated  beam  llltninates  the  plane  102  x  127  mm 
diffraction  grating  which  disperses  the  spectral  components.  The  dispersed 


beam  is  focused  by  a  10.2  cm  diameter,  81  cm  focal  length,  collecting 
mirror  to  form  a  dispersed  image  of  the  slit  at  the  output  plane  of  the 
dispersion  unit,  which  is  12  cm  outside  the  dispersion  unit  housing.  This 
output  serves  as  the  input  to  the  streak  camera. 


Streak  Camera  Optics 

The  first  optical  element  in  the  streak  camera  is  a  5.1  cm  diameter, 

33  cm  focal  length,  quartz  field  lens.  This  lens  causes  the  cones  of  light 
for  different  wavelengths  to  converge  in  such  a  manner  as  to  minimize  the 
size  of  the  convex  mirror  needed  to  avoid  vignetting.  This  is  accomplished 
by  constructing  the  field  lens  to  focus  an  image  of  the  collecting  mirror 
onto  the  convex  mirror. 

The  next  element  in  the  streak  camera,  the  convex  mirror,  spreads  the 
dispersed  beam  to  illuminate  a  30.5  cm  diameter  concave  primary  mirror. 

This  mirror  focuses  the  light  to  a  spectrally  dispersed  image  of  the  slit 
112  cm  from  the  primary  mirror.  The 'rotating  mirror  is  located  28  cm  ahead 
of  the  focused  image.  It  causes  the  image  to  sweep  along  a  240  degree  arc 
of  film  as  the  mirror  is  rotated.  The  rotating  mirror  does  not  intercept 
all  of  the  light  in  the  beam.  Some  passes  above  the  rotating  mirror  to  an 
aperture  in  the  back  of  the  camera  casing.  An  RCA  type  4831  photomultiplier 
tube  is  Installed  in  this  aperture  to  monitor  the  total  intensity  of  the 
light  within  the  spectral  range  actually  reaching  the  film.  The  rotating 
mirror  has  three  reflecting  faces  forming  an  equilateral  triangle  so  that 
one  face  is  always  in  position  to  reflect  light  to  the  film.  The  mirror  is 
long  enough  to  intercept  nearly  the  entire  light  cone  in  the  horizontal 
direction.  Each  face  is  7.6  cm  long  and  3.3  cm  high. 


Camera  Assembly  Slit  and  Resolution  Limits 


The  slic  which  controls  both  the  height  and  width  of  the  Image  on 
the  film  Is  located  at  the  Input  to  the  dispersion  unit.  The  streak 
camera  does  not  have  a  separate  slit  of  its  own.  The  width  of  the  slit 
has  a  micrometer  adjustment  indicating  slit  width  in  units  of  10  um. 

This  width  is  parallel  to  the  direction  of  wavelength  dispersion  and  so, 
together  with  the  line  density  on  the  grating,  can  be  used  to  adjust  the 
spectral  resolution  o£  the  apparatus.  The  maximum  slit  width  possible  is 
3  mm,  but  the  design  of  the  optical  system  outside  the  camera  assembly 
assumed  that  neither  the  height  nor  the  width  of  the  slit  would  exceed  1  mm. 

The  height  of  the  slit  aperture  is  adjustable  by  a  stepped  mask  with 
heights  of  0.274  mm,  0.482  mm,  1.978  mm  and  larger  up  to  about  12  mm.  Only 
the  smallest,  0.274  mm,  opening  has  been  used.  The  temporal  dispersion  is 
obtained  by  sweeping  the  image  of  the  slit  in  a  direction  parallel  to  the 
height  of  the  slit  image;  so  the  time  resolution  of  the  instrument  can  be 
adjusted  by  the  slit  height  and  the  rotating  mirror  speed.  The  magnification 
between  the  slit  and  the  film  is  0.99  and  the  spectral  dispersion  spreads 
the  image  out  across  the  25  mm  usable  width  of  the  film. 

The  time  resolution  which  can  actually  be  achieved  is  affected  by  several 
factors  which  include  the  maximum  operating  speed  of  the  rotating  mirror,  the 
spatial  resolution  of  the  film,  and  brightness  of  the  light  source.  The 
Kodak  2475  recording  film  used  is  rather  coarse  grained,  having  a  spatial 
resolution  of  about  10  line  pairs  per  millimeter.  This  together  with  the 
2500  revolutions  per  second  maximum  mirror  speed  and  a  28  cm  sweep  arm 
implies  a  best  time  resolution  of  11  nanoseconds.  At  a  typical  mirror  speed 
of  2480  revolutions  per  second,  the  0.274  mm  slit  yields  a  time  resolution 
of  about  31  nanoseconds. 
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With  a  limited  light  intensity  available  there  is  a  trade  off  to  be 
made  between  spectral  and  temporal  resolution.  To  obtain  enough  light  on 
the  film  when  using  the  600  lines/mm  grating  blazed  for  3000  %t  it  was 
necessary  to  open  the  spectral  slit  to  0.75  mm.  Since  this  grating  gives 
a  dispersion  of  41  S/mm  on  the  film,  and  the  magnification  between  the  slit 
and  the  film  is  nearly  unity,  the  spectral  resolution  is  about  30  X. 
Similarly,  for  the  300  lines /mm  and  150  lines/mm  gratings  the  slit  widths 
were  0.50  and  0.35  mm,  respectively,  giving  the  spectral  resolutions  of 
40  and  55  %.  If  more  light  were  available  or  less  time  resolution  were 
acceptable,  greater  spectral  resolution  could  be  achieved  by  closing  the 
slit  width  down  as  narrow  as  0.1  mm.  For  any  smaller  slit  width  the  spatial 
resolution  of  the  film  would  limit  the  spectral  resolution. 


B.  PHOTOGRAPHIC  FILM  AND  DEVELOPMENT 


The  photographic  film  selected  for  use  in  the  streak  camera  was  Kodak 

recording  film  2475  (ESTAR-AH  Base).^’^  This  film  was  selected  for  its 
very  high  speed  and  its  spectral  response.  This  is  a  coarse  grained  film 

to  obtain  the  high  film  speed.  Its  spatial  resolution  is  only  about  10 

line  pairs  per  millimeter  but  this  is  adequate  considering  the  slit  sizes 

which  are  used. 

Several  other  films  were  tried.  One  of  these,  Kodak  pan  film  2484,  was 
nearly  equivalent  to  the  2475  film.  The  film  developing  process  was  adjusted 
to  obtain  the  maximum  contrast  and  sensitivity  without  excessive  increase  in 
the  film’s  fog  level.  The  process  that  was  used  for  all  film  records  is: 

1.  Develop  Kodak  2475  film  in  Kodak  DK-50  developer  at 
20°C  for  7  min.  Agitate  for  15  sec  after  adding  the 
developer  and  every  2  min  thereafter. 

2.  Stop  development  using  Kodak  indicator  stop  bath 
at  20°C.  Pour  in  the  stop  bath,  agitate  for 

30  sec,  and  drain. 

3.  Fix  using  Kodak  fixer  at  20°C  for  10  min.  Agitate 
for  15  sec  every  2  min. 

4.  Wash  in  20°C  running  tap  water  for  25  min. 

5.  Rinse  with  distilled  water  and  hang  to  dry. 
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C.  THE  OPTICAL  SETUP  AND  ITS  DESIGN 

The  most  important  consideration  in  the  design  of  the  optical  system 
between  the  light  source  and  the  dispersion  unit  is  to  maximize  the  total 
amount  of  light  which  can  be  transmitted  through  the  target  and  into  the 
usable  light  cone  of  the  camera,  subject  to  the  constraints  imposed  by  the 
physical  setup  of  the  experiment.  One  of  these  constraints  is  that  only  a 
6  mm  diameter  area  in  the  center  of  the  target  cell  remains  free  from  side 
rarefactions  for  the  duration  of  the  experiment,  about  1  us.  Another  is 
that  the  closest  the  dispersion  unit  can  be  mounted  to  the  target  places 
the  objective  mirror  at  an  optical  path  length  of  2.9  m  from  the  target. 

This  2.9  m  optical  path,  together  with  the  46  cm  focal  length  of  the 
objective  mirror  implies  a  magnification  between  the  target  and  slit  of 
1/5.3  so  that  a  5.3  mm  dimension  on  the  target  is  imaged  into  1  mm  on  the 
slit.  Since  1  mm  is  as  large  a  slit  opening  as  is  needed  and  5.3  mm  is 
still  smaller  than  the  rarefaction-free  center  of  the  target  cell,  no 
additional  optics  are  needed  between  the  target  and  the  object  mirror. 

For  reflection  experiments  a  lens  must  be  selected  to  focus  the  light 
source  onto  the  target.  The  lamp  in  this  case  is  mounted  outside  the  target, 
so  a  convenient  location  for  the  lens  is  in  a  target  tank  access  port 
10  cm  behind  and  69  cm  to  the  side  of  the  target.  The  diameter  of  this 
lens  must  be  great  enough  to  provide  a  cone  of  light  that  will  completely 
illuminate  the  objective  mirror.  At  79  cm  from  the  image  on  the  target 
this  requires  a  diameter  of  2.10  cm  plus  an  additional  0.66  cm  for  off- 
axis  rays  focused  at  the  edge  of  a  1  mm  wide  slit.  Some  additional  lens 
diameter  is  desirable  to  make  alignment  less  critical  and  to  allow  for  some 
tilt  in  the  target  mirror.  A  6.35  cm  diameter  lens  was  selected  with  0.63  cm 
around  the  perimeter  covered  by  the  lens  mount,  which  must  be  rather  strong 


because  che  lens  serves  as  a  vacuum  and  pressure  boundary  for  Che  target 
tank.  The  final  effective  diameter  is  5.08  cm. 

The  lens  to  flashlamp  distance  is  selected  so  that  the  image  of  the 
flashlamp  will  completely  cover  the  rarefaction-free  center  of  the  target, 
and  thus  completely  cover  the  slit.  Using  a  15.2  cm  focal  length  lens, 
the  lens  to  flashlamp  distance  is  18.8  cm.  Illuminating  the  6  mm  diameter 
rarefaction-free  region  thus  requires  a  flashlamp  with  1.4  mm  bore  diameter. 
With  a  2  or  3  mm  bore  flashlamp  this  setup  allowed  reasonable  leeway  in 
alignment  so  a  commercially  available  15.2  cm  focal  length  plano-convex 
ultraviolet  lens  and  a  commercially  available  3  mm  bore  flashlamp  were 
selected. 

When  it  was  decided  to  mount  the  flashlamp  in  the  projectile,  a  3  mm 
bore  lamp  was  selected.  Although  this  is  not  as  large  as  the  6  mm 
rarefaction-free  region,  it  is  large  enough  for  its  0.57  mm  image  to 
completely  cover  the  0.274  mm  height  of  the  slit.  The  0.57  mm  image  will 
not  cover  the  maximum  slit  width  used;  the  flashlamps  in  the  projectile 
are  oriented  so  that  the  long  axis  of  the  image  of  the  lamp  falls  along 
the  longest  dimension  of  the  slit.  With  this  arrangement  no  optical  elements 
other  than  a  flat  relay  mirror,  are  needed  between  the  projectile  light 
source  and  the  objective  mirror  of  the  dispersion  unit. 


D.  THE  SAMPLE  CELL 


Design  Considerations  for  Reflection  Experiments 

The  first  considerations  were  that  the  cell  be  transparent  to 
ultraviolet  light  at  pressures  above  100  kbars  and  that  its  shock  properties 
be  known  in  this  pressure  range.  These  criteria  led  to  the  selection  of 
Z-cut  sapphire  as  the  principal  cell  material.  Barker  and  Hollenbach^ 
describe  Z-cut  sapphire  properties  and  note  that  it  is  useful  as  a  window 
material  up  to  pressures  as  high  as  130  to  150  kbars. 

Another  consideration  was  the  amount  of  light  that  would  be  transmitted 
by  the  sample  material  in  the  cell.  The  brightness  available  from  the 
flashlamp  was  only  marginally  adequate  for  the  time  resolution  needed  so 
it  was  necessary  to  use  a  very  short  optical  path  within  the  cell.  This 
limitation  was  compounded  by  the  fact  that  the  light  would  make  two  passes 
through  the  cell  and  by  the  high  absorption  of  carbon  disulfide  within 
some  wavelength  regions  of  interest.  Consideration  of  these  limitations  led 
to  the  selection  of  a  nominal  cell  thickness  of  0.1  mm. 

It  was  necessary  to  prevent  light  reflected  from  the  back  surface  of 
the  cell  from  reaching  the  camera  assembly.  About  0.08  of  the  incident  light 
is  reflected  from  the  first  vacuum-sapphire  interface  at  the  back  of  the 
cell.  This  reflected  light  has  not  reached  the  carbon  disulfide  layer  so 
carries  ro  useful  information..  If  the  back  surface  were  parallel  to  the 
mirror  behind  the  impactor,  the  back  surface  reflection  would  be  parallel  to 
and  essentially  coincident  with  the  reflection  from  the  cell  mirror.  Both 
reflected  beams  would  enter  the  camera  assembly.  To  prevent  this,  the  back 
piece  of  sapphire  was  made  wedge  shaped  with  its  back  side  tilted  six  degrees 
from  parallel  to  its  front.  The  Z  crystal  axis  is  perpendicular  to  the  front 


face. 


The  question  of  where  to  place  the  mirror  in  the  cell  assembly 

received  considerable  attention.  The  mirror  could  have  been  vapor  plated 
onto  the  sapphire  at  the  front  of  the  cell  cavity.  However,  it  was  feared 
that  the  specular  reflective  quality  of  the  mirror  would  be  lost  when  the 
shock  reached  mirror.  In  fact  even  a  small  effect  on  the  mirror  could 
not  be  tolerated  because  it  would  be  impossible  to  tell  whether  a  change 
observed  in  the  light  beam  was  due  to  the  response  of  the  sample  to 
the  shock  wave  or  due  to  the  response  of  the  mirror.  To  avoid  these  problems 
it  was  decided  to  place  the  mirror  behind  the  sapphire  impactor  in  Che 
projectile  where  it  would  remain  unshocked  for  the  significant  part  of  the 
experiment.  There  was  still  some  question  as  to  the  transmission  quality  of 
the  sapphire-sapphire  impact  interface.  Subsequent  results  showed  a  slight 
loss  of  light  at  impact  but  not  enough  to  cause  concern. 

The  elastic  shock  velocity  in  sapphire  is  about  11  mm/usec.  To  obtain 
at  least  a  6  mm  rarefaction  free  volume  in  the  center  of  the  cell  for  the 
duration  of  a  1  usee  experiment,  a  cell  diameter  of  31.8  mm  was  selected. 

Also,  with  this  elastic  shock  velocity  in  mind  a  12.7  mm  chick  cell  back  was 
selected  to  prevent  the  shock  from  reaching  the  back  of  the  cell  within  1  usee 
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C ell  Construction 

Having  decided  on  Che  desired  features  of  Che  sample  cell  1C  was 
necessary  Co  develop  Che  final  design  and  Che  conscrucdon  cechnlques. 

•  ?  shows  Che  cell  design  used. 

The  fronc  of  Che  cell  served  only  Co  confine  Che  concencs,  so  ic 

was  made  as  chin,  2.07  mm,  as  possible  wlchouc  making  1c  Impossible  Co 

mounC  in  a  supporting  sCrucCure.  Mouncing  Che  cell  fronc  was  complicaced 
by  Che  need  Co  allow  Che  projecCile  Co  lmpacc  ic  wlchouc  impedimenc  on  one 
side  and  Che  need  Co  mounC  ic  only  0.1  mm  from  Che  back  sapphire.  This  lefc 

only  che  edge  of  che  cell  fronc  Co  use  for  holding  ic  in  place.  The  firsc 

cells  used  epoxy  co  hold  Che  fronc  sapphire  in  a  brass  ring.  The  brass  was 
made  slighcly  chinner  chan  che  fronc  sapphire  and  one  side  was  lapped  flac. 

The  brass  ring  and  sapphire  were  clamped  againsc  a  lapped  flac  steel 
place  as  che  epoxy  cured.  This  ensured  chac  che  sapphire  and  brass  surfaces 
were  in  che  same  plane. 

The  sapphire  cell  back  was  also  mounced  wich  epoxy  in  a  brass  ring. 

Afcer  che  epoxy  cured,  che  assembly  was  chucked  in  a  lathe  and  the  brass 
was  faced  off  to  form  a  surface  parallel  Co  che  sapphire  surface  but  displaced 
outward  by  an  amount  equal  to  Che  desired  cell  Chickness,  usually  about  0.1  am. 
Then,  the  flat  side  of  the  cell  front  assembly  was  placed  againsc  this  surface. 
A  cavity  about  0.1  mm  thick  remained  between  che  front  and  the  back  sapphires. 
To  allow  for  filling  Che  completed  cell,  two  0.25  mm  holes  were  drilled  in 
the  brass  at  opposite  edges  of  che  back  sapphire.  The  outer  ends  of  the 
holes  were  enlarged  to  accept  1.6  mm  copper  cubing;  and  the  inner  ends  were 
connected  to  the  cell  cavity  by  milled  slots  in  the  brass.  To  seal  the 
assembled  cell,  an  "0"  ring  groove  was  turned  into  che  face  of  che  back 
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CELL  BACK  Scale  1.5*  I 

Figure  7  —Cell  Design.  The  sapphire  cell  windows 
were  mounted  in  these  brass  rings. 


brass.  Eight  machine  screws  just  outside  the  "0"  ring  held  the  cell 
together.  The  6-32  flathead  screws  were  threaded  into  the  back  brass 

with  their  heads  countersunk  into  the  front  of  the  cell. 

The  narrowness  and  brittleness  of  the  epoxy  joint  between  the  front 

sapphire  and  its  brass  ring  caused  some  problems.  As  the  front  sapphire 
was  tightened  over  the  "0"  ring,  the  brass  would  flex  slightly  and  the 
epoxy  would  crack,  causing  leaks.  Also,  it  was  feared  that  sample 
materials  such  as  CS^  would  Leach  contaminants  out  of  the  epoxy.  So, 
on  all-metal  thermal-shrink  method  of  mounting  the  sapphire  in  the  brass 
rings  was  developed. 

The  hole  in  the  brass  ring  was  made  0.05  mm  smaller  than  the 
sapphire  and  the  inner  brass  surface  was  coated  with  a  thick  layer  of  tin- 
lead  solder.  Flux  residue  was  removed;  then  this  ring  was  heated  on  a  flat 
lapped  steel  plate  until  the  hole  had  expanded  enough  to  allow  the  sapphire 
to  be  dropped  in.  The  sapphire  and  brass  were  then  firmly  held  against  the 
flat  plate  as  the  plate  was  cooled  by  setting  it  on  a  thick  cold  metal 
surface.  This  technique  developed  a  strong  bond  between  the  sapphire  and 
brass.  At  first,  two  or  three  tries  were  needed  to  produce  a  Leak  free 
assembly  that  was  acceptable  in  all  aspects.  Production  efficiency  was 

improved  by  placing  a  tight  specification  of  13  urn  total  indicated 
runout  on  the  roundness  of  the  sapphire,  and  by  improved  skill.  Both  front 
and  back  sapphires  were  mounted  with  this  heat-shrink  technique.  The 
copper  fill  tubes  were  soldered  into  these  cells. 

The  "0"  ring  material  must  be  compatible  with  the  liquid  in  the  cell. 

_  7  7 

Teflon  was  used  in  early  cells  with  carbon  disulfide.  Viton  was  used  in 

later  cells  because,  being  a  more  rubbery  material,  it  seals  better  with 

less  force,  thereby  making  accurate  control  of  cell  thickness  easier.  Teflon 

and  Viton  are  resistant  to  carbon  disulfide.  Some  other  "0"  ring  materials 


are  not. 


Thin  Cell  Modifications 

Several  cells  in  the  range  of  0.6  to  1.2  ym  thickness  have  been  constructed. 

These  extremely  thin  cells  were  achieved  by  modifying  the  manner  in 

which  the  cell  front  was  attached  to  the  back.  Instead  of  facing  off  the 

o 

back  brass  flat  with  the  sapphire,  the  brass  was  faced  off  at  a  1/2 
angle  tilted  up  away  from  the  sapphire.  Also,  the  inner  edge  of  the 
brass,  next  to  the  sapphire,  was  recessed  25  um  below  the  sapphire  surface 
as  shown  in  Fig.  S. 

When  the  cell  front  assembly  was  screwed  onto  the  tilted  brass,  it 
took  a  dish  shape.  The  heat-shrink  assembly  of  the  front  piece  is  strong 
enough  to  force  some  curvature  in  the  front  sapphire .  The  eight 
screws  were  carefully  adjusted  so  that  the  front  sapphire  was  symmetrically 
curved  toward  the  back  sapphire  and  touched  the  back  sapphire  in  the  center. 

It  was  easy  to  tell  when  the  cell  had  been  adjusted  to  this  condition  by 
doing  the  work  under  a  monochromatic,  long-coherence  length  light.  A 
beautiful  set  of  Newton's  rings  became  apparent.  When  che  cell  was  prepared 
for  the  experiment  in  the  target  tank,  it  had  a  vacuum  around  it  and  a  known 
pressure  inside.  The  known  pressure  was  achieved  by  having  the  cell 
connected  to  a  reservoir  of  the  sample  liquid  outside  the  target 
tank  at  atmospheric  pressure.  To  complete  adjustment  of  the  cell 
thickness,  at  the  bench,  this  operating  pressure  in  the  cell  was 
simulated  by  using  compressed  nitrogen.  The  screws  were  adjusted 
until  the  fringes  in  the  cell  showed  it  to  be  the  desired  thickness 
when  there  was  one  atmosphere  of  differential  pressure  between  the 
Inside  and  the  outside  of  the  cell.  Thickness  was  checked  by  pressing 
on  the  front  of  the  cell  until  the  front  and  back  sapphire  made 
contact,  then  counting  the  fringes  as  the  force  was  removed. 
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Figure  8  .-“Machining  of  Thin  Cell  Brass.  Dimensions  have  been 
exaggerated  here  to  show  how  the  back  brass  of  1  ym  thick  cells  was  machined. 
The  front  sapphire  and  brass  took  a  very  slight  dish  shape  when  screwed  onto 
the  tilted  brass  back.  Adjusting  the  torque  on  the  screws  allowed  submicron 
control  over  the  spacing  between  the  sapphire  windows. 
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When  che  cell  could  not  be  adjusted  with  reasonable  torque  on  the  screws. 
Chin  pieces  of  wire  were  used  as  shins  between  the  front  and  back  brass, 
under  che  outside  of  some  or  all  of  the  screws.  Considerable  effort  was 
needed  to  achieve  final  adjustment  of  these  thin  cells,  but  the  results 
were  stable,  reproducing  themselves  well  when  the  internal  pressure  in  the 
cells  was  removed  and  reapplied.  Cells  thinner  than  1  um  in  the  center 
region  probed  by  che  light  beam  were  produced  using  this  technique. 


E.  TARGET  ASSEMBLY 


When  Che  sample  cell  was  completed  it  was  epoxied  into  a  steel  and 
epoxy  target  assembly  as  shown  in  Fig.  9.  This  was  used  to  mount  the 
cell  at  the  muzzle  of  the  gun.  The  steel  ring  part  of  the  assembly  was 
12.7  cm  i.d.,  15.2  cm  o.d.,  and  about  2.5  cm  from  front  to  back.  The 
outer  6  mm  of  the  ring  was  a  lip  about  6  mm  Chick.  Brass  strips  were 
clamped  against  this  lip  to  hold  the  assembly  firmly  against  the  target 
holder  at  the  muzzle.  After  impact  the  brass  strips  would  bend  back 
allowing  the  entire  assembly  to  break  away  without  damage  to  the  target 
holder.  The  interior  of  the  ring  was  filled  to  a  depth  of  1.6  to  1.9  cm 
with  an  epoxy  potting  compound  containing  powdered  calcium  carbonate  filler. 
The  front  or  impact  face  of  this  assembly  was  lapped  flat  before  installing 
the  cell. 

In  addition  to  Che  sample  cell,  several  other  components  were  mounted 
in  the  target  assembly.  These  included  the  trigger  pin,  the  high  voltage 
contacts,  and  the  rubber  insulating  ring.  The  high  voltage  contacts, 
the  insulating  ring,  and  the  groove  in  which  the  insulating  ring  was 
mounted  were  omitted  from  targets  used  in  the  reflection  experiments. 

Brackets  for  mounting  the  mirrors  behind  the  target  were  also  attached 
to  the  target  assembly  after  it  was  mounted  on  the  target  holder. 

After  the  face  of  the  target  assembly  was  lapped,  a  hole  5.79  cm 
in  diameter  was  bored  through  it  for  mounting  the  cell.  In  transmission 
experiment  targets,  a  groove  6.8  mm  deep,  6.03  cm  i.d.,  and  10.48  cm  o.d. 
was  turned  into  the  front  face  surrounding  the  hole  for  the  cell.  Next, 
three  holes  were  drilled  in  the  epoxy  for  the  high  voltage  pins  and  the 
trigger  pin.  The  two  high  voltage  pin  holes  were  tapped  with  8-32  threads. 
The  trigger  pin  hole  was  1.67  cm  in  diameter  and  the  brass  insert  for 
mounting  the  trigger  pin  was  epoxied  into  it.  The  three  holes  were  120° 
apart  and  centered  on  a  circle  of  8.10  cm  diameter. 


The  black  rubber  ring,  6.35  cm  i.d.,  9.86  cm  o.d.,  and  5.3  mm 
Chick,  was  glued  into  Che  groove  in  Che  face  of  Che  cargec  assembly  using 
silicone  rubber  sealanc  and  Eascman  adhesive.  A  small  amounc  of  silicone 
sealanC  was  used  under  che  ring  around  Che  high  volcage  pin  holes  co  ensure 
good  insulacion.  The  Eascman  910  was  used  around  che  resc  of  che  ring. 

Six  millimecer  diamecer  holes  were  cue  in  che  rubber  ring  for  che  high  volcage 

pins  and  a  one  cencimecer  hole  for  che  crigger  pin. 

To  mounc  che  cell  in  che  cargec  assembly,  che  assembly  was  firsc 
clamped  face  down  on  a  lapped  flac  sceel  place.  The  edges  of  che  cell 
and  che  hole  in  which  ie  was  mounced  were  coaced  wich  a  chin  layer  of  Epon 
815  epoxy.  Then  che  cell  was  inserced  inco  che  hole  uncil  che  fronc  sapphire 
was  rescing  flac  againsc  che  sceel  place.  The  cell  was  clamped  in  chis 
posicion  uncil  Che  epoxy  cured. 

The  high  volcage  pins  were  8-32  chreaded  sceel  rods  about  3.8  cm 
long.  A  narrow  capered  poinc,  30°  cocal  included  angle,  was  machined  onCo 
each  pin.  These  pins  were  chreaded  chrough  cheir  holes  uncil  che  poinced 
end  projected  0.25  mm  beyond  che  rubber  ring.  This  scill  left  them  about 
1.17  mm  below  che  surface  of  Che  cargec.  Quick  setting  epoxy  and  a  nut 
were  applied  to  che  pins  on  the  back  side  of  the  target  to  hold  them  in 
position. 

The  high  volcage  wires  were  connected  Co  Che  pins  wich  two  nuts.  A 
plastic  Cube  was  attached  around  the  wire  connecCion  and  filled  wich 
epoxy  so  chat  no  uninsulated  conductor*  were  exposed  on  che  back  of  che 
cargec.  The  high  volcage  wires,  which  were  18  gauge  high  volcage  test  probe 
wire,  were  led  out  of  the  cargec  tank  by  epoxy ing  them  into  holes  in  a  PMMA 
feed-chrough  in  a  flange  on  an  access  port. 
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F.  THE  LIGHT  SOURCE 


Reflection  Experiment  Flashlamps 

The  light  sources  chosen  for  the  reflection  experiments  were  xenon 
flashlamps  of  one  inch  arc  length,  3  mm  bore,  450  torr  xenon  fill,  and 
1  mm  thick  Suprasil  envelope  obtained  from  ILC  Technology,  Inc. 
under  model  number  L-3958.  The  lamps  were  driven  by  a  3.3  uf  capacitor 
bank  charged  to  6000  V.  The  total  stored  energy  was  60  joules.  A 
triggered  spark  gap  controlled  the  timing  of  the  discharge.  The  timing 
of  the  spark  gap  trigger  signal  was  obtained  from  a  charged  shorting  pin 
in  the  target  assembly.  To  help  ensure  rapid  and  consistent  triggering 
of  the  flashlamps,  the  high  voltage  trigger  signal  to  the  spark  gap  was 
also  connected  to  a  trigger  wire  on  the  flashlamp. 


Flashlamps  for  the  Projectiles 

Commercial  flashlamps  were  found  to  be  highly  unreliable  when  mounted 
in  the  projectile.  The  reasons  why  the  commercial  flashlamps  did  not  work 
are  somewhat  speculative  but  two  possibilities  seem  likely.  First,  the 
flashlamp  is  fairly  long  so  the  envelope  may  fail  when  subjected  to  stresses 
from  shrinking  of  the  epoxy  it  was  potted  in,  and  from  stresses  induced 
by  the  acceleration  in  the  gun.  One  lamp  failed  from  epoxy  shrinking  stresses 
alone.  This  failed  lamp  was  potted  in  a  single  thick  pour  of  epoxy.  Subse¬ 
quent  lamps  were  potted  using  three  smaller  pours.  The  second  possibility 
comes  from  the  fact  that  the  electrodes  extend  17  to  20  mm  without  support 
into  the  cavity  of  the  flashlamp,  perpendicular  to  the  direction  of  accelerat 
in  the  gun.  It  is  possible  that  these  o’.ectrodes  broke  during  acceleration 
possibly  breaking  the  envelope  too. 
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A  special  flashlamp  design  was  developed  for  use  in  Che  projectiles. 
Keeping  the  possible  failure  modes  in  mind,  the  new  lamps  were  designed 
as  short  as  possible  and  with  sturdy  electrodes  projecting  only  about  2  mm 
into  the  lamp  cavity  as  shown  in  Figure  10.  Since  the  lamps  were  for  one 

g 

shot  service,  the  quartz  envelope  to  electrode  seals  were  made  with  Torr  Seal 
high  vacuum  epoxy  rather  than  quartz  to  metal  seals.  Two  3.2  mm,  2Z  thoriated 
tungsten  welding  rods  were  used  for  electrodes  and  a  1.6  ms  copper  tube  was 


epoxied  in  next  to  one  of  the  electrodes  for  evacuating  and  back  filling 
with  xenon  to  450  torr.  This  design  proved  perfectly  reliable. 


1/16  COPPER 
FILL  TUBE 


FUSED  QUARTZ 
TUBEv 

/"  = 


ELECTRODES  (2) 


TORR  SEAl 


TRIGGER 

WIRE 


Scale  1.5: 1 

Figure  10  .—Flash lamp  for  Mounting  in  Projectile. 

Torr  Seal  is  a  high  vacuum  epoxy  with  high  viscosity  which  made 
it  convenient  for  sealing  the  electrodes.  The  electrodes  were 
1.4  inch  sections  of  1/8  inch  thoriated  tungsten  welding  electrodes 


G.  THE  PROJECTILE 


Projectiles  for  Reflection  Experiments 

Projectiles  for  these  experiments  were  made  by  modifying  the  standard 
aluminum  projectile,^*  which  is  a  one  piece  hollow  aluminum  cylinder  closed  at 
one  end.  It  is  10.16  cm  in  diameter  and  20.32  cm  long.  A  1.9  cm  thick 
aluminum  plate  was  mounted  on  the  open  end  and  the  impactor  was  mounted  in 
the  plate.  Prior  to  mounting  the  impactor,  the  face  of  the  projectile  was 
lapped  flat  and  perpendicular  to  the  sides  of  the  projectile  to  better 
than  0.01  mill ir ad Ians .  To  ensure  that  the  impactor  face  would  be  parallel 
to  the  face  of  the  projectile,  an  "0"  ring  was  placed  in  the  bottom  of  the 
hole  in  which  the  impactor  was  to  be  mounted.  This  held  the  surface  of  the 
impactor  slightly  above  the  surrounding  projectile  face  until  a  lapped 
flat  steel  plate  was  placed  across  the  face  of  the  projectile  and  weighted 
to  compress  the  "0"  ring.  This  held  the  impactor  flat  against  the  steel  plate 
until  the  surrounding  epoxy  cured.  In  this  way  the  Impactor  was  mounted  with 
its  surface  flush  with,  and  parallel  within  0.1  milliradlan  to  the  face  of 
the  projectile. 


To  provide  insulation,  the  f las hi amp  and  all  ocher  projectile  flash 
head  parts  were  potted  into  a  4 . 1  cm  thick  section  of  polymethylmethacrylate 


(PIMA) ,  which  was  bolted  and  epoxled  onto  the  front  end  of  a  standard 
aluminum  projectile.  Figures  11  and  12  show  the  flash  head  construction. 
Contact  with  the  high  voltage  pins  in  the  target  was  made  by  two  brass  contact 
strips  mounted  in  a  PMMA  contactor  ring,  so  that  when  the  contactor  ring  hit 
the  rubber  ring,  the  brass  strips  and  the  high  voltage  were  electrically 
insulated  from  the  surrounding  vacuum,  thereby  preventing  any  unwanted 
discharges . 

In  addition  to  the  two  high  voltage  contact  strips  in  the  contactor 

ring,  there  was  a  third  contact  strip  which  extended  to  the  edge  of  the 

projectile  diameter.  This  strip  served  to  make  contact  with  the  velocity 

pins  and  the  trigger  pin  in  the  target.  Each  of  the  three  contact  strips 

covered  a  projectile  rotational  angle  of  90°  so  that  projectile  rotation 
o 

up  to  45  could  be  tolerated.  Projectile  rotation  was  usually  less  than 
10°  and  always  less  then  20d.  Six  flathead  6-32  brass  screws  passing 
through  the  contact  strips  and  the  contactor  ring  were  used,  along  with 

quick  setting  epoxy,  to  mount  the  contactor  ring  to  the  projectile. 

To  make  electrical  contact  with  the  lamp,  one  of  the  bolts,  passing 
through  each  high  voltage  contact  strip,  was  threaded  into  a  brass  plug 
potted  into  the  PtMA  projectile  flash  head.  A  12  gauge  copper  wire  was 
connected  to  the  other  end  of  each  brass  plug  by  a  screw.  The  other  end 
of  each  wire  was  soldered  to  a  brass  fitting  clamped  to  a  flashlamp  electrode. 
The  final  part  of  the  lamp  assembly  was  a  thin  nickel  trigger  wire  which  was 
wrapped  around  the  flashlamp  envelope  just  outside  the  tip  of  the  anode,  run 
dose  dong  the  envelope  on  che  side  sway  from  the  Impactor,  wrapped  around 
the  envelope  dose  over  the  cathjd^^then  connected  to  the  cathode . 


Figure  11  .--Projectile  Bead  Face.  After  Che  flashlamp 
waa  potted  into  the  projeceilc  head,  the  projectile  head  was  completed 
by  installing  the  inpactor  and  the  RIKA  contactor  ring  with  its  brass 
contactor  strips. 
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Figure  12  .—Projectile  Head  Cross  Section.  This  figure  shows 
how  the  flashlamp  sits  under  the  impactor  and  how  electrical  connection 
is  nade  between  the  flashlamp  and  the  brass  contactor  strips. 
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H.  ELECTRONIC  FUNCTIONS  AND  EQUIPMENT 


Electronics  were  used  to  perform  a  number  of  tasks  in  conducting 
these  experiments.  These  tasks  Included  the  following: 

1.  Provide  high  voltage  power  and  a  trigger  signal  to 
the  flashlasp. 

2.  Provide  timing  signals  to  synchronize  the  flashlamp 
and  oscilloscope  traces  with  the  experiments. 

3.  Record  the  signal  from  the  photomultiplier  in  the 
streak  camera. 

4.  Measure  projectile  velocity  and  streak  camera  mirror 
speed. 

Figures  13  and  14  show  the  discharge  circuits  used  to  drive  the 
flashlamps  for  the  reflection  experiments  and  the  transmission  experiments, 
respectively. 

Figure  15  shows  the  circuit  used  to  charge  the  trigger  pin  which 
detected  arrival  of  the  projectile  at  the  target.  Shorting  the  center 
conductor  of  the  trigger  pin  to  ground  produced  a  negative  voltage  transient 
which  was  uaed  to  trigger  the  oscilloscopes  and  the  flashlamp. 

The  photomultiplier  which  was  added  to  the  streak  camera  was  an 
RCA  type  4831.  Since  the  photomultiplier  is  Inherently  much  more  sensitive 
than  the  photographic  film,  It  was  necessary  to  limit  its  sensitivity.  This 
was  accomplished  by  operating  at  a  low  cathode  potential  of  -  650  V,  and  by 
placing  a  neutral  filter  of  density  1  In  the  light  path  to  the  photomultiplier. 
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Figure  13  — Flashlamp  Discharge  Circuit.  This  circuit  was  used 

to  drive  the  flashlaap  in  reflection  experiaents.  The  300  kfi  resistance 
was  located  in  the  Cordin  control  unit.  The  other  coaponents  were  in  a 
box  near  the  flashlaap. 


Figure  14  — Flashiemp-in-projectile  Discharge  Circuit 


Figure  15  —Trigger  Pulse  Circuit. 
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Projectile  velocity  was  determined  by  measuring  the  time  required 
for  the  projectile  to  traverse  the  accurately  measured  Interval  between 
two  small  brass  pins.  There  were  four  such  pins  forming  three  measured 
intervals.  Associated  with  each  measured  pin  was  a  similar  grounded  pin 
which  made  contact  with  the  projectile  just  before  the  measured  pin  to 
insure  that  the  projectile  was  a  good  electrical  ground  when  it  hit  the 
measure  pin.  Tektronix  DC  505-A  universal  counter  timers  were  used  to 
measure  the  time  intervals. 

The  speed  of  the  rotating  mirror  in  the  streak  camera  was  determined 
by  measuring  the  frequency  of  a  signal  from  a  pickup  coil  mounted  near  the 
magnetized  rotor  shaft  of  the  streak  camera  mirror. 


5.3  RESULTS  ANO  DISCUSSION 

Reflection  Experiments 

Pertinent  data  for  each  of  the  experiments  are  summarized  in  Table  I. 
The  600  lines/mm  grating  was  used  in  all  the  reflection  experiments.  It 
was  set  to  record  data  only  in  the  UV,  as  shown  in  Figs. 16-18.  No  reflectio 
measurements  were  made  in  the  visible. 

There  were  four  reflection  experiments  with  peak  pressures  of  99, 

123,  83  and  127  kbars.  Fig.  16  shows  the  streak  camera  and  photomultiplier 
(PM)  records  for  the  123  kbar  experiment.  Before  the  first  shock  enters 
the  CS2  the  records  show  light  transmitted  through  a  bandpass  from  2670  to  2' 
between  two  absorption  regions  in  the  CS^  spectrum.  At  the  time  of  impact 
a  slight  decrease  in  light  intensity.  For  the  first  180  ns  after  impact, 
a  123  kbar  shock  was  traveling  through  the  front  face  of  the  cell.  When 
it  reached  the  CS2-filled  cavity,  all  light  through  the  cell  abruptly  cut 
off  for  a  period  of  about  30  or  40  ns.  Light  then  reappeared  and  increased  ^ 
in  intensity  as  the  pressure  in  the  cell  rang  up.  The  light  which  reappeared 
had  a  broad  band  spectrum.  After  reaching  its  peak  intensity  at  C  in  the  fi 5 
the  light  intensity  recorded  by  the  photomultiplier  decayed  in  the  interval  Zl 
according  to  the  relation 

*(t)  *  ^xf0,73  +  0,27  expf28~0n^ 

while  the  sample  was  held  at  peak  pressure. 

Fig. 17  shows  the  streak  camera  and  PM  records  of  the  83  kbar  experiment 
The  same  features  are  evident,  except  for  the  decay,  as  were  seen  In  the  123 
kbar  experiment.  These  Included  pre-impact  transmission  between  2850  and 
2670  A,  a  slight  decrease  in  light  at  impact,  extinction  of  nearly  all  light 
as  the  first  shock  tranversed  the  CS2,  and  increasing  reflection  as  the 


pressure  In  the  sample  rang  up. 
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There  was  ambiguity  In  the  reflection  setup  in  that  light  reflected 
from  either  the  mirror  or  the  CSg  sample  could  reach  the  camera  assembly 
and  be  recorded.  The  127  kbar  experiment  shown  in  Fig.  18  was  an  attempt  to 
determine  where  the  reflection  was  occurring.  The  mirror  was  omitted  in  this 
experiment.  This  experiment  produced  a  record  that,  after  arrival  of 
the  shock  at  the  CSg  sample,  was  very  similar  to  the  123  kbar  experiment 
indicating  that  the  broad  band  light  was  reflected  from  the  shocked  CS£ 

layer,  not  the  mirror.  The  record  before  impact  shows  some  light  in  the 

Q 

2700  A  bandpass.  This  was  reflected  from  various  sapphire-vacuum  surfaces 

in  the  cell  and  impactor.  The  record  from  shot  number  79-057  was  very  weak  and 

is  not  shown,  but  it  displayed  the  same  general  features  as  the  other  reflection 
experiments. 

Reflection  Experiment  Summary 

The  observations  made  in  these  reflection  experiments  can  be 
summarized  as  follows.  The  sample  absorbed  nearly  all  light  for  a 
period  which  correlates  well  with  the  time  the  first  shock  was  traversing 
the  sample  cell.  Following  this,  the  sample  became  progressively  more 
reflective  while  the  pressure  In  the  cell  was  ringing  up.  In  the  higher 
pressure  experiments,  123  and  127  kbars,  the  reflectivity  reached  a  peak 
then  decayed  while  the  sample  was  held  at  high  pressure.  This  decay  in 
the  reflectivity  did  not  occur  In  the  83  kbar  experiment.  The  values  of 
reflectivity  were  estimated  to  be  in  the  range  of  0.2  to  l.o.  Table  II 
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Tab la  II  ~Shockad  Carbon  Disulfida  Raflactivicy  Data 


WiTHi*""  *a  Carbon  Disulfida  Raflactivicy 

Shoe  Praaaura  1  Max. ,  Eq.  (3.9)  Min.«  Eq.  (3.10) 

P 


Transmission  Experiments 

Transmission  experiments  are  summarized  In  Table  II.  All  were 
done  with  a  flashlamp  in  the  projectile.  The  range  of  wavelengths 
was  varied  to  cover  both  the  near  UV  and  the  visible. 

The  first  three  transmission  experiments  had  peak  pressures  of  55, 

60  and  57  kbars,  respectively,  and  used  progressively  thinner  cells,  145,  10, 
and  0.8  ym,  in  attempts  to  obtain  more  light  transmission.  The  first  cell 
was  thick  enough  that  the  transit  times  for  the  first  few  shocks  could 
be  resolved  by  the  streak  camera.  Fig. 19  shows  the  resulting  streak 
camera  record  in  which  both  absorption  regions  are  seen  to  have  expanded 
toward  the  red  as  the  first  shock  traversed  the  cell.  The  blue  edge  of 
the  central  absorption  band  did  not  change.  As  the  second  shock  started 
through  the  cell,  the  absorption  regions  expanded  again,  and  the  short 
wavelength  absorption  region  merged  with  the  longer  wavelength  region. 

At  the  time  of  the  third  shock,  all  light  was  cutoff  to  the  edge  of  the 
film  at  4070  A. 

In  the  next  two  experiments  the  cells  were  too  thin  for  the 
streak  camera  to  resolve  the  shock  transit  times,  so  the  light  was  cut 
off  abruptly,  except  for  a  slight  amount  at  the  long  wavelength  edge  of 
the  film  In  the  experiment  with  the  0.8  ym  thick  cell.  Records  for  this 
shot  are  shown  In  Fig.  20.  The  Increase  In  light  Intensity  at  Impact  results 
from  the  elimination  of  two  reflecting  surfaces  as  the  impactor  and  target 
sapphire  make  contact.  Light  cutoff  occurs  first  In  the  absorption 

O 

band  at  3200  A.  Ring-up  time,  given  In  Table  II,  Is  2  nsecs,  and  cutoff 

O 

In  the  lower  pass  band  (^3050  A)  occurs  approximately  one  nsec  after  shock 
arrival.  The  slit  Image  on  the  film  Is  30  nsecs  wide,  so  there  may  be 
considerable  error  In  the  apparent  9  nsec  delay  between  cutoff  at  3200  A 
and  3050  A. 
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The  records  from  shot  number  81-017  are  shown  In  Fig.  21  Here  the  150  11ne/H 
grating  has  been  used  to  Include  the  visible.  Cell  thickness  Is  again 
very  small,  approximately  0.6  urn,  and  the  same  delay  as  In  the  previous 
shot  exists  between  shock  arrival  at  the  CS,,  as  measured  by  extinction  to  3300  A 
and  cutoff  at  3600  A  In  the  pre-shock  transmission  region.  The  combined  absorption 
regions  expand  abruptly  to  about  4100  A  then  recede  slightly  to  3940  X,  during 
the  next  200  nsec.  The  peak  pressure  in  this  experiment  was  58  kbars. 


The  final  experiment  In  this  series  was  also  with  a  thin  cell,  but 
with  Impact  pressure  »  121  kbars;  spectral  range  Is  again  near  UV  and 
visible.  Records  are  shown  In  Fig. 22.  The  photomultiplier  record  shows 
that  Integrated  light  Intensity  dropped  to  0.08  of  Its  preshock  value 
within  20  nsec  as  the  shock  reached  the  carbon  disulfide.  The  Integrated 
Intensity  then  Increased  slightly  for  about  30  nsec  to  about  0.09  of  Its 


preshock  value.  This  small  remaining  Intensity  decayed  to  an  undetectable 
level  over  the  next  500  nsec.  On  the  streak  camera  record,  the  0.08 
Intensity  level  caused  only  a  very  slight  darkening  In  the  center  of  the 
film.  Most  of  the  light  signal  was  abruptly  cutoff.  Some  densitometer 
scans  In  the  time  direction  are  shown  In  Fig.  23.  Rise  times  range  from 
about  30  nsec  to  about  90  nsec.  Taking  the  shortest  values  to  result  from 
streak  camera  resolution  and  densitometer  response,  there  remains  a  fast 
relaxation  time  of  the  order  of  50  nsecs  to  opacity  in  the'cfhtral  region  In 
addition  to  the  500  nsec  time  for  the  Integrated  Intensity  recorded  by  the 
photomultiplier. 
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DISCUSSION 


The  broad  features  of  the  effect  of  shock  waves  on  transmission  In  the 
liquid  at  an  Impact  pressure  of  58  kbars  are  shown  In  Fig. 21.  Ring-up  time 
to  near  Impact  pressure  In  this  experiment  was  about  2  nsec,  much  less  than 
the  resolution  of  the  recording  system.  Cutoff  of  transmitted  light  to 
4100  A  was  completed,  l.e.  Intensity  was  reduced  by  at  least  a  factor  of 

O 

ten,  within  the  time  resolution  of  the  system.  A  density  scan  at  4750  A 

shows  no  discontinuity  when  the  shock  enters  the  CSg.  The  abrupt  drop  In 

0 

Intensity  at  wavelengths  greater  than  6000  A  Is  partially.  If  not  entirely, 
due  to  the  cutoff  of  UV  contributions  appearing  In  the  second  order. 

Thickness  of  the  cell  In  shot  number  81-010,  Fig. 19  was  145  um,  which 
was  large  enough  that  the  effects  of  successive  reverberations  could  be 
recorded.  Three  reverberations  and  the  corresponding  edge  shifts  of  Region  A 
are  required  to  extinguish  the  upper  transmission  region;  the  lower  one  requires 
only  two.  The  steps  are  blurred  because  of  the  30  nsec  window  and  limited  spectral 
resolution  system.  With  wavelengths  of  the  successive  steps  chosen  as  shown  In  the 
figure,  the  dependence  of  shift  frequency  on  calculated  pressure  and  temper¬ 
ature  Is  given  In  Fig.  22.  The  upper  edge  of  the  film  in  Fig. 19  Is  at  the 
upper  wavelength  cutoff  limit  shown  In  Fig.  21.  The  linear  dependence  on 
pressure  Is  remarkable  and  unexpected.  Shown  for  comparison  are  a  static 
pressure  measurement  at  room  temperature  of  CS2  dissolved  In  ethanol  made  by 
J.  Schnur  and  colleagues,^  (Fig.  25),  and  a  measurement  at  123°C  and  near  ambient 
pressure  made  In  this  laboratory  (Fig.  24  ).  The  last  mentioned  point  Involves 
a  subjective  Judgement,  as  can  be  seen  by  examining  the  displacement  between 
the  solid  curves  of  Fig.  24.  The  eye.  In  scanning  a  film,  picks  out  a  point 
at  a  particular  wavelength  and  reports  It  as  the  band  edge.  If  this  point 
corresponds  to  very  small  absorbance  or  very  large  absorbance,  the  edge 
displacement  due  to  the  Increase  In  temperature  will  appear  to  be  quite  large. 

If  It  corresponds  to  Intermediate  absorbance,  displacement  Is  small.  In  the 


present  case  two  Independent  observers  determined  band  edge  positions  by 
eye  directly  from  film  records  at  293°K  and  396°K  and  compared  them  with 
densitometer  traces.  Results  were  quite  similar  for  the  two  observers  and, 
when  corrected  for  lag  In  the  densitometer  recorder,  lay  In  the  region  of 
intermediate  absorbance. 

Neither  static  pressure  nor  temperature  alone  Is  adequate  to  explain  the  shifts 
produced  by  the  shock  waves,  nor  Is  their  sum.  It  Is  possible  that  a  combination  of 
pressure  and  temperature  would  produce  a  much  larger  shift  than  either  by  Itself,  or 
that  Isolation  of  CS^  molecules  In  solution  removes  some  cooperative  effect 
that  Is  significant.  It  Is  Important  to  note.  In  this  regard,  that  Raman 
measurements  of  the  symmetric  stretch  vibration  In  CSg  are  not  significantly 
affected  by  the  application  of  pressure,  even  above  12.6  kbars,  where  It 
transforms  to  a  solid;1  This  supports  the  view  that  there  Is  little  molecular 
Interaction  In  condensed  CSg  and  that  the  static  pressure  measurements  In 
solution  can  legitimately  be  compared  with  pure  CS2  measurements.  It  Is 
obviously  Important  to  explore  the  effects  of  static  pressure  more  thoroughly. 

The  third  step  In  the  motion  of  the  upper  band  edge  In  Fig.  19  extends 
the  edge  to  4070  A  or  beyond.  Comparing  this  with  Fig.  21  we  see  that  this 
Is  near  the  ultimate  limit  of  the  edge  for  pressures  In  the  neighborhood  of 
55-60  kbars.  The  third  reverberation  In  Region  A  of  Fig. 19  which  produces 
this  last  edge  motion.  Is  seen  from  Table  II  to  generate  a  temperature  of 
552°K  (shot  number  81-010).  This  Is  remarkably  close  to  the  temperature  of 
580 °K  at  which  Kleman  found  the  absorption  of  R  and  S  systems  In  the  gas  phase 
- —to-be-comparable  to  that  In  the  T  and  V  systems  (Fig.  5.1-1).  The  4300  A  upper 
limit  of  the  R  system  lies  beyond,  but  Is  comparable  to  the  4100  A  limit  shown 
In  Fig.  24.  These  near  coincidences  suggest  that  temperature  may  Indeed  play  a 
significant.  If  not  exclusive,  role  In  shifting  the  band  edges. 


If  temperature  Is  significant,  then  It  should  be  expected  that 
cooling  of  the  CS2  In  shot  number  81-017  (Fig.  21  )  might  play  a  role  In 
regression  of  the  band  edge  after  the  Initial  cutoff  to  4100  A.  Cell 
thickness  Is  only  0.6  um.  The  thermal  diffusion  coefficient  for  sapphire  Is 
a  hundred  times  greater  than  that  of  CS2,  so  the  sapphire  acts  approximately 
like  a  constant  temperature  sink  Insofar  as  the  CS2  layer  Is  concerned. 

Using  handbook  values  for  thermal  conductivity  and  specific  heat  of  CS2*  and 
correcting  for  density,  the  1/e- thing  time  for  the  CS2-sapph1re  temperature 
difference  Is  315  nsec.  The  band  edge  In  Fig.  21  after  Initial  cutoff  Is 
quite  fuzzy,  but  it  looks  as  If  the  1/e-thlng  time  might  be  closer  to  100  nsec. 
This  Is  not  unreasonable  agreement,  since  diffusion  may  be  considerably  enhanced 
by  pressure  and  temperature.  However,  regression  of  the  band  edge  seems  to  stop 
at  about  3940  A.  But,  cooling  continues  and  even  with  the  1/e-thlng  time 
of  315  nsec,  temperature  In  the  C$2  should  drop  to  about  325*K  after  1000  nsec. 
The  record  of  Fig.  21  shows  nothing  like  this. 

Cutoff  of  the  lower  passband  In  Fig.  19  gives  less  Information  than  the  upper 
The  apparent  Initial  band  edge  at  2660  A  is  Influenced  by  rapidly  diminishing 
sensitivity  of  the  film  with  decreasing  \,  by  diminishing  reflectivity  of 
the  grating,  and  by  decreasing  sensitivity  of  the  camera  as  film  edges  are 
approached.  The  "true"  band  edge  may  be  considerably  below  2660  A. 

The  most  we  can  say  then  Is  that  the  first,  5  kbar,  shock  shifts  the  band 
edge  toward  the  red  by  at  least  150  A  and  the  second  shock  closes  the  pass 
band  entirely.  If  edge  shift  Is  proportional  to  pressure  for  this  lower  band, 
as  for  the  upper,  the  amount  of  pressure  required  to  dose  the  2660-2900  A 
band  gap  does  not  exceed  8  kbars.  Referring  to  Table  I,  we  see  that  this  is 
Just  the  pressure  of  the  first  shock  In  shot  number  80-010,  so  the  occurrence 
of  the  cutoff  region  before  reflection  Is  fn  accord  with  the  observations  on 


transmission  shot  number  81-010.  Since  first  shock  pressure  in  the  other 
reflection  shots  is  greater  than  8  kbars,  the  cutoff  preceding  reflection 
Is  consistent  with  the  transmission  experiment  In  all  cases. 

Of  the  four  reflection  shots  in  Table  I,  three  are  shown  In  Figs. 16-18 
the  fourth,  shot  number  79-057,  produced  only  a  very  faint  record  because  of 
a  filter  in  the  light  path  and  Is  not  shown.  All  were  in  the  UV. 

Reflection  coefficients  for  the  experiments  shown  In  Figs. 16  and  17 
were  defined  In  terms  of  the  ratio  PM(E)/PM(0)  •  R,  where  PM(E)  Is  photo¬ 
multiplier  output  in  the  experiment,  e.g.,  the  lower  parts  of  Figs. 16  and  17 
and  PM(0)  Is  photomultiplier  output  measured  In  a  preshot  test  with  empty 
cell.  There  Is  some  variation  In  flashlamp  output  from  shot  to  shot; 
this  produced  a  correspondi ng  uncertainty  in  R.  There  are  six  reflecting 
Interfaces  in  the  experiment,  and  these  must  be  considered  in  calculating  R. 

Transmission  and  reflection  coefficients  for  normal  incidence  at  the 
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Interface  between  two  refracting  media  are 
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Light  signals  In  each  section  of  the  assembly  can  be  divided 
signals,  IL,  and  right-travelling  signals,  IR,  as  in  Fig.  25. 
internal  sections  are 
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into  left-travelling 
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At  the  right-hand  boundary  the  reflected  light  Is  lost  because  of  the  bevel 
on  the  back  sapphire.  Then 
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Define  Rj  *  XR  j/IL  th«n*  E<»s-  and  <3)  can  be  "ritten 


(1  ~  T)  +  (2T  -  1 )Rj 
1  -  (1  -  T)Rj 
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Since  R1  *  0.85  Is  the  reflectivity  of  the  mirror  at  the  far  left,  Eqs.  (1) 
and  (5)  can  be  combined  to  yield  a  value  of  R  for  each  region. 

For  the  empty  cell,  denoted  by  superscript  "e",  there  are  only 
sapphire-air  boundaries,  so  T  »  0.924.  Then  R|  *  0.857.  In  the  reflection 
experiment  the  cell  Is  filled  with  CSg.  Let  r c  be  the  reflectivity  of  the 
CS2  sample,  and  let  IR>4  ■  0.  Ift  =  IR  g  for  the  experiment,  and  Ip  »  IR^5 
for  the  empty  cell.  Then  r|  ■  r#  and  IR  5  Is  the  light  reflected  from  the 
CS2-sapph1re  Interface  or  from  the  shock  front.  Thus 
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There  Is  some  uncertainty  In  the  determination  of  r|  because  the  angle 
of  Incidence  was  not  exactly  normal  and  some  of  the  reflected  light  did  not 
fall  on  the  entrance  si  It  of  the  spectrograph.  If  no  reflections  are  lost, 

Eq.  (6)  applies;  If  all  are  lost. 
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Eqs.  (6)  and  (7)  bound  the  CS2  reflectivity.  Values  obtained  for  shots 
80-009  and  80-010  are  given  In  Table  II.  These  values  of  reflectivity  are 
extraordinarily  large  and  are  more  nearly  characteristic  of  a  poor  conductor 
than  of  a  dielectric. 

The  phenomena  of  light  extinction  In  the  transmission  experiments  and 
of  broad  band  reflection  In  the  reflection  experiments  cannot  be  Independent. 
Reflection  measurements  did  not  extend  Into  the  visible,  but  It  Is  reasonable 

O 

to  suppose  that,  at  pressures  up  to  60  kbars,  reflection  would  cease  at  4100  A 
just  as  extinction  does.  Since  reflection  results  from  scattering,  and  extinction 
results  from  scattering  plus  absorption.  It  Is  evident  that  absorption  Is  small 
compared  to  scattering.  The  magnitude  of  the  reflection  Is  such  that  It  could 
not  be  generated  by  Incoherent  scattering;  the  solid  angle  subtended  by  the 
detector  Is  too  small.  In  order  to  be  an  effective  reflector  at  a  given 
frequency,  the  scattering  must  be  temporally  coherent,  l.e.,  the  phase  lag 
between  absorption  and  emission  must  be  constant,  and  spatially  coherent  over 
an  area  of  diameter  d  such  that  X/d  «  1.  The  thickness  of  the  scattering 
layer  must  also  be  a  small  fraction  of  a  wavelength  In  order  that  signals  from 
front  and  back  of  the  layer  be  additive.  If  we  take  scattering  cross-section 
per  unit  area  to  be  unity  and  thickness  of  the  scattering  layer  to  be 
a  <  0.1  ym,  then  scattering  cross  section  per  molecule  Is 

o$  ■  1/na 

23 

where  n  •  number  of  molecules/cc  •  0.17  x  10  Then 

cs  >  10“17cma 

This  Is  at  least  an  order  of  magnitude  greater  than  the  absorption  cross-section 
of  CSg  at  the  center  of  region  A  estimated  from  the  film  before  shock  arrival 
In  Fig.  21. 
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The  existence  of  a  thin  scattering  layer  of  large  cross  section  Is  In 
accord  with  the  observation  that  light  transmitted  through  a  cell  0.6  urn  thick 
Is  cut  off  as  completely  by  the  shock  as  In  a  thicker  layer.  This  suggests 
that  the  scattering  layer  Is  associated  wl th— perhaps  takes  place  In— the 
shock  front,  but  this  Is  speculation.  Two  possible  mechanisms  which  cou\d 
account  for  coherent  scattering  of  the  observed  magnitude  are  formation  of 
conducting  complexes  and  resonance  fluorescence,  possibly  from  excited  states. 
Further  Investigations  will  be  necessary  before  we  can  develop  a  more  detailed 
understanding. 

The  spectrograms  shown  In  Figs.  16  and  22  contain  some  evidence  for  the 

supposed  decomposition  reaction  which  Initiated  this  work.  The  reflection 

experiment  at  127  kbars  Impact  pressure  given  in  Fig. 16  shows  a  decay  of 

reflected  intensity  with  a  time  constant  of  the  same  magnitude  as  those 

13 

reported  by  Sheffield  at  the  same  temperature.  The  extended  cutoff 

O 

shown  In  Fig.  23  for  the  4800  A  region  has  a  shorter  relaxation  time  and 
a  somewhat  smaller  pressure,  121  kbars.  Both  these  phenomena  may  be  mani¬ 
festations  of  decomposition.  The  transmission  record  of  Fig. 21  Is  in 

14  IS 

accord  with  the  observation  by  Yakushev,  et  al.,  ’  of  yellowing  of 
visible  light  as  decomposition  approached. 

Much  of  the  foregoing  discussion  Involves  pressures  and  temperatures 
recorded  in  Tables  II  and  III.  The  nature  of  these  experiments  requires  that 
these  temperatures  be  calculated  from  an  equation  of  state.  The  one  used 
here  was  developed  by  Sheffield.**  It  has  some  limitations,  as  do  all 
-such- creations,  but  pressures  should  be  quite  accurate.  Temperature  errors 
have  not  been  evaluated. 
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Fig.  16  Spectrogram  and  photomultiplier  record  for  reflection  experiment. 
Shot  number  80-009,  ■  123  kbars. 
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Fig.  17  Spectrogram  and  photomultiplier  record  for  reflection  experiment 
Shot  number  80-010,  P<  *  83  kbars. 
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Fig  18  Reflection  experiment  with  no  mirror  In  the  projectile. 
Shot  number  80-013,  P1  «  127  kbars. 
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Shot  nuaber  81-014,  Pj  >  57  kb«rs. 
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Ffg.  25  continued 


Fig.  26  Schematic  for  calculation  of  light  Intensities  In 
. cell  and  Impac-tor. 


VI.  DISCUSSION 


There  are  basically  two  questions  one  wishes  to  address  In  the  kinds 
of  shock  wave  studies  described  here: 

1.  Do  shock  wave  techniques  provide  a  useful  tool  for 
studying  the  kinetics  of  chemical  reactions? 

2.  Are  shock-induced  chemical  reactions  explicable  In 
terms  of  conventional  chemical  kinetics,  or  are 
there  hitherto  unexplored  areas  of  chemistry  Involved? 

The  answer  to  the  first  Is  certainly  yes,  as  demonstrated  by  Sheffield's 
work  In  Section  IV  and  to  a  lesser  extent  by  Ogllvle's  work  in  Section  V. 
These,  however,  are  but  demonstrations  of  feasibility.  The  breadth  of 
application  and  the  utility  of  the  Information  obtained  remains  for  further 
investigation. 

Answers  to  the  second  are  not  yet  forth  coming.  It  has  several 
subdivisions:  are  the  high  pressures  and  temperatures  produced  by  the  shock 
sufficient  to  produce  the  observed  effects?  Are  P  and  T  along  with  the 
anticipated  creation  of  defects  by  the  shear  accompanying  the  shock  front 
enough?  Are  there  non-equilibrium  energy  transfers  from  shock  to  molecule 
which  Increase  Its  reactivity?  Is  there  shear  deformation  of  molecules  which 
contribute  to  the  observed  effects?  Answers  to  these  questions  will  not  be 
easily  obtained.  Spectroscopy  Is  clearly  a  useful  tool.  As  long  as  It 
produces  unexpected  results.  Its  potential  cannot  be  fully  evaluated.  The 
spectra  are  unusual;  they  are  not  the  kind  described  In  textbooks,  but  they 
clearly  contain  Information. 
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Cumulative  evidence  suggests  that  the  shock  front  is  locally  very 
thin— probably  a  few  atomic  spaclngs  or  less.  This  certainly  makes  It 
possible  to  think  of  Impulsive  energy  transfer  to  vibrational  states  as  a 
method  of  non-equilibrium  excitation.  The  most  clear-cut  piece  of  data 

o 

that  relates  to  this  Is  obtained  from  shot  number  81-010,  where  a  120-150  A 
shift  of  the  red  absorption  edge  Is  obtained  from  the  Initial  5  kbar  shock. 

This  appears  to  be  out  of  agreement  with  static  pressure  and  temperature 
measurments  by  a  wide  marglne.  But  the  static  data  are  equivocal  and 
additional  experimental  support  Is  required.  A  new  result  obtained  since 
Ogilvie's  work  was  completed  Is  surprising.  Rapid  adiabatic  compression  of 
the  CS^  sample  by  a  high  amplitude  compressive  ramp  shows,  within  experimental 
error,  the  same  rate  of  shift,  approximately  25  S/kbar,  as  the  shock.  In 
entirely  different  experiments,  ultronsonlc  absorption,  it  has  been  found 
that  liquid  CS2  has  relaxation  times  of  5-50  nsec  associated  with  vibrational 
states.  Rise  time  In  the  ramp  experiment  was  about  *s  ki lobar  per  nsec,  which 
suggests  a  relation  among  the  experiments. 

Our  research  Is  continuing  along  the  lines  described  In  Section  V,  with 
support  and  cooperation  of  other  laboratories,  particularly  J.  Schnur  and 
colleagues  at  NRL,  and  Kaufman  and  Koskl  at  John  Hopkins.  Efforts  are 
underway  to  do  Raman  spectroscopy  In  Shaner's  group  at  Los  Alamos,  and  Graham 
and  colleagues  at  Sandla,  Albuquerque,  are  preparing  to  do  optical  spectroscopy. 
With  this  concentration  of  effort,  we  may  reasonably  expect  some  definitive 
answers  to  Question  2  In  the  fairly  near  future.  With  these  will  come 
further  clarification  of  the  role  of  spectroscopy  In  shock  wave  chemistry. 
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Some  Notes  on  the  Transition  from  Deflagration  to  Detonation 


Insufficient  attention  has  been  paid  to  the  logical  framework  provided  by  the 
Chapman-Jouguet  theory  of  detonation  and  deflagration  for  understanding  of  the 
deflagration  to  detonation  transition.  A  brief  outline  of  its  application  to  the 
problem  is  presented  here. 

The  C-J  Is  a  steady  state  theory,  applying  strictly  only  to  those  hydrodynamic 
and  reaction  processes  which  are  plane,  one-dimensional,  and  unchanging  In  time  when 
viewed  from  the  correct  moving  coordinate  system.  But  we  may  expect  that  the  theories 
will  apply  approximately  to  unsteady  processes,  provided  significant  changes  do  not 
occur  during  times  the  order  of  [reaction  zone  thickness/propagation  velocity]. 
Moreover,  the  C-0  theories  are  independent  of  reaction  time  and  reaction  zone  thickness 
provided  only  that  the  flow  Into  and  out  from  the  reaction  zone  Is  one-dimensional 
and  steady.  This  means  for  example,  that  combustion  Involving  propellant  fracture 
falls  within  the  framework  of  the  theory  as  long  as  flow  Into  and  out  from  the  fracture 
combustion  region  is  one-dimensional  and  steady.  Deviations  from  steady  flow  can  be 
treated  as  perturbations.  If  they  are  not  too  large.  With  this  In  mind  one  can 


proceed  as  follows. 

AB  In  Fig.  1  Is  the  unreacted  R-H  curve 
of  the  propellant,  centered  at  (p0,VQ).  CDEF 
Is  the  locus  of  states  of  the  fully-reacted 
products  reached  from  the  Initial  state  p0,VQ 
by  detonation  or  combustion  In  a  plane  reaction 
front.  The  line  00  Is  parallel  to  the  V-axIs; 


216 


V 


t 


-2- 


OE  Is  parallel  to  the  p-axls.  States  on  CO  are  reached  from  0  by  combustion;  those 
on  EF  are  reached  by  detonation.  The  speed  of  propagation  of  the  plane  reaction 
front  with  respect  to  the  material  ahead  of  It  Is 


U2 


,2  prpo 


o  '1 


On  CD,  p^  <  PQ,  VQ  <  Vj.  On  EF,  p^  >  pQ,  VQ  >  V^.  In  order  for  deflagration  to 
progress,  pg  must  be  greater  than  p1 .  If  pQ  ■  p^,  U  ■  0  and  the  reaction  zone  does 
not  move.  Since  unbumed  material  must  be  consumed  continually  In  order  for  the 
reaction  to  exist,  such  a  constant  pressure  deflagration  cannot  exist,  but  It 
provides  a  bound  to  possible  deflagrations.  In  order  for  the  deflagration  to 
accelerate  Pj-PQ  must  get  larger,  VQ-V^  smaller,  or  both.  The  reacted  products  are 
normally  gaseous;  If  the  unbumed  material  Is  solid,  the  two  curves  AOB  and  CDF 
can  be  expected  to  converge  at  high  pressures  and  diverge  at  low  pressures.  It  Is 
not  possible,  then,  for  U  to  become  very  large  If  (P0>VQ)  remains  fixed  and  Pi 
decreases.  U  may.  In  fact,  decrease  as  pj  decreases  because  of  the  rapid  Increase 

of  V,  at  low  pressures. 

*  —  \ 

It  follows,  then,  that  If  a  deflagration  Is  to  accelerate,  the  point  (P0*v0) 
must  migrate  upward,  and  It  must  do  so  more  rapidly  than  (Pj,Vj).  That  Is,  the  mean 
pressure  In  the  deflagrating  system  must  Increase  and  pQ  must  remain  greater  than 
P].  I*  P0-Pj  remains  constant  as  pQ  Increases,  the  wave  will  continue  to  accelerate 
because  -V0  will  grow  smaller.  This  requirement  that  p  ■  (P0*P-j)/2  Increase 
Implies  that  the  region  of  deflagration  must  somehow  be  protected  from  the  Influence 
of  pressure-free  boundaries.  Otherwise  the  pressure-relief  waves  being  continually 
produced  by  these  boundaries  will  overtake  the  region  of  deflagration  and  will 
reduce  p,  thus  preventing  the  acceleration. 

This  Isolation  from  the  effects  of  pressure-free  boundaries  can  be  accomplished 
In  two  ways:  by  confining  the  explosive  In  a  fixed  volume  so  that  the  released  energy 
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of  combustion  can  drlvt  tho  pressure  up,  or  by  causing  the  boundaries  to  be 
far-removed  fro*  the  region  of  deflagration;  l.e.,  by  using  a  very  large  piece  of 
explosive.  "Large"  In  this  context  will  depend  on  propagation  velocities  and 
reaction  rates. 

The  conditions  which  must  be  created  In  order  for  the  deflagration  to  transform 
to  a  detonation  can  now  be  described  more  precisely.  If,  by  judicious  adjustment 
of  the  mean  pressure  and  the  pressure  difference,  we  can  arrive  at  the  situation 
Indicated  In  Fig.  2,  where  the  straight  line  connecting  the  "Initial  state",  (pn*Va), 
and  the  undisturbed  or  reference  state, 

Intersects  or  Is  tanget  to  the  reacted 
Hugonlot,  we  will  have  arrived  at  a  condition 
In  which  a  non-reactlve  shock  wave  can  produce 
the  transition  from  p^.V^  to  po,VQ.  This, 
combined  with  subsequent  deflagration  to 
p<l ,Vj ,  constitutes  a  detonation  wave. 

These  Ideas  are  elementary,  and  for  that  reason  may  be  useful.  It  Is  probably 
reasonable  to  assume  that.  In  one-dimensional  flow,  we  can  always  find  a  plane  ahead 
of  the  reaction  zone  where  no  reaction  has  occurred  and  a  plane  behind  where  all 
reaction  Is  complete,  or  nearly  complete.  If  we  make  the  separations  between  these 
planes,  l.e.,  the  reaction  zone,  as  small  as  possible,  we  should  be  able  to  direct 
our  attention  to  what  Is  happening  outside  the  reaction  zone  without  too  much 
concern  about  the  details  of  what  Is  happening  Inside.  To  be  sure,  one  of  the  things 
that  will  occur  Is  that  the  reaction  zone  thickness  will  change  as  the  deflagration 
accelerates,  but  that  effect  may  not  be  of  primary  significance. 

I  have  made  some  preliminary  attempts  to  map  the  course  of  an  accelerating 
reaction  In  the  p-u  plane  without  notable  success.  It  will  require  careful  thought 
and  probably  the  working  out  of  examples  with  various  boundary  conditions.  But  I 
think  It  will  be  worth  the  effort. 
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